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Kinematic Wave Models of Network Vehicular Traffic 

Abstract 



The kinematic wave theory, originally proposed by (jLighthill and Whitham 



Richard; 



^^^^^ 



19561 ) ■ has been a good candidate for studying vehicular traffic. In this 
dissertation, we study kinematic wave models of network traffic, which are expected 
to be theoretically rigorous, numerically reliable, and computationally efficient. 

For inhomogeneous links, we reformulate the Lighthill-Whitham-Richards model 
into a nonlinear resonant system. In addition to shock and rarefaction waves, standing 
(transition) waves appear in the ten basic w ave solutions. The solutions are co nsistent 



mm). 



with those by the supply-demand method fpaganzd Il995al : ILebacqud . 

For merging traffic, we examine existing supply-demand models and, particularly, 
distribution schemes. Further, we propose a new distribution scheme, which captures 
key merging characteristics and leads to a model that is computationally efficient and 
easy to calibrate. 

For diverging traffic, we propose an instantaneous kinematic wave model, con- 
sisting of nonlinear resonant systems. After studying the seven basic wave solutions, 
we show that this model is equivalent to a supply-demand model with modified defi- 
nitions of traffic demands. 

For traffic with mixed-type vehicles, we show the existence of contact waves. 
Using simulations by the developed Godunov method, we demonstrate that First-In- 
First-Out (FIFO) principle is observed in this model. 



11 



For network traffic fiow, we propose a multi-commodity kinematic wave (MCKW) 
model, in which we combine kinematic wave models of different network components 
and a commodity-based kinematic wave theory. We also propose an implementation 
of the MCKW simulation and carefully design the data structure for network topology, 
traffic characteristics, and simulation algorithms. The solutions are consistent with 
FIFO principle in the order of a time interval. 

For a road network with a single origin-destination (0/D) pair and two routes, 
we first demonstrate the formation of an equilibrium state and find multiple equilib- 
rium status for different route distributions. We then show the formation of periodic 
oscillations and discuss their structure and properties. 

Finally, we summarize our work and discuss future research directions. 
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Chapter 1 



Introduction 



1.1 Background 
1.1.1 Traffic congestion 

As the backbone of the intermodal transportation network in the United States, road 
networks - consisting of highways, arterial roads, surface streets, and other kinds of 
roadways - connect air, transit, rail, and port facilities and terminals. In particu- 
lar, highways carry 90 percent of passenger travel and 72 percent value of freight. 
Therefore, the performance of its road networks largely defines the mobility of the 
nation and affects economic and social activities in the United States. However, traffic 
conditions on road networks in many metropolitan areas are becoming increasingly 
congested: in 68 major American urban areas, the percentage of un-congested periods 
was two third s of the whole peak period in 1982, while the percentage drops to one 



third in 1997 ( Schrank and Lomax 



19991 ) 



In a road network, traffic congestion can be recurrent and non- recurrent. Recur- 
rent traffic congestion is generally caused by limited physical infrastructure, increasing 
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travel demand, rush hours, and toll booths. Non-recurrent congestion is associated 
with accidents, work zones, and weather. Traffic congestion increases travel delay 
and fuel consumption and adversely affects safety, mobility, productivity, the human 
and natural environment. As a result, increasing delay has seriously damaged the 
speed and reliability of road networks, which are vital to emerging industries, like 
warehousing and logistics. 

It has been a major goal of transportation scientists and engineers to alleviate 
traffic congestion. Many strategies have been pursued to achieve this goal. One 
method is to expand existing road facilities by adding new roads or lanes, rebuilding 
key network components, and enhancing the physical condition of roadways. Another 
method is to apply traffic management and operation technologies, including auto- 
matic highways, travel demand management, freeway management, incident manage- 
ment, emergency response management, weather response management, value pricing, 
arterial signal control, on-ramp metering, traveler information, changeable message 
signs, and so on. 

All these approaches are intended to either increase the capacity of road networks 
or reduce their load. Among them, the method of expanding existing road networks 
has seen limited use due to huge construction costs and the difficulty of addressing 
public and environmental concerns. Therefore, the idea of using existing infrastruc- 
ture more efficiently by building Intelligent Transportation Systems (ITS) is a more 
preferable option. 

1.1.2 The role of traffic models 

Whichever of the aforementioned strategies are taken to improve the mobility of 
a road network, their success relies on a better understanding of the properties of 
congestion and the overall performance of the road network. Widely available sensing. 
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information and communication technology has had a major influence in detecting 
traffic conditions in a timely manner. However, collecting meaningful traffic data is 
not an easy job. Moreover, many management and control decisions are grounded on 
an estimate of traffic conditions in the near or long-term future. Thus, it is essential 
to accurately estimate traffic conditions in a road network during certain time periods 
and understand the evolution pattern of traffic conditions, i.e., the traffic dynamics. 
This calls for the development of traffic models. 

A traffic model is a function which relates the movement of a vehicle to driver's 
behavior, vehicle type, network characteristics, weather conditions, traffic signals, 
guidance information, and interaction with other vehicles. The movement of a vehicle 
can be represented by its position at any time, i.e., its trajectory, from which its 
speed and acceleration rate can be obtained. Given all vehicles' trajectories, one 
can measure the performance of a road network, e.g., travel time, level of service, 
congestion level, etc. As shown in Figure II. H traffic models can considered as a 
theoretical substitution of a real traffic system. 

Besides the constraints in Figure traffic models are also subject to both 
space and time limitations, i.e., initial traffic conditions and boundary conditions. 
Moreover, all traffic models are subject to calibration and validation. In calibration, 
parameters in a traffic model are adjusted so that it acts as closely to the real traffic 
system as possible. In validation, the output of a traffic model is compared with 
observation of real traffic. 

The role of traffic models in transportation engineering is two- fold. First, they 
provide better understanding of traffic dynamics, in particular the formation and 
propagation of traffic congestion. Hence, traffic researchers can use traffic models 
to identify possible bottlenecks. Second, they can serve as a simulation platform, 
on which different strategies for improving mobility, as shown in Figure II. If can be 
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Figure 1.1: The definition and role of traffic models in traffic studies 



developed and evaluated. For example, in the plan for expanding a road network, 
traffic models can be used to simulate proposed expanded networks and help chose the 
most cost-effective strategy. For another example, they are also helpful in determining 
the best location of tolling booths, which are designed to divert traffic away from busy 
roads by charging fees. Finally, traffic models can be used to evaluate previously 
implemented strategies. 
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1.1.3 Traffic models and simulation packages 



In a vehicular traffic system, a number of trips - defined by their origin/destination, 
departure time/arrival time and travel route - interact on the road network and 
generate various dynamics and phenomena. To study these traffic phenomena and 
the corresponding applications, many traffic models and simulation packages have 
been proposed in literature. 

From a microscopic point of view, a stream of vehicular traffic is the combination 
of the movements of all vehicles. In a microscopic model of traffic, the movement 
of each vehicle and interactions between the vehicle pairs are studied. Three types 
of microscopic models have been suggested. O ne approach is the CM family of car- 



following models developed in the 1960's (e.g., fjGazis et al. 



196l|)). In these models. 



the movement of a vehicle is described by an ordinary differential equation in time. 



In another approach, known as coupled-map l attice models, a ve 



licle's movement is 



20001). Yet another 



defined by equations that are discrete in time (jChowdhurv et al. 
approach, cellular automata (CA) models are based on the framework of statistical 
physics, in which not only independent variable s (space and time) but als o descriptive 



variables (speed and acceleration) are discrete (jChowdhurv et al 



20001 ) 



From a macroscopic point of view, one approach to modeling traffic is to treat 
the traffic as it were a gas of interacting particles, in which each particle represents 
a vehicle. The resulting model is called a kinetic model of vehicular traffic. Another 
approach is to consider traffic flow as a compressible fluid (continuum). Such models 
are known as continuum models, or fluid-dynamical models. In continuum models, 
the basic characteristics are flow rate g, traffic density p, and travel speed f , which are 
all functions in time and space. The first-order continuum models are generally called 
kinematic wave models, in which traffic dynamics are regarded as a combination of 
"kinematic" waves in these quantities. 
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Based on the aforementioned theoretical models of vehicular tra ffic, many simula- 



ti on p ackages have been developed: for example, the PASSER series (ITexas Transportation Institute 



1991 



), TRANSYT dCourage and Wallace . 



1991 



1995D, DYNAS MART iTavakrishnan et all 

1"^ 



1998D, FRESIM (ISmith and Noel 



I INTEGRA^ 



mi), NETSIM 



1993), FREFLO (iPaynel. 



ION (IVan Aerde and the Transportat 



Federal Highwav Administration 



1979D. MITSIM flMassachusetts Institutt 
^ I99J 



19991) . TRANSIMS flTRANSIMSil . and PARAMICS (ICameron and Dimca 
In some simulation packages, however, the underlying traffic flow theories lack math- 
ematical rigor, particularly for traffic dynamics at highway junctions. Furthermore, 
numerical methods used in some simulations have not been well justified and may 
cause numerical instability. Another common drawback with many simulation pack- 
ages is their enormous computational cost. Therefore, many of these simulation 
packages have serious limitations in applications. That said, continuing efforts to 
develop theoretically rigor, numerically sound, and computationally effective models 
are still necessary, in particular for applications to dynamic traffic assignment and 



other advanced traffic engineering strategies. 



1.2 Continuum models 
1.2.1 Kinematic wave models 

In the kinematic wave models, traffic is viewed as a continuous media and charac- 
terized by traffic density (p), travel speed [v), and flow-rate (g). The movements 
of vehicles on a road network are considered as combinations of kinematic waves in 
either of these three quantities. 

The different types of kinematic waves are associated with different components 
of a traffic system. For example, on a homogeneous link with uniform conditions 
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in vehicles, drivers, weather, etc., one can observe two basic waves: decelerating 
shock waves, generally seen when lights turn red, and accelerating rarefaction waves, 
generally seen when lights turn green. In a traffic system with inhomogeneous links, 
merges, diverges, or different vehicle types, more complicated waves can be observed, 
such as standing transitional waves, contact waves, and periodic waves. 

In kinematic wave models, a hyperbolic conservation law is derived from traffic 
conservation: 



Pt + Qx 



0. 



1.11 



A fundamental assumption in kinematic wave models is that the flow-rate g is a 
function of the traffic density p; i.e., q = Q{p), which is called the fundamental 
diagram. Hence, v = V{p) = Q{p)/p. Generally, flow-rate is a concave function 
in density and retains its maximum, the capacity, at the critical density pc- When 
traffic density is higher than the critical density, it is in the over-critical region and in 
the under-critical region, otherwise. The fundamental diagram, or the speed-density 
rela tion, varies with link characteri stics, vehicle types, a nd so on. 



Lighthill and WhithamI (|l955br i and 



Richardsl (|l956| ) ffist proposed and analyzed 



kinematic waves on homogeneous links. The corresponding model is known as the 
LWR model and is described by a ffist-order, nonlinear PDE: 



pt + Q{p), 



0. 



;i.2) 



This model is an hyperbolic conservatio n law, whose Riemann pr oblems ^ can be 



analyzed using well developed tools (e.g.. 



Lax, 



1972; 



SmoUer . 



19831 ). Analysis shows 



that solutions to Equation ll.2[ where /(x, p) = pV^(p), have wave properties analogous 
^In a Riemann problem, the initial trafRc conditions are described by a Heaviside function or 
step function. 
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to those of water flow in channels (jLighthill and Whitham . 



1955ar i. In other words, 



the Riemann problem consists of either shock or rarefaction waves. 

Nume r ically , the LWR model can be solved with a first-order Godunov method 



( Godunov 



19591), in which a link is partitioned into a number of cells, a time dura- 



tion into a number of time steps, and traffic conditions in each cell at a time step are 
uniform. In a Godunov method, traffic conditions are updated according to the con- 
servation equation, Equation ll.21 i.e., during a time interval, the increasing number of 
vehicles in a cell are the difference between the in-fiow through its up stream boundar y 



and the out-fiow through its downstream boundary. Traditionally (jLeVeque . 



mm 



flows through cell boundaries are computed from wave solutions of Riemann problems. 



Daganzo 


1995a 


Lebacquf . 


199fS) 



19961 ). In this intuitive, engineering 



method, the flow through a boundary equals the minimum of the traffic demand of its 
upstream cell an d the traffic supply of th e downstream cel l. Here the traffic demand 



( Lebacque 



19961), called sending flow in (jDaganzo, 



1995ar i. of a cell is defined as its 



flow-rat e when tra f fic is u nder-critical or its c apacity when ove r-critical, and its traffic 



supply (|Lebacque 



1995a| ). is the capacity when 



19961 ). receiving flow in (jDaganzo , 
traffic is under-critical or the flow-rate when under-critical. Since one does not have 
to understand wave solutions of Riemann problems, this method has been widely 
applied in traffic studies. 

Many attempts have been made to model more complicated traffic systems under 



Lighthill and WhithamI ()l955bt ) 



the framework of the LWR model. For example, 
discussed the kinematic wave theory of traffic dynarn ics on an inhomogeneous link 



with lane-drops or curvatures. Lately, 



Lebacqud (|l996| ) gave a detailed analysis of the 



kinematic wave solutions of suc h a traffic system and summarized the solutions into 



the supply-demand method. In (jDaganzo , 



1995ar ) , the concepts of sending flow (traffic 
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supply) and receiving flow (traffic demand) was intuitively extended for simulating 
traffic on an inhomogeneous link. 

Another attempt is to model traffic dynamics on a highway network, where the 



focus is on highway 



unct ions, including merges, diverges, and other intersections. 



Holden and Risebrd (|l995[ ) studied kinematic waves initiated at highway junctions. 



by assuming the existence of an optimization problem at each junction and excluding 



route choice behavior. Without route choice, th is model sees limited app 



Kuhne and Michalopoulod (|l992r i 



ications 



in reality. In the network traffic flow model by 
on-ramps and off-ramps are considered as sources and sinks respectively. Although 
including the influence of ramps on mainline freeways, this models omits the other 
side: the influence of mainline freeways on ramps. Thus it also has limitations with- 
out giving a full picture of traffic dynamics. To cap ture overall traffic phenomena in 



a traffic system, (iDaganzo 



Lebacquf . 



199&} proposed some discrete models. 



in which traffic demands and supplies for cells around an intersection are defined the 
same as those for cells inside a link, and route choice behavior and certain optimiza- 
tion rules for fiows are incorporated in order to determine unique fiows through the 
intersection. In literature, there have been little progress in analyzing the kinematic 
waves initiated at a highway junction considering route choice behavior, mainly due 
to the difficulties in formulating them into a system of continuous partial differential 
equations as in the LWR model. 



for ad dressing different concerns in traffic system s. In 



19971 ). traffic on special lanes is i nvest igated. In (jPaganzo 



incorporated. 



s. In ( 


Dasranzo . 


1997 


Daeranzo et al.. 


(Daganzo. 


2002 


), driver behavior are 



Wong and Wond (120021) discuss differentiated vehicle 



teresting studies can b e found in (|Vaughan et al 



\2m 



Lebacquf . 



1984 : 



Newell 



t ypes. Other in- 



1993, 



1999; 



LeVeque 



2m 
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Compared to other traffic models, kinematic wave models have the following 
appealing features. First, they have inherit compliance with many applications in 
large-scale road networks, in which aggregate quantities such as traffic counts, flows, 
and space-mean travel speed are more important than characteristics of individual 
vehicles. Second, kinematic wave models can generally be written into a system of 
hyperbolic partial differential equations, or hyperbolic conservation laws. Thus, one 
can better understand the formation and structure of a traffic phenomenon on a road 
network through theoretical analysis of these equations. Finally, there exist many 
sound, efficient numerical methods for solving hyperbolic conservation laws, and one 
can carry out efficient and trustful simulations of large-scale road networks. 



1.2.2 Higher-order models 



Among many efforts to extend the LWR theory to capture instabilities in pr actical 



traffic flow , one direction 



f|l97ll) and 



eads to higher-order, or nonequilibrium, models. 



Payne 



WhithamI (|1974( ) introduced a momentum equation to capture the change 



in travel speed in addition to the traffic conservation equation. Equation 11.1 



vt + vv^ H 

P 



V{p) - V 



(1.3) 



Here the constant Cq is the traffic sound speed, the source term ^^p^'"" is called a 
relaxation term, and r is the relaxation time. With the momentum equation, the 
PW model attempts to model driver behavior by accounting for drivers' anticipation 
and inertia. On e can show that the LWR model is an asymptotic approximation of 



the PW model (|Schochetl . 



19881) . 



The PW model, as a second-order system of hyperbolic c onservation laws with a 



source term, can be numerically solved by Godunov methods (j.Tin and Zhang 



2001b[ l. 



With simulations it was shown that, in addition to modeling stable traffic like the 
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LWR model, th e PW m odel is capable of modeling the formation of vehicle clusters 



( Jin and Zhang . 



2003a|). However, the PW model has drawn some criticism since 



it allows wave solutions with speed highe r than vehicle tr avel speed and may yield 



back-traveling (or negative-speed) results (|Daganzo 



1995c 



Another nonequilibrium traffic flow model is due to (jZhane 



1998 . 



1999 



200G 



2001a|). In this model, a modifled momentum equation is included, 

vt + vv^ + ^^^^^p^ = . 1.4 

P r 

This model bears shortcomings similar to the PW model, yet differs from the latter 
in that the sound speed c = pV'{p) varies with respect to traffi c den sitv p- Moreover, 



it is always stable and therefore acts like the LWR model 



In order to avoid the negative-speed drawbacks of the PW model, 



Aw and Rascle 



( 20001 ) identified a number of principles and proposed a satisfactory momentum equa- 
tion in the following form: 

V{p) - V 



vt + {v - pp'{p))v^ 



T 



;i.5) 



where p{p) is a pressure law and is increasing. Further, the Riemann problems were 
discussed for the model. However, how p{p) is related to driver-behavior was not 
specified in their study. 



In the same spirit, IZhand (|2002t ) derived a model similar to Aw and Rascle's 
model from a car-following model: 

Vip) - V 



Vt + {v + pV'{p))v, 



;i.6) 



which is in the framework of Equation 11.51 with p{p) = —V{p). In this model, the 
definition of p{p) is derived from a car-following model. Both models. Equation 11.51 
and Equa tion 11.61 no longer admit wave solutions faster than traffic and avoids back- 



traveling (jAw and Rascle 



2000; 



Zhang, 



20021 ). However, they are always stable 
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2003^ and, therefore, lose the PW model's abihty to simulate unstable traffic and 



vehicle clusters. 



Also in the framework of Equation ll.5( another model was proposed in (jjiang et al 



2002) as 



Vt + (V - co)v^ 



V{p)-v 



;i.7) 



where the pressure function p{p) = cq In p. Like the PW model, this model is unstable 
under some traf fic conditions, but yie lds non-physical solutions when traffic is in 



unstable region 
More comp. 



Tin and Zhang. 



2003dl ). 



i cated traffic systems, for example, with highway junctions ()Liu et al. 
, l2000[ ). have been studied with non-equilibrium models, primarily 
the PW model. Moreover, several simulation packages are also based on the PW 
model. In this dissertation, however, we focus on the LWR model, i.e., the first-order 



continuum model. 



1.3 Fundamental diagrams 



In the kinematic wave models, fundamental diagrams capture constraints on a traffic 
system such as road characteristics, vehicle type, driver's behavior, weather condi- 
tions, and traffic rules. Therefore, the success of such models rely on the accuracy of 
the fundamental diagram. 

One typical relationship that has been obs erved between fiow rate and occupancy 



is shown in Figure IT!^ from (jHall et al 



198(j ). It's generally assumed that the equi- 



librium travel speed V{p) is decreasing with respect to traffic density; i.e., V'{p) < 0, 
and the fundamental diagram Q{p) = pV{p) is concave; i.e., Q"{p) < 0. 

Due to their importance, many fundamental diagrams have been proposed since 
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Figure 1.2: A typical occupancy-flow rate relationship (jHall et al 



?D ^-O Id » CO TH 



the early era of traffic engineering practice. Some traditional speed-density relations 
are listed in Table ITTTl 

More rece nt examples are the following: the triangular fundamental diagram 



( Newell 



19931) . 



when under-critical, 



Vf (— — 1) otherwise, 



(1.8) 



a non - convex fundamental diagram (jKerner and Konhausei 



19981 ) 



1994; 



Herrmann and Kernel . 



V{p) = 5.0461 [(1 + exp{[p/pj - 0.25]/0.06})"^ - 3.72 x 10"^] //r, (1.9) 
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Functions 




Greenshields fl935) 


Vfd O/ On) 


Greenbere: (1959) 


vq \n(pi/p) 

^ \r Jit/ 


Newell ri961) 


Vf [l-exp (M(i_p^./p))j 


Underwood fl961') 


t;/exp(-p/pj) 


Drake et al. (^9(]7) 


^/exp -|(^) 


Pines (^967) 





Table 1.1: Traditional speed-density relationship functions 



and the exponential form (|Del Castillo and Benitez 
V{p) = Vf ll — exp 




1995 



l-exp^(^-l, 

Vf p 



;i.io) 



The aforementioned fundamental diagr ams are all co ntinuous functions. In addi- 



tion, to capture "two-capacity" phenomena (jBank; 



19911 ) , a disco ntinuous fundamen - 



tal diag ram in "r e verse- lambda" shape has been proposed (e.g., 



Koshi et al. 



19831) . 



Finally, iDaganzol (|l997t ) has proposed a novel, two-regime fundamental diagram for 
differentiated road and vehicle types. 

In this dissertation, we focus on studying network inhomogeneities such as merges, 
diverges, and other junctions. Thus, we will consider continuous fundamental dia- 
gram, mostly of the triangular shape. 



1.4 Motivation for the dissertation research 

While driving in "stop and go" traffic, many people may have wondered how such 
congestion is formed, propagated, and diffused. As we know, the network structure 
is an important factor in determining the characteristics of traffic congestion. For 
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example, congestion generally forms at merging junctions and propagates through 
diverges. Therefore, traffic conditions on a road network have to be considered as a 
whole. 

The kinematic wave theory is a good candidate for studying vehicular traffic 
dynamics due to its theoretical rigor, numerical soundness, and computational ef- 
ficiency. However, its potential has yet to be fully explored. Theoretically, there 
has been little progress i n findi r ig kinematic wa ves in network traffic systems since 



( Lighthill and Whitham 



1955b : 



Richards 



19561 ). and traffic congestion in road net- 



works has not been well understood in the framework of kinematic waves. Numer- 
ically, little attention is paid to the stability and convergence of solution methods. 
Computationally, many simulation models based on kinematic wave theory still keep 
track of individual vehicles and lose a part of the kinematic wave models' efficiency. 
Due to the aforementioned limitations, kinematic wave theories of network traffic flow 
have yet to be improved so that they can be applied in solving practical problems, 
e.g., in Intelligent Transportation Systems. 

In the framework of kinematic wave theories, we study network traffic dynam- 
ics and provide a theoretically rigorous, numerically reliable, and computationally 
efficient simulation model for understanding and mitigating congestion in large-scale 
road networks. In order to fully explore the advantages of the kinematic wave the- 
ory, we consider the following four aspects of traffic modeling in this dissertation. 
First, we study traffic dynamics at crucial components of a road network, such as 
link bottlenecks, merges, and diverges, as well as different types of vehicles. In ad- 
dition, we consider the dynamics of additive mult i- commodity traffic. These studies 
will offer a more rigorous understanding of traffic dynamics in terms of kinematic 
waves. Second, we discuss systematically with numerical methods for computing 
commodity-specified fluxes through junctions in the framework of supply-demand 
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method (jDaganzo 



1995a : 



Lebacque 



19961). Third, we develop a macroscopic simu- 



lation platform of multi-commodity traffic on a road network. Finally, we apply this 
simulation model to study traffic on a road network and probe traffic phenomena 
related to network topology and route choice behavior. 

The structure of this dissertation is as follows. In Chapter 2, we study the 
kinematic waves on inhomogeneous links. In Chapter 3, we discuss kinematic wave 
models of merging traffic, which are based on the supply-demand method. In Chapter 
4, we propose a new kinematic wave theory for diverging traffic. In Chapter 5, we 
discuss theoretical and numerical solutions for mixed-type traffic. In Chapter 6, we 
develop a the multi-commodity kinematic wave simulation platform for network traffic 
flow. In Chapter 7, we apply this simulation model to study some traffic phenomena 
in road networks. Finally, we summarize our research results and present future 
research plans in Chapter 8. 



Chapter 2 



Kinematic wave traffic flow model 
of inhomogeneous links 



2.1 Introduction 



he kin emat ic wave traffic flo w model of LWR was introduced bv lLighthill and Whitham 



(1955a) and 



Richards 



(19561) for modeling dense traffic flow on crowded roads, where 



the evolution of density p{x, t) and flow-rate q{x, t) over time is described by equation, 



pt + Qx 



0. 



{21] 



This equation follows conservation of traffic that vehicles are neither generated nor 
destroyed on a road section with no entries and exits. 

The conservation equation alone is not sufficient to describe traffic evolution, be- 
cause it does not capture the unique character of vehicular flow — drivers slow down 
when their front spacing is reduced to affect safety. The LWR model addresses this 
issue by assuming a functional relationship between local flow-rate and density, i.e., 
q = f{x,p). This flow-density relation, also known as the fundamental diagram of 
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traffic flow, is often assumed to be concave in p and is a function of the local charac- 
teristics of a road location, such as the number of lanes, curvature, grades, and speed 
limit, as well as vehicle and driver composition. When a piece of roadway is homo- 
geneous; i.e., the aforementioned characteristics of the road are uniform throughout 
the road section, the fundamental diagram is invariant to location x and the LWR 
model becomes 

Pt + f{p). = 0. (2.2) 

In contrast, if a section of a roadway is inhomogeneous, the LWR model is 

pt + f{x,p), = 0. (2.3) 

Here we introduce a more explicit notation, an inhomogeneity factor a(x), into the flux 
function /(x, p) and obtain the following equivalent LWR model for an inhomogeneous 
road 

A + /(a,p). = 0. (2.4) 

This equation is particularly suited for our later analysis of the LWR model for 
inhomogeneous roads (We shall hereafter call Equation 12.21 the homogeneous LWR 
model and Equation 12.41 the inhomogeneous LWR model). 

Both the homogeneous and inhomogeneous LWR models have been studied by 
researchers and applied by practitioners in the transportation community. Note that 
the homogeneous version Equation 12.21 is nothing more than a scalar conservation 
law. Therefore, its wave solu tions exist and are unique under the so-called "Lax 
entropy condition" (|Laxl . Il972i ). These solutions are formed by basic solutions to the 
Riemann problem of Equation l2.21 in which the initial conditions jump at a boundary 
and are constant both upstream and downstream of the jump spot. Nevertheless, 



CHAPTER 2. INHOMOGENEOUS LINK MODEL 



19 



because analytical solutions are difficult to obtain for Equation 12.21 with arbitrary 
initial/boundary conditions, numerical solutions have to be computed in most cases. 
The most often used approximation of Equation 12.21 is perhaps that of Godunov. In 
the Godunov method, a roadway is partitioned into a number of cells; and the change 
of the number of vehicles in each cell during a time interval is the net inflow of vehicles 
from its boundaries. The rate of traffic flowing through a boundary is obtained by 
solving a Riemann problem at this boundary. Besides the Godunov method, there 
are other types of approximations of the homoge neous LWR mod el, and some of them 




are shown to be variants of Godunov's method ffLebacqud . [l 

In contrast to the well researched homogeneous LWR model, the inhomogeneous 
model is less studied and less understood. Of the few effor ts to rigorously solve 



Daganzol ()l995ar ) and iLebacqud ()l99fir ) 



the inhomogeneous LWR model, the works of 
should be mentioned. In his cell transmission model, Daganzo started with a discrete 
form of the conservation equation and suggested that the flow through a boundary 
connecting two cells of a homogeneous road is the minimum of the "sending flow" 
from the upstream cell and the "receiving flow" of the downstream cell. The "sending 
flow" is equal to the upstream flow-rate if the upstream traffic is undercritical (UC) or 
the capacity of the upstream section if the upstream traffic is overcritical (OC); on the 
other hand, the "receiving flow" is equal to the capacity of the downstream section if 
the downstream traffic is UC or the downstream flow-rate if the downstream traffic 
is OC. In the homogeneous case, the boundary flux computed from the "sending 
flow" and the "receiving flow" is the same as that computed from solutions of the 
associated Riemann problem. Since the deflnitions of "sending flow" and "receiving 
flow" can be extended to inhomogeneous sections, Daganzo's method can also be 
applied to the inhomogeneous LWR model. Different from Daganzo, Lebacque started 
his method with the solution of the "generalized" Riemann problem for Equation 



CHAPTER 2. INHOMOGENEOUS LINK MODEL 



20 



12.31 In this work, Lebacque came up with some rules for solving the "generalized" 
Riemann problem. These rules play the same role as entropy conditions. Moreover, 
Lebacque found that the boundary flux obtained from solving the Riemann problem 
is consistent with that from Daganzo's method, and he called Daganzo's "sending 
flow" demand and "receiving flow" supply. 

The methods of Daganzo and Lebacque are streamlined versions of Godunov's 
method for the inhomogeneous LWR model. They hinge upon the definitions of the 
demand and supply functions, which can be obtained unambiguously when /(a, p) is 
unimodal. When /(a, p) has multiple local maximum, or when the traffic flow model 
is of higher order, it is yet to be determined if equivalent demand/supply functions 
exist. Thus, these two methods may not be applicable to solve the LWR model that 



has multiple critical points on its funda mental 



traffic flow , such as the 



model by (|Zhang . 



1998, 



Pavne - Whitham 



199S 



2000 



diagram, nor higher-or der models of 



2001a 



mm 



Whitham 



19741 ) model and a 



. Note that, however, these high er-order 



mode 



mm 



s for homogeneous roads can still be solved with Godunov's method (jZhane . 



In this chapter, we present a new method for solving the Riemann problem for 
Equation 12.41 which can be extended to solve higher-order models. By introducing an 
additional conservation law for a{x), we consider the inhomogeneous LWR model as 
a resonant nonlinear system and study its properties (Section 12. 2j) . We also solve the 
Riemann problem for Equation 12.41 and show that the boundary flux at the location 
of the inhomogeneity is consistent with the one given by Lebacque and Daganzo for 
the same initial condition (Section EH)) • Finally, we demonstrate our method through 
solving an initial value problem on a ring road with a bottleneck, and draw some 
conclusions from our analyses. 
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2.2 Properties of the inhomogeneous LWR model 
as a resonant nonlinear system 

Instead of directly study the inhomogeneous LWR model described by Equation 12.41 
we augment Equation 12 .41 into a system of conservation laws through the introduction 
of an additional conservation law aj = for the inhomogeneity factor a{x), which leads 
to 



Ut + F{U), 



0, 



(2.5) 



where U = {a, p), F{U) = (0,/(a,p)),x E R,t > 0. Without loss of generality, 
we assume the inhomogeneity is the drop/increase of lanes at a particular location, 
and write the fundamental diagram as f{a,p) = pV{^), where v = V{^) is the 
speed-density relation. The results obtained hereafter apply to other types of inho- 
mogeneities, such as changes in grades. 

The inhomogeneous LWR model Equation 12.51 can be linearized as 



Ut + dF{U)U, = 0, 
where the differential dF{U) of the flux vector F{U) is 



(2.6) 



dF 



The two eigenvalues of dF are 







The corresponding right eigenvectors are 

^ a' a ^ a' 



Rn 








^ a ' a 


(-) 






a a a 









, Ri = 






1 



(2.7) 



(2i 
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and the left eigenvector of df /dp as li = 1. 

System Eg nation 12 .51 is a non-strictly hyperbolic system, since it can happen that 
Ai = Aq. We consider a traffic state = {a^,p^) in this system as critical if 



Ai(f/*) = 0; 



(2.9) 



i.e., at critical states, the two wave speeds are the same and system Equation 12.51 is 
singular. For a critical traffic state f/* we also have 



fpp<0 



(2.10) 



since flow-rate is concave in traffic density, and 
d 



da 



^)V'(^)|^, 

a a 



a a 



(2.11) 



A consequence of properties Equation I2.1UI and Equation 12.111 is that the lin- 
earized system Equation 12.61 at f/^, has the following normal form 



Sa 









5a 




+ 










t 


1 







0. 



(2.12) 



System Equation 12 . 1 21 has the solution 6p{x,t) = 6a'{x)t + c, and the solution goes to 
infinity as t goes to infinity. Therefore Equation 12. 121 is a linear resonant system, and 
the original inhomogeneous LWR model Equation 12.51 is a nonlinear resonant system. 

For Equation 12. 5[ the smooth curve F in ?7-space formed by all critical states 
are named a transition curve. Therefore F is defined as 

F = {U\Xi{U)=0}. 

Since Ai(f/) = V{^) + f l^'(^), we obtain 

F = |(a,p)|- = a, where a uniquely solves V{a) + aV'{a) = o| ; (2.13) 
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i.e., the transition curve for Equation 12. 51 is a straight hne passing through the origin 
in [/-space. In Equation I2.13t a is unique since f{a,p) is concave in p. 



The entropy solutions to a nonlinear resona nt system are different 



a strictly hyperbolic system of conservation laws. 



Isaacson and Templd (|l992[ ) proved 



rom those to 



that solutions to the Riemann problem for syste m Equ ation 12 .5j exist and are unique 



Lin et al.l (|l995| ) presented solutions 



with the conditions Equation 12. 91 -Equation l^m 
to a scalar nonlinear resonant system, which is similar to our system Equation 12.51 
except that / is convex in their study. In the next section we apply those results to 



solve the Riemann problem for the inhomogeneous LWR model. 



2.3 Solutions to the Riemann problem 

In this section we study the wave solutions to the Riemann problem for Equation 12.51 
with the following jump initial conditions 



U{x,t = 0) = 




if X < 

(2.14) 



where the initial values of Ul, Ur are constant. For computational purposes, we are 
interested in the average flux at the boundary a; = over a time interval At, which 
is denoted by /q and defined as 

/o = fiUix = 0,t))dt. (2.15) 

The augmented inhomogeneous LWR model Equation l2.5l has two families of basic 
wave solutions associated with the two eigenvalues. The solutions whose wave speed is 
Ao are in the 0-family, and the waves are called 0-waves. Similarly the solutions whose 
wave speed is Ai are in the 1-family, and the waves are called 1-waves. The 0-wave is 
also called a standing wave since its wave speed is always 0. The 1-wave solutions are 
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Figure 2.1: Integral curves 

determined by the solutions of the scalar conservation law pt + f{a,p)x = 0, where a 
is a constant. Corresponding to the two types of wave solutions, the integral curves 
of the right eigenvectors Rq and Ri in ?7-space are called 0- and 1-curves respectively. 
Hence the 0-curves are given by /(t/)=const, and the 1-curves are given by a =const. 
A 0-curve, a 1-curve, and the transition curve F passing through a critical state 
are shown in Figure ITTl where a is set as the vertical axis and p the horizontal axis. 

As shown in Figure ITTI the 0-curve is convex, and the 1-curve is tangent to the 
0-curve at the critical state U^. The transition curve F intersects the 0- and 1-curves 
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at ?7*, and there is only one critical state on one 0-curve or 1-curve. For any point 
U, only one 0-curve and one 1-curve pass it. In Figure 12.11 the states left to the 
transition curve are undercritical (UC) since p/a < a; and the states right to the 
transition curve are overcritical (OC) since p/a> a. 

The wave solutions to the Riemann problem for Equation EUl are combinations of 
basic 0-waves and 1-waves. Since Equation 12.51 is a hyperbolic system of conservation 
law, its wave solutions must satisfy Lax's entropy condition that the waves from 
left (upstream) to right (downstream) should increase their wave speeds so that they 
don't cross each other. For Equation l2.5l as a resonant nonlinear system, an additional 
entropy condition is introduced by Isaacson and Temple, 

The standing wave can NOT cross the transition curve F. (2.16) 

This entropy condition is equivalent to saying that, relative to the apexes of the fun- 
damental diagrams, traffic conditions upstream and downstream of inhomogeneities 
are on the same side. That is, they should be either both UC or both OC. 

With the two entropy conditions, the solutions to the inhomogeneous LWR model 
exist and are unique. The wave solutions for UC left state Ul are shown in Figure 
12. 2( and those for OC left state Ul are shown in Figure 

In the remaining part of this section, we discuss wave solutions to the Riemann 
problem for Equation 12.51 present the formula for the boundary flux /q related to 
each type of solution, summarize our results and compare them with those found in 
literature. 

2.3.1 Solutions of the boundary fluxes 

When Ul = {cll,Pl) is UC; i.e., Pl/cll < «, where a is defined in Equation I2.13| 
we denote the special critical point on standing wave passing Ul as f/*. Thus, as 
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500 1000 1500 



P 

Figure 2.2: The Riemann problem for Ul left of F 

shown in Figure I2.2| the ?7-space is partitioned into three regions by DU^, OU-^ 
and f/=,C, where DU^ = {{a, p)\a = a^,p < p*}, Of/* = F fl {0 < p < p*} and 
U^C = {{a, p)\f{a, p) = fiUij^p > p*}. Related to different positions of the right 
state Uji in the f/-space, the Riemann problem for Equation 12 . 51 with initial conditions 
Equation 12^3 has the following four types of wave solutions. For each type of solutions 
we provide formula for calculating the associated boundary flux /q. 

Type 1 When Ur is in region ABUlU^DA shown in Figure 1221 i-e., 

f{UR) < f{U,) = /(Ul), PR/an < a and ur > a„ (2.17) 
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1500 



Figure 2.3: The Riemann problem for Ul right of F 



wave solutions to the Riemann problem are of type 1. These solutions consist 
of two basic waves with an intermediate state Ui = Pl|/(afl,pl)=/(c/,)=/(c/L))• 
Of these two waves, the left one {Ul, Ui) is a standing wave, and the right one 
{Ul, Ur) is a rarefaction wave with characteristic velocity Ai(a, p) > 0. 

From Figure 12.21 we can see that the Riemann problem may admit this type 
of solutions when > aji or a/, < a^, i.e., when the road merges or diverges 
at X = 0. Here we present an example of this type of solutions in Figure 12. 4[ 
where the roadway merges at a; = 0. In the case when the roadway diverges at 
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Figure 2.4: An example for wave solutions of type 1 for Equation 12.51 with initial 
conditions Equation 12.141 

a; = 0, we can find similar solutions. 

From Figure ITU we obtain the boundary flux /q = /{Ul) = fiU^) for wave 
solutions of type 1. 

Type 2 When Ur is in region BUlU^^CB shown in Figure YI7]\ i.e., 

f{Ua) > = /(f/L), (2.18) 

wave solutions to the Riemann problem are of type 2. These solutions consist 
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Figure 2.5: An example for wave solutions of type 2 for Equation 12.51 with initial 
conditions Equation 12.141 

of two basic waves with an intermediate state Ui = Pl|/(afl,pl)=/(^7,)=/(t/L))• 
0f these two waves, the left {Ul, Ui) is a standing wave, and the right (f/i, Ur) 
is a shock wave with positive speed a = f^^^'>~f^^*'> > Q. 

From Figure 12. 2[ we can see that the Riemann problem may admit this type 
of solutions when the downstream traffic condition Ur is UC or OC, or the 
roadway merges or diverges at a; = 0. Here we present an example of this type 
of solutions in Figure where the downstream traffic condition is OC and the 
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roadway merges at x = 0. Similar solutions can be found for other situations 
that satisfy Equation 12.181 

From Figure 1231 we obtain the boundary flux /q = /(f/i) = f{U*) for wave 
solutions of type 2. Here we have the same formula as that for wave solutions 
of type 1. 

Type 3 When Urisvo. region OU^CO shown in Figure i.e., 

f{UR)<f{U,)=f{UL), PR/aR>a, (2.19) 

wave solutions to the Riemann problem are of type 3. These solutions consist 
of two basic waves with an intermediate state Ui = (^l, Pi|/(a£,pi)=/{c/H))- Of 
these two waves, the left one {Ul,Ui) is a shock wave with negative speed 
a = '^^'^plZp^^'^^ < O5 ^iid the right one (Ui, Ur) is a standing wave. 

From Figure we can see that the Riemann problem may admit this type of 
solutions when the roadway merges or diverges at x = 0. Here we present an 
example of this type of solutions in Figure 12.61 where the roadway merges at 
a; = 0. In the case when the roadway diverges at a; = 0, similar solutions can 
be found. 

From Figure EEl we obtain the boundary flux /q = /{Ur) for wave solutions of 
type 3. 

Type 4 When t/ij is in region OU^,DO shown in Figure I2I21 i-e., 

/{Ur) < f{U^) = JXUl), pR/aR < pja^ and aR < a„ (2.20) 

wave solutions to the Riemann problem are of type 4. These solutions consist 
of three basic waves with two intermediate states: Ui = {^l, Pi\f{aL,pi)=f{U2)) 
and U2 = {ttR, P2\p2/aji=a)- Of these three waves, the left one {Ul, Ui) is a shock 
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Figure 2.6: An example for wave solutions of type 3 for Equation 12.51 with initial 
conditions Equation 12.141 

wave with negative speed a = '^^^^J^Z'^J'^^^ < 0; the middle one {Ui,U2) is a 
standing wave with zero speed, and the right one {U2, Ur) is a rarefaction wave 
with characteristic velocity Ai(a, p) > 0. 

From Figure I2I21 "we can see that this type of solutions are admitted only when 
the roadway merges at x = 0. Here we present an example of this type of 
solutions in Figure EH 

From Figure EZl we obtain the boundary flux /q = /(t/2) for wave solutions of 



CHAPTER 2. INHOMOGENEOUS LINK MODEL 



32 




Figure 2.7: An example for wave solutions of type 4 for Equation 12.51 with initial 
conditions Equation 12.141 

type 4. 

When Ul = (a^, Pl) is OC; i.e., Pl/il > we denote the special critical point on 
1-wave curve passing Ul as U^; i.e., = {cil, P*\p,/aL=a)- Thus, as shown in Figure 
12.31 the f/-space is partitioned into three regions by three curves DU^, = {a = a* = 
aL,0 < p < p*}, OU^, = {0 < a < a^, p = aa} and U^C = {a > a^,f{a,p) = f{U^). 
Related to different positions of the right state Un in the t/-space, the Riemann 
problem for Equation 12.51 with initial conditions Equation 12.141 has the following six 
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types of wave solutions. For each type of solutions we provide formula for calculating 
the associated boundary flux /q . 

Type 5 When Ur resides in region ABU^^DA shown in Figure 1231 i-e., 



wave solutions to the Riemann problem are of type 5. These solutions con- 
sist of three basic waves with two intermediate states: Ui = f/* and U2 = 
{o-Ri p2\f{U2)=f(u,))- Of these three waves, the left one {Ul,Ui) is a rarefaction 
wave with negative characteristic wave velocity Ai(a, p), the middle one (f/i, U2) 
is a standing wave and the right one {U2, Ur) is a rarefaction wave with positive 
characteristic velocity \i{a,p). 

From Figure EISl we can see that this type of solutions are admitted only when 
the roadway diverges at x = 0. Here we present an example of this type of 
solutions in Figure I2IH1 

From Figure ITHl we obtain the boundary flux /q = /(t/2) for wave solutions of 
type 5. 

Type 6 When Ur resides in region BU^CB shown in Figure IT^ i.e., 



solutions to the Riemann problem are of type 6. These solutions consist of three 
basic waves with two intermediate states: Ui = and U2 = {dR, P2\f{U2)=f{Ut))■ 
0^ these three waves, the left one {Ul, Ui) is a rarefaction wave with negative 
characteristic velocity Xi{a,p), the middle one {1/1,1/2) is a standing wave and 
the right one {U2, Ur) is a shock wave with positive speed a = . 



fiUR)<f{U^), pR/aR< a and aR>a^ = ql, 



(2.21) 



f{UR) > f{U,), 



(2.22) 
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Figure 2.8: An example for wave solutions of type 5 for Equation 12.51 with initial 
conditions Equation 12.141 

From Figure 12. 3[ we can see that this type of solutions may be admitted when 
the downstream traffic condition is UC or OC; However, they are admitted only 
when the roadway diverges at a; = 0. Here we present an example of this type 
of solutions in Figure IT^ where the downstream traffic condition is OC. In the 
case when the downstream traffic condition is UC, we can find similar solutions. 

From Figure we obtain the boundary flux /q = /(t/2) for this type of wave 
solutions. Here we have the same formula as that for wave solutions of type 5. 
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Figure 2.9: An example for wave solutions of type 6 for Equation 12.51 with initial 
conditions Equation 12.141 

Type 7 When Uji resides in region CU^FUlEC shown in Figure ESI i-e., 

f{UL) < JXUr) < f{U,) and PR/an > a, (2.23) 

wave solutions to the Riemann problem are of type 7. These solutions consist of 
two basic waves with an intermediate state Ui = {cll, Pi\f(Ui)=f{Uii))- Of these 
two waves, the left one {Ul, Ui) is a rarefaction with negative characteristic 
velocity Ai(a,p), and the right one (f/i, Ur) is a standing wave. 

From Figure 1231 "we can see that the Riemann problem may admit this type of 
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Figure 2.10: An example for wave solutions of type 7 for Equation 12.51 with initial 
conditions Equation 12.141 

solutions when the roadway merges or diverges at x = 0. Here we present an 
example of this type of solutions in Figure 12.101 where the roadway diverges at 
a; = 0. In the case when the roadway merges, we can find similar solutions. 

From Figure EUni we obtain the boundary fiux /q = /(t/fi) for wave solutions 
of type 7. 




Figure 2.11: An example for wave solutions of type 8 for Equation 12.51 with initial 
conditions Equation 12.141 

Type 8 When Un locates in region FUlEOF shown in Figure 1231 i-e., 

fiUa) < JXUl) < and pa/an > a, (2.24) 

wave solutions to the Riemann problem are of type 8. These solutions consist of 
two basic waves with an intermediate state Ui = {cll, Pi\f{Ui)=f{UR))- Of these 
two waves, the left one {Ul, Ui) is a shock with negative speed a = '^^'^^jZ.p^^^K 
and the right one {Ui, Ur) is a standing wave. 

Like in the previous case, the Riemann problem may admit this type of solutions 
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when the roadway merges or diverges at x = 0. Here we present an example of 
this type of solutions in Figure E.lll where the roadway diverges at a; = 0. In 
the case when the roadway merges, we can find similar solutions. 

From Figure 1^.111 we obtain the boundary flux /q = /(fZ/j) for wave solutions 
of type 8. The formula is the same as that for wave solutions of type 7. 

Type 9 When Ur resides in region DU^FGD shown in Figure 1231 i-e., 

fiUi) < fiUn) < f{U^), pa/aR < a and a^j < a, = a^, (2.25) 

wave solutions to the Riemann problem are of type 9. These solutions consist 
of three basic waves with two intermediate states: Ui = {(^l, Pi\f(Ui)=f{U2)) 
and U2 = {ciR, P2\p2/aii=a)- Of these three waves, the left one {Ul,Ui) is a 
rarefaction with negative characteristic velocity Ai(a, p), the middle one {Ui, U2) 
is a standing wave, and the right one {U2, Ur) is a rarefaction with positive speed 
Ai(a,p). 

From Figure ESI we can see that this type of solutions are admitted only when 
the roadway merges at x = 0. Here we present an example of this type of 
solutions in Figure 12.121 

From Figure 0.12| we obtain the boundary flux /q = /(t/2) for wave solutions 
of type 9. 

Type 10 When Ur resides in region GFOG shown in Figure EHl i-e., 

/(Ur) < /(Ul) < fiU^), pR/aR<a and aR<a^=aL, (2.26) 

wave solutions to the Riemann problem are of type 10. These solutions consist 

of three basic waves with two intermediate states: Ui = {cll, Pi\f{Ui)=f{U2)) 

U2 = {dR., P2\p2/ait=a)- Of thcsc three waves, the left one {Ul,Ui) is a shock 
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Figure 2.12: An example for wave solutions of type 9 for Equation 12.51 with initial 
conditions Equation 12.141 

with negative speed, the middle one (f/i, U2) is a standing wave, and the right 
one {U2, Ur) is a rarefaction wave with positive characteristic velocity Xi{a,p). 

Like in the previous case, this type of solutions are admitted only when the 
roadway merges at a; = 0. Here we present an example of this type of solutions 
in Figure ITTSl 

From Figure 1^.131 we obtain the boundary flux /q = f{U2) for wave solutions 
of type 10. Here we have the same formula as that for wave solutions of type 9. 
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Figure 2.13: An example for wave solutions of type 10 for Equation 12.51 with initial 
conditions Equation 12.141 

2.3.2 Summary 

In each of the 10 cases discussed above, the boundary flux /q is equal to one of 
the following four quantities: the upstream flow-rate /{Ul), the downstream flow- 
rate fiUji), the capacity of the upstream roadway and the capacity of the 
downstream roadway f^°-^. For wave solutions of type 1 and 2, the boundary flux 
is equal to the upstream traffic flow-rate; i.e., /q = fiUij- For wave solutions of 
type 3, 7 and 8, the boundary flux is equal to the downstream traffic flow-rate; i.e.. 
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No. 


left state Ul 


right state Ur 


fo 


1 


UC 


/{Ur) < fiUi), aR > a*, Pr/ur < a 


fiUL) 


2 


UC 


/{Ur) > /{Ul) 


fiUL) 


3 


UC 


IXUr) < /{Ul), pR/aR > a 


fiUR) 


4 


UC 


/{Ur) < /{Ul), PR/aR < a, gr < a* 


£max 
JR 


5 


oc 


f{UR) < Jl"'', aR> ai, pR/aR < a 


X'max 

Jl 


6 


oc 


/(Ur) > fr" 


Jl 


7 


oc 


/{Ul) < /{Ur) < /r^ PR/aR > a 


JXUr) 


8 


oc 


/(Ur) < fiUL), Pr/ur > a 


IXUr) 


9 


oc 


/{Ul) < f{UR) < f'^'''', PR/aR <a,aR< 


fmax 
JR 


10 


oc 


/{Ur) < /(Ul), pR/aR < a, aR < aL 


fmax 
JR 



Table 2.1: Solutions of the boundary fluxes /g 



/o — fi^R)- For wave solutions of type 4, 9 and 10, the boundary flux is equal to the 
capacity of the downstream roadway; i.e., /q = Jr"'^. For wave solutions of type 5 and 
6, the boundary flux is equal to the capacity of the upstream roadway; i.e., /q = /™"^. 
In Table 123 the boundary fluxes are listed for the 10 types of wave solutions to the 
Riemann problem, as well as the conditions when the Riemann problem admit those 
solutions. 

Note that when = aR, i.e., when Equation 12.41 becomes a homogeneous LWR 
model, wave solutions and the solutions of the boundary fluxes provided here are the 
sam e as those for th e homogeneous LWR model. 



Lebacqud ()l99fir ) studied the Riemann problem of the inhomogeneous LWR for 
Equation 12.31 He classified the problem according to two criteria. The first criterion 
is to compare capacity of the upstream cell and that of the downstream cell. For the 



roadway with variable number of lanes, it is equivalently to compare the number of 
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Conditions 


Solutions by Lebacque 


Types 


Our solutions 


ai < an, Ul UC, Ur UC 


fiUL) 


1 


fiUL) 


aL < aR, Ul UC, Ur OC 


mm{ f{UL),f{UR)} 


2 or 3 


/{Ul) or f{UR) 


flL < aR, Ul OC, Ur UC 


£max 

Jl 


5 or 6 


£max 

Jl 


ciL < aR, Ul OC, ?7/? OC 




6, 7 or 8 




aL > aR, Ul UC, Ur UC 


min{/--,/(f/i) 


1 or 4 


/(Ul), rr 


flL > a^, Ul UC, f/,? OC 


min{/(?7i),/(t/^,)} 


2 or 3 


/(Ul) or /(Ur) 


flL > aR, Ul OC, ^7^ UC 


/max 


9 or 10 


fmax 
JR 


aL > aR, Ul OC, Ur OC 




7 or 8 


fiUR) 



Table 2.2: Comparison with Lebacque 's results 

lanes of the upstream cell and that of the downstream cell. The second criterion is 
to consider whether the upstream and downstream traffic conditions are UC or OC. 
With these criteria, he discussed 8 types of waves solutions to the Riemann problem 
and obtained the formula for the boundary flux related to each type of solutions. The 
conditions for those types of wave solutions as well as the formulas related to those 
types of solutions are listed in Table IT^ Under each of those conditions, the Riemann 
problem may admit different types of solutions discussed in Section [2.3. II The types 
of solutions and our related formulas for the boundary flux are also presented in Table 
12.21 From this table, we can see that our results are consistent with those provided 
by Lebacque, although the Riemann problem is solved through different approaches. 

The consistency of our results with existing results can also be shown by intro- 
ducing a simple formula for the boundary flux. If we define the upstream demand 
as 



/(f^L), PL/aL<a 
fr". Pi/aL > a 



(2.27) 
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and define the downstream supply as 

/« - I " (2.28) 

then the boundary flux can be simply computed as 

/* = min{/2,/*}. (2.29) 

I 1 I 

Not e that f* =fiU^ )- Formula Equation 1221 was also provided bv iDaganzd ()l995ar ) 



and 



Lebacqud (|l996| ). 



2.4 Simulation of traffic flow on a ring road with 
a bottleneck 

2.4.1 Solution method 

The augmented inhomogeneous LWR model, expressed in conservation form Equa- 
tion EIH can be solved efficiently with Godunov's method under general initial and 
boundary conditions. In Godunov's method, the roadway is partitioned into N cells 
and a duration of time is discretized into M time steps. In a cell i, we approximate 
the continuous equation Equation 12.51 with a finite difference equation 

At Ax 

whose component for p is 

m+l m f* , _ f* 

Pi - Pi Ji+i/2 ^ Q ^2 

At Ax ' 

where denotes the average of p in cell i at time step m, similarly p^"*"^ is the 
average at time step m + l; /*_i/2 denotes the flux through the upstream boundary 
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of cell i, and similarly f*_^i/2 denotes the downstream boundary flux of cell i. In 
Equation 12.311 the boundary flux f*__i/2 related to solutions to a Riemann problem 
for Equation 12.51 with the following initial conditions: 

(2.32) 

X> Xi-i/2 

which have been discussed in Section 12.31 
2.4.2 Numerical results 

We use the approximation developed earlier to simulate traffic on a ring road. The 
length of the ring road is L = 800/ = 22.4 km. The simulation time is T = 500r 
= 2500 s = 41.7 min. We partition the road [0,L] into = 100 cells and the time 
interval [0, T] into K = 500 steps. Hence, the length of each cell is Ax = 0.224 km 
and the length of each time step is At = 5 s. Since |A*| < = 5//r, we find the CFL 



(Courant et al 



192a ) condition number 



maxlAJ^ < 0.625 < 1. 
' ' Ax ~ 

Moreover, we adopt in this simulation the fundamental diagram used in (Kerner and 
Konhauser, 1994; Herrmann and Kerner, 1998) with the following parameters: the 
relaxation time r = 5 s; the unit length / = 0.028 km; the free fiow speed Vf = 5.01/t 
= 0.028 km/s = 100.8 km/h; the jam density of a single lane pj = 180 veh/km/lane. 
The equilibrium speed-density relationship is therefore 

-1 



v^{p,a{x)) = 5.0461 



fl+exp|[— 4 0.251/0.06}^ -3.72x10" 



where a{x) is the number of lanes at location x. The equilibrium functions V{p, a{x)) 
and f{p,a{x)) are given in Figure imi 
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Figure 2.14: The Kerner-Konhauser model of speed-density and flow-density relations 

The first simulation is about the homogeneous LWR model. Here we assume that 
the ring road has single lane everywhere; i.e., a{x) = 1 for any x G [0, L] , and use a 
global perturbation as the initial condition 

p{x,0) = ph + Aposin^^, sG[0,L], 
v{x,0) = t;,(p(a;,0), 1), xe[0,L], 

with ph = 28 veh/km and Apo = 3 veh/km (the corresponding initial condition 
Equation 12.331 is depicted in Figure ITT^j) . 

The results are shown in Figure 12. 16^ from which we observe that initially wave 




Figure 2.15: Initial condition Equation 12.331 with ph = 28 veh/km and Apo = 3 
veh/km 

interactions are strong but gradually the bulge sharpens from behind and expands 
from front to form a so-called A^-wave that travels around the ring with a nearly fixed 
profile. 

In the second simulation we created a bottleneck on the ring road with the fol- 
lowing lane configuration: 

' 1, xe [320/,400n, 
a{x) = { L , ^2.34) 

2, elsewhere . 
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p/Pj v/(i/T) q/lPjiz-c) 




0.14 0.15 0.16 0.17 4 4.1 4.2 4.3 0.62 0.64 0.66 0.68 



Figure 2.16: Solutions of the homogeneous LWR model with initial condition in Figure 

Em 

As before, we also use a global perturbation as the initial condition 

p(x,0) = a{x){ph + Apo sin ^^), xG[0,L], 
f(x,0) = f*(p(x, 0), a(a;)), xG[0,L], 

with ph = 28 veh/km/lane and Apo = 3 veh/km/lane (the corresponding initial 
condition Equation 12.351 is depicted in Figure ITTTj) . 

The results for this simulation are shown in Figure I?. 18^ and are more interesting. 
We observe from this figure that at first flow increases in the bottleneck to make 




Figure 2.17: Initial condition Equation 12.331 with ph = 28 veh/km/lane and Apo = 3 
veh/km/lane 

the bottleneck saturated, then a queue forms upstream of the bottleneck, whose 
tail propagates upstream as a shock. At the same time, traffic emerging from the 
bottleneck accelerates in an expansion wave. After a while, all the commotion settles 
and an equihbrium state is reached, where a stationary queue forms upstream of 
the bottleneck, whose in/out flux equals the capacity of the bottleneck. Similar 
situations can be observed in real world bottlenecks, although queues formed at such 
bottlenecks rarely reach equilibrium because, unlike in the ring road example, their 
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traffic demands change over time. Therefore, we observe queues forming, growing, 
and dissipating at locations with lane drops, upward slopes, or tight turns. Sometimes 
queues formed at a bottleneck can grow fairly long, to the extent that they entrap 
vehicles that do not use the bottleneck. Under such situations, we can implement 
various types of control strategies, such as ramp metering, to control the extent of 
the bottleneck queues so that they do not block vehicles that wish to exit upstream 
of the bottleneck. For this purpose the numerical method presented here can be used 
to help model and design effective control. 

2.5 Concluding remarks 

We studied the inhomogeneous LWR model as a nonhnear resonant system. The non- 
linear resonance arises when the two characteristics of the augmented LWR model 
coalesce. Critical states and a transitional curve F can be defined in the U space 
based on the behavior of these characteristics, which are in turn used to solve the 
Riemann problem for the inhomogeneous LWR model. It is found that, under the 
entropy conditions of Lax and of Isaacson and Temple, there exist ten types of wave 
solutions. Formulas for computing the boundary fiuxes related to different types of 
wave solutions were also obtained. These formulas, after translated into the sup- 
ply/demand framework, are found to be consistent with those found in literature. 
For problems with general initial/boundary conditions, the method of Godunov was 
applied to solve the inhomogeneous model numerically. 

The method presented here can be extended easily to model more complicated 
situations, such as multiple inhomogeneities. Suppose at location x, there are i = 
1, • • • , n types of inhomogeneities, such as changes in number of lanes, grade, and cur- 
vature. We introduce an inhomogeneity vector a{x) — {ai{x), a2{x), • • • , a„(x))'^, and 
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200 400 600 200 400 600 200 400 600 

x/l x/l x/l 



Figure 2.18: Solutions of the inhomogeneous LWR model with initial condition in 
Figure ITT7I 

express the flow-density function as f{a{x),p). Then the conservation law becomes 

pt + f{aix),p)^ = 0, 
a{x)t = 0, 

and this higher-dimensional nonlinear resonant system can be solved in a similar way. 
It is worth mentioning that the augmentation approach taken in this chapter also 
applies to higher-order traffic flow models for inhomogeneous roads. 



Chapter 3 

Kinematic wave traffic flow model 
of merging traffic 



3.1 Introduction 

For developing advanced traffic control strategies, dynamic traffic assignment (DTA) 
algorithms, and other technologies in Advanced Traffic Management Systems (ATMS) 
and Advanced Traveler Information Systems (ATIS), traffic engineers need the assis- 
tance of network traffic flow models that can capture system-wide features of traffic 
dynamics and are computationally efficient for a network of realistic size. The kine- 
matic wave model is a promising candidate for these tasks since it provides a realistic 
description of dynamic traffic phenomena in the aggregate level in terms of expansion 
and shock waves and as such is highly efficient for simulating traffic dynamics in a 
large network. 



In the sem inal kinematic wave model by 



Lighthill and WhithamI ()l955b| ) and 



Richardsl (jl956l ). a.k.a. the LWR model, how a disturbance in traffic propagates 



through a link was thoroughly studied. To model traffic dynamics on a network with 
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the kinematic wave model, however, one needs to carefully study traffic dynamics at 
a merge, a diverge, or other compone nts of a ne t work. T he kinematic wave m o dels o f 



mer ging traffic have been studied by 



and 



Daganzd (|l995al) 



Holden and Risebrd (|l995r ). 



Lebacqud (|l996| ). In the model by Holden and Risebro, traffic flows through a 



merge are determined by an optimization problem. However, the physical meaning 
and the objective function of the optimization problem are not known or supported 
by observations. On the other hand, the models by Daganzo and Lebacque are based 
on the definitions of the local traffic supply and demand and can be considered as 
reasonable extensions of the kinematic wave model of link traffic flow. In this chapter, 
we will examine the latter models so that they can be better understood, more easily 
calibrated, and more efficiently applied in simulation. 

As we know, the LWR model, in which the evolution of traffic density t), flow- 
rate q{x, t), and travel speed v{x, t) is studied in space x and time t, can be written as a 
partial differential equation based on the fact of traffic conservation and the adoption 
of a fundamental diagram. For the purpose of simulation, the LWR model is generally 
written in a discrete form: a link and a duration of time are partitioned into a number 
of cells and time steps respectively, and the increment of the number of vehicles in 
a cell at each time step equals to the difference between the inflow into and outflow 
from that cell during the time step. In the discrete LWR model, to solve the flow 
through a boundary (i.e., the inflow into the downstream cell and the outflow from the 
upstream cell), two equivalent approaches can be use d: in the mathernatical approach, 



one so lves the Riemann problem at that boundary ([Lebacque . 



■Tin and Zhang. 



20n8br ) : in the engineering approach, a.k.a. the supply-demand method, the supply of 
the upstream cell and the demand of the downstream cell are computed ffist and the 
boundary flow is taken as their minimum . Here the concepts of supply and demand 



were first introduced by 



Daganzd ()l995a[ ). but using the terms of "sending flow" and 
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'receiving flow" instead; the terms of "supply" and "demand" were first used by 



Lebacqud (|l996r i. The definitions of demand and supply are as follows: the demand 



of a cell is equal to its fiow-rate when the traffic condition is under-critical (i.e., free 
flow) and its flow capacity when overcritical (i.e., congested); the supply is equal 
to the flow capacity of the cell when the traffic condition is under-critical and the 
fiow-rate when overcritical.^ 

For computing fiows through a merge, including the outfiows from the upstream 
cells and the infiow into the downstream cell, Daganzo (1995) extended the supply- 
demand method as follows: the outfiow from an upstream cell is smaller than or 
equal to its demand, the infiow to the downstream cell is smaller than or equal to its 
supply, and the infiow is equal to the sum of the outfiows in order to preserve traffic 
conservation. In this supply-demand method, the infiow is unique since it is equal to 
the minimum of the supply and total demand. But the outfiow from each upstream 
cell may not be unique. Thus one has to find a way to distribute to each upstream 
cell a fraction of the total outfiow, which is equal to the infiow. Here we call such a 
way of determining the distribution fractions the distribution scheme. 

Lebacque proposed another extension of the supply-demand method: the supply 
of the downstream cell is first distributed as a virtual supply to each upstream cell, 
the outfiow from each upstream cell is the minimum of its demand and virtual supply, 
and the infiow into the downstream cell is equal to the sum of the outfiows. Thus the 
distribution scheme in Lebacque's method is used to determine the fractions of virtual 
supplies, and is more general since more feasible solutions of fiows can be found in 
this method. 

Both Daganzo and Lebacque provided general formulations of the kinematic wave 



^From the definitions of demand and supply, we can see that the flow through a boundary is 
bounded by the capacity. 
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model of merges. Here, we do not intend to extend these formulations. Rather, we are 
interested in the distribution schemes used in these models since a distribution scheme 
is the key to uniquely determine flows through a merge. Since in possible applications 
of a merge model one wants to obtain unique flows under a given situation, the 
distribution schemes are worth a thorough examination. 

At a first glance, the determination of distribution fractions seems to be com- 
plicated since they may be affected by travelers' merging behavior, the geometry 
of the studied merge, traffic capacities, differences between the upstream cells, traf- 
fic conditions, and traffic control. Considering part of these factors, both Daganzo 
(1995) and Lebacque (1996) provided some suggestions on the distribution fractions: 
Lebacque suggested that the distribution fraction of an upstream cell is proportional 
to its number of lanes; Daganzo considered that upstream cells bear different priori- 
ties and hence introduced parameters for priorities in his distribution fractions. Both 
suggestions have their limitations: Lebacque's distribution scheme is very coarse and 
fails under certain situations, while Daganzo's scheme becomes very complicated for 
a merge with three or more upstream links. Moreover, priorities in Daganzo's distri- 
bution scheme vary with fiow levels, which seems to be counter intuitive. Therefore, 
we devote this study to the better understanding of various distribution schemes in 
a merge model, and propose a new distribution scheme which is well-defined, com- 
putationally efficient, and capable of capturing the characteristic differences between 
different branches of the merge. 

In this chapter, we first review the discrete kinematic wave model of merges 
and discuss different formulations of the supply-demand method for computing flows 
through a merge f Section 13. 2|) . In Section [3.3.21 after discussing existing distribution 
schemes, we propose a simple distribution scheme, which incorporates the "fairness" 
condition. In this scheme, the distribution fractions are proportional to traffic de- 
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mands of upstream cells. This scheme is shown to work well in simulations due to its 
many merits: 1) it is capable of capturing the characteristic differences between up- 
stream cells (e.g. the speed difference between the upstream freeway and on-ramps); 
2) it is easy to calibrate because additional parameters such as priorities do not need 
to be explicitly introduced; and 3) it is computationally efficient. In Section we 
present an example of two merging flows and demonstrate in numerical simulations 
that the discrete kinematic wave merge model incorporating the "fairness" condition 
is well-deflned and converges in flrst order. In the conclusion part, we present the 
supply-demand method for computing flows through a diverge and a general junction 
for single- commodity traffic flow, and discuss related future research. 



3.2 The discrete kinematic wave model of merges 
with the supply-demand method 

In the kinematic wave traffic flow model of a road network with a merge, the LWR 
model can be used to describe traffic dynamics of each branch, for which flows through 
the merge can be considered as boundary conditions. Thus, in this section, we 
flrst review the discrete LWR model, the deflnitions of supply and demand, and the 
supply-demand method for computing flows through link boundaries. After reviewing 
the models of merges under the supply-demand framework by Daganzo (1995) and 
Lebacque (1996), we then demonstrate the importance of distribution schemes. At 
the end of this section, we will discuss the properties of existing distribution schemes. 
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3.2.1 The discrete LWR model in the supply-demand frame- 
work 

In the LWR model for each branch of a merge, traffic dynamics are governed by a 
traffic conservation equation, 



Pt + Qx 



0, 



(3.1) 



and an equilibrium relationship between p and q, also known as the fundamental 
diagram. 



q = Qia,p), 



(3.2) 



where a{x) is an inhomogeneity factor, depending on road characteristics, e.g., the 
number of lanes at x. Since q = pv, we also have a speed-density relation: v = 
V{a,p) = Q{a,p)/p. For vehicular traffic, generally, v is non-increasing and q is 



concave in p. Examples o 



be found in (^Newell 



empirical models of speed- a nd flow-density relations can 



Kerner and Konhauser 



19941 ). Related to the fundamental 



diagram, the following deflnitions are used in this chapter: the maximum flow-rate 
at X is called the traffic capacity, and the corresponding density is called the critical 
density; traffic flow is overcritical when its density is higher than the critical density, 
and under-critical conversely. 

From Equation 13.11 and Equation 13. 2t the LWR model can be written as 



pt + Q{a,p), 



0, 



(3.3) 



where < p < pj {pj is the jam density). When a{x) is uniform with respect to loca- 
tion X, the LWR model is called homogeneous. Otherwise it is called inhomogeneous 
Both the homogeneous and inhomogeneous models are hyperbolic systems of conser- 



vation laws. Actually the former, which is a strict hyperbolic conservation law (jLax, 
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1972), is a special case of the latter, a n on-strictly hyperbolic system of conseryation 



1992). Therefore, the 



laws and a resonant nonlinear system (jlsaacson and Temple , 
following discussions for the inhomogeneous LWR model are yalid for any kind of 
links. 

With jump initial conditions, the LWR model Equation 13.31 is solyed by shock 
wayes, expansion wayes, and standing wayes. These waye solutions are unique under 
the so-called "entropy" conditions. Howeyer, solutions of the LWR model with general 
initial and boundary conditions can not be expressed in analytical form, which calls 
for approximate solutions with n umerical method s. One efficient numerical method 



for solying Equation 13.31 is due to 



GodunoyI ()l959( ). In the Godunoy method, the link 



is partitioned into N cells, a duration of time is discretized into M tim e steps, and the 



discre tization of space and time satisfies the Courant-Friedrichs-Lewy (ICourant et al 



19281 ) (CFL) condition so that a yehicle is not allowed to cross a cell during a time 



interyal. Assuming that the spacing Ax and the time step At are constant, is the 
ayerage of p in the cell i at time step m, ({^^^12 ^^"-^ ^+1/2 infiow into and the 

outfiow from cell i from time step m to m + 1 respectiyely, the LWR model Equation 
13.31 for cell i can be approximated with a finite difference equation: 

m+l 



Pi 



Pi 



+ 



m+1/2 _ m+1/2 
1i-l/2 ~ 1i+l/2 



0. 



(3.4) 



At Ax 

In Equation 13. 4| the flow through the hnk boundary Xj_i/2, i.e., q'j'^^y2^, can be 
computed in two approaches. One is from the waye sol utions of the Riemann p roblem 



for Equation 13.31 with the following initial conditions (jjin and Zhang . 



2003b(): 



U{X = Xi_i/2,t = tm) 



Ull^ X < Xi_i/2 



(3.5) 



X > X- 



wher e U = (a, p). Another is the supply-demand method (jDaganzo 



i-1/2 



1995a; 



Lebacque 



19961), in which the flow through a link boundary is the minimum of the traffic demand 
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of its upstream cell and the t raffic supply of i ts dow nstream cell. The two approaches 



were shown to be equivalent (|Jin and Zhang , 



2003br i. However, the method of solving 



the Riemann problem and the supply-demand method have different fates for studying 
traffic dynamics through a merge: there has been no formulation of the Riemann 
problem for merging dynamics in literature, but the supply-demand method has been 
extended and applied in the discrete kinematic wave models of merges. 

In the following, we describe in detail the supply-demand method. Considering 
the link boundary at Xj_i/2, whose upstream and downstream cells are respectively 
denoted by cell i — 1 and cell i, supposing that the traffic densities of the two cells 
are pY^i and pf^ at time step m, Daganzo (1995) and Lebacque (1996) suggested the 
following supply-demand method for computing q^^^l^ ■ First, traffic demand of the 
cell i — \ (called "sending flow" by Daganzo), D™"^^''^ , and traffic supply of the cell 
i (called "receiving flow" by Daganzo), S'™^^''^, are defined by 

^_Li /o QT-^ , cell i — 1 is under-critical, 

I QT^i, otherwise; 

/o QT""^, celH is under-critical, 

= <^ (3.7) 

I Q^, otherwise; 

where Q™"^ is the capacity of cell i, and Qf^ the flow-rate of cell i at time step m. 
The demand can be considered as the maximum flow that can be discharged by the 
cell i — 1 from time step m to m + 1; the supply S^~^^^^ is the maximum flow that 
can be received by the cell i. Thus, the boundary flow satisfles (all superscripts will 
be suppressed hereafter) 

- ""^-^ (3.8) 
qi-i/2 < Si. 

Note that Equation 13.81 admits multiple solutions. To identify the unique boundary 
flow, an additional "optimality" condition, that the actual boundary flow always 
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reaches its maximum, is assumed. Hence, the boundary flow can be simply computed 

by 

gj_i/2 = mm{Di_i,Si}. (3.9) 

Here the "optimahty" condition can be considered as an entropy condition, which 
helps to choose a physical solution out of all feasible solutions. 

3.2.2 The kinematic wave model of merging traffic in the 
supply-demand framework 

In this subsection, we review the kinematic wave model of merging traffic in the 
supply-demand framework. In this type of models, the supply-demand method is 
used to compute flows through a merge. Without loss of generality, we consider 
a merge that connects two upstream cells to one down stream cell. Furthermore, 
we assume that, at time step m, traffic demands of the two upstream cells and the 
traffic supply of the downstream cell are Di, D2, and Sd respectively. We denote the 
outflows from the upstream cells by gi and 52 and the inflow into the downstream 
cell by q from time step m to m -|- 1. Then, according to traffic conservation, we have 

q = qi + q2- 

The basic assumption in the supply-demand method for computing the flows 
through a merge is that the flows, qi, q2, and q, are determined by traffic conditions 
Di, D2, Sd, and/or other characteristics of the merge. Another assumption, as in 
the supply-demand method for computing the flow through a link boundary, is the 
optimality condition. Two types of optimality conditions have been proposed: one is 
that the total flow q reaches its maximum, and the other is that both qi and q2 reaches 
their individual maximums. Following the flrst assumption leads to Daganzo's merge 
model (1995), and following the second leads to Lebacque's (1996). 
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In Daganzo's supply- demand method, we have the following optimization prob- 
lem: 

max q = qi + q2 
s.t. 

(3.10) 





< 






< 




qi + q2 


< 






> 


0, 



from which we can find the total flow, 

q = mm{Di + D2,Sd}. 

However, (gi, ^2) may have multiple feasible solutions. This can be shown with Figure 
13.11 when Sd > -D1+-D2, the solution is unique and at point Q; i.e., (gi, ^2) = (-Di, -D2); 
but when Sd < D1+D2, the solution can be any point on the line segment AB. For the 
latter situation, Daganzo defined two (non- negative) distribution fractions ai and 02, 
which satisfy ai + a2 = I and may be related to Di, D2, Sd, and other characteristics 
of the merge. Then, the total flow q is distributed by qi = aiq {i = 1, 2). One example 
when Sd < Di + D2 is depicted in the figure, with given fractions ai and 0:2. Figure 
13.11 also shows that ai or 02 are restricted by Di, D2, and Sd- For instance, for Sd 
given in the figure, ai can not be 1. 

Lebacque suggested another supply-demand method: the supply of the down- 
stream cell is first distributed to the two upstream cells with two fractions ai and 02, 
and it is assumed that the flows qi and q2 reaches their individual maximums, i.e., 
we can compute the flow qi {i = 1, 2) as the following: 

Si = OiiSd, 

(3.11) 

qi = mm{Di,Si}. 
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Figure 3.1: Feasible solutions in Daganzo's supply-demand method 

The feasible solutions of Lebacque's method without fixed fractions are shown in 
Figure As shown, in this model, when Sd > Di + D2, {qi, ^2) can be any point on 
D1BQAD2] when 5*^ < D1 + D2, {qi, ^2) can be any point on D1BAD2. In Lebacque's 
formulation, therefore, ai and a2 are not restricted by Di, D2, or Sd, and the total 
flow q may not reach its maximum minj/^i + D2, Sd} in this method. 

Comparing Daganzo's and Lebacque's methods, we can see that: 1) when the frac- 
tions are the same, the two methods give the same flows; 2) for given Di, D2, and Sd, 
the feasible solution domain of Daganzo's method is contained by that of Lebacque's 
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Figure 3.2: Feasible solutions in Lebacque's supply- demand method 

since the distribution fractions in Daganzo's method (but not in Lebacque's) are 
confined by the supplies and the demand. 

From the above analysis, we can see that both Daganzo's and Lebacque's models 
in the supply-demand framework are based on reasonable assumptions, and Lebacque's 
method Equation 13.111 yields a larger set of feasible solutions than Daganzo's. In 
addition, we think that both formulations are clear and general enough to contain 
physical solutions. Thus, in this chapter, we do not intend to investigate further the 
formulations. Instead, we are interested in the distribution schemes used in these 
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models. 

The reasons why distribution schemes are worth further, deeper discussions are 
as follows. First, we can see from Figure IXD and Figure that distribution schemes 
play a key role in uniquely determining flows through a merge. Thus whether solutions 
of flows are physical is highly dependent on the distribution scheme used. Therefore, 
in order to apply these models to simulate traffic dynamics at a merge, we need a 
better understanding of their distribution schemes. Second, the distribution fractions 
can be affected by travelers' merging behaviors, the geometry of a merge, differences 
between the upstream cells, traffic conditions, and possible control strategies imposed 
on an on-ramp. On the surface, a distribution scheme that models all these factors 
may be extremely complicated. A closer look at various distribution schemes is needed 
to find a simple yet physically meaningful one. Third, it is possible that many valid 
distribution schemes are available. When this happens, a distribution scheme that is 
easy to calibrate and computationally efficient is always preferred. 



3.3 Investigation of various distribution schemes 

In this section, we take a closer look at various distribution schemes and see how their 
distribution fractions are affected by traffic conditions, i.e., Di, D2, and Sd, and other 
characteristics of a merge. We start with a review and a discussion on the existing 
distribution schemes of Daganzo and Lebacque, then propose a simple distribution 
scheme and demonstrate that the supply-demand method incorporating this scheme 
is capable of addressing all factors that we concern about. 
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3.3.1 Discussion of existing distribution schemes 

As we know, different types of links have different characteristics. As a result, even 
when an upstream highway and an on-ramp have the same number of lanes and 
traffic density, the downstream link usually receives more vehicles from the upstream 
highway than from the on-ramp due to differences in design speeds and geometry. For 
example, when vehicles queues up on both a L-lane highway and 1-lane on-ramp that 
merge toget her, the ratio of flow from the on-ramp to th at from th e highw av is about 



1/(2L - 1) (IDaganzo 



19961 1. From these observations, IDaganzol f)1995al ) suggested 
that different upstream links bear different priorities and proposed a distribution 
scheme including parameters for priorities. 

Figure EISl shows how Daganzo's distribution scheme is defined. In the figure, the 
priorities of the highway and the on-ramp are denoted as pi and p2 {pi + P2 = 1); 
respectively. Here the upstream link ui is assumed to have higher priority than U2] 
i.e., pi/p2 > Di/D2- Then the solution (gi,g2) can be shown to be one of three 
cases: i) when Sd < -Di/pi; i.e., 5*^ is x- and y-intercept of line i, the solution 
('?i,'?2) = {PiSd,P2Sd) is at point 1; ii) when Sd G {Di/pi,Di + D2); i.e., Sd is the 
X- and y- intercept of line ii, the solution (^1,^2) = {Di, Sd — -Di) is at point 2; 
iii) when Sd > Di + D2] i.e., Sd is the x- and y- intercept of line iii, the solution 
(?i5?2) = {Di,D2) is at point 3. 

Thus in Daganzo's distribution scheme, we can find that the fraction ai is defined 

as 

Pi, Sd < Di/pi, 



f , Djp, <Sd<D^ + D2, (3.12) 



Sd>D, + D2; 



D1+D2 ■ 

and a2 = 1 — ai. From this definition and Figure ESI "we can see that the priorities pi 
and p2 have to satisfy P2/P1 < D2/D1. i.e., they have to change with respect to traf- 
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Figure 3.3: Solutions of flows in Daganzo's distribution scheme 

fic conditions D2 and Di. Therefore, such "priorities" are not uniquely determined 
by road characteristics, as one would expect. This property makes this distribution 
scheme less attractive. Moreover, even we allow non-constant priorities, the distribu- 
tion scheme with fractions defined in Equation 13. 121 becomes quite complicated when 
considering a merge with more than two upstream links. 

Lebacque (1996) suggested another distribution scheme, in which ai equals to the 
ratio of the number of lanes of link Ui to that of link d. When all branches of a merge 
are highways with the same characteristics, and the traffic conditions of upstream 
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links are overcritical, the demand of each upstream hnk is equal to its lane capacity 
times the number of lanes. In this case, it is reasonable that outflow from each 
upstream link is proportional to the number of lanes; i.e., the distribution scheme 
by Lebacque works well. However, when the upstream links are not similar, e.g., 
one is highway and the other on-ramp, the fractions are obviously not proportional 
to the number of lanes. Lebacque's distribution scheme fails in another case when 
ai + 0:2 > 1 and Si < Di, because it may yield invalid solutions of g = gi + ^2 > Sd- 

3.3.2 A simple distribution scheme and its interpretation 

Our above analysis has revealed certain drawbacks of the existing distribution schemes, 
here we propose a simple distribution scheme, which, as we will see later, removes 
these drawbacks yet is capable of capturing characteristics of a merge. In this dis- 
tribution scheme, the distribution fractions are only related to the demands Di and 
D2, as defined in Equation 13.131 

a, - 



Di+D2- 

As shown in Figure combining the distribution scheme Equation l3.13l with models 
Equation 13.101 or Equation 13.111 we are able to uniquely determine the flows: the 
solution (gi, ^2) with this distribution scheme is simply the intersection of qi + q2 = Sd 
and qi/q2 = D1/D2 when Di + D2 > Sd, and the point Q otherwise. 

A distribution scheme is in fact equivalent to an additional entropy condition, 
which helps identify qi and g2- Thus we also call the distribution scheme Equation 
13.131 the "fairness" condition, because the distribution fractions are proportional to 
traffic demands of upstream links; i.e., the upstream cell with more "sending" flow 
is given more chances. This "fairness" condition, to some extent, is supported by 
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Figure 3.4: Solutions of flows in the simple distribution scheme 

observations at crowded merges, e.g., vehicles from an on-ramp generally wait until 
there is a big enough gap to merge when traffic is fluid. When many vehicles from the 
on-ramp cannot merge and queue up, they may squeeze in and force vehicles from 
the upstrea m mainline freeway to slow down or switch lanes to give way to them 



(Kita 



19991 ). These observations show that vehicles from the upstream cells compete 
"fairly" with each other for admission into the downstream cell. 

From Equation l3.13| we can see that the distribution scheme is not directly related 
to capacities, number of lanes, the difference between upstream links, or control of 
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on-ramps. Thus we can say that this distribution scheme uses the fewest parameters. 
Therefore, it will be easy to calibrate and efficient to compute. Indeed, it is the 
simplest distribution scheme that we can have. Besides, note that the fractions are 
independent of the downstream traffic supply 5*^, although the flows are related to 
Sd- As a mathematical exercise, the following theorem shows the fractions in this 
scheme are in fact the only fractions that are independent of S^- 

Theorem 3.3.1 Suppose that the fractions ai and a2 are independent of the down- 
stream supply Sd, then the fractions will be as in Equation \'J.l'A 

Proof. We have qi = aiq and q2 = ot2q- From Equation 13.111 we then obtain 

«ig < D^, (3.14) 
a2q < D2. (3.15) 

Since ai and a2 are independent of Sd, set Sd = Di + D2, we obtain q = Di + D2- 
Thus, both Equation 13. 141 and Equation 13. 151 have to take the "=" sign, and we have 
Equation 13.131 ■ 
On the other hand, characteristics of a merge can be indirectly captured in the 
simplest distribution scheme Equation 13.131 since capacities, number of lanes, design 
speeds, and on-ramp control can be included in the definition of traffic demands. 
Equation ESI For example, when upstream links have the same per lane capacity and 
are congested. Equation 13.131 will give fractions proportional to the number of lanes. 
Thus this distribution scheme coincides with Lebacque's scheme based on the number 
of lanes. As to Daganzo's "priorities" , they are embedded in the simple distribution 
scheme as follows: when the freeway and the on-ramp have the same number of lanes 
and density, the freeway generally admits higher free flow speed, has higher flow-rate 
and higher demand, and therefore has larger outflow that reflects its higher priority. 
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In addition, the resultant supply-demand method can be ap plied to determine flows 
through a merge with a controlled on-ramp ((Daganzoi . il995ai) : when the metering rate 
of the on-ramp, whose real traffic demand is D2, is r, we can apply the controlled 
traffic demand of the on-ramp min{r, D2} in the supply-demand method. Therefore, 
although the distribution scheme plays "fairly" , the resultant supply-demand method 
of merges can address the characteristics of a merge by incorporating them into the 
computation of traffic demands and traffic supply. 

From discussions above we can see that, in this simple distribution scheme, char- 
acteristics of upstream links, the control of on-ramps, and other properties of a merge 
are captured in the definitions of demand and supply. This is why distribution frac- 
tions depends only on demands in this scheme. Although demand functions of up- 
stream cells are related to many factors and may not be easily obtained, they have 
to be found in all supply-demand methods. In this sense, the simple distribution 
scheme and, therefore, the supply-demand method with this scheme, are the easiest 
to calibrate and the most computationally efficient. 



3.3.3 Properties of the discrete kinematic wave model of 
merges with the simplest distribution scheme 

With the distribution scheme Equation 13.131 the supply-demand method has the 
following further properties: 

Equivalence of models by Daganzo and Lebacque: With the distribution scheme 
Equation 13.131 as shown in Figure EiH the solution (gi,g2) will be on the line 
segment OQ. We can see that Daganzo's model Equation 13.101 and Lebacque's 
modelEquation lH.lll are equivalent with these fractions. Both yield the following 
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fluxes through the merge: 

q = mm{Di + D2,Sd}, 

^1 = ^-15^^ (3-16) 

Extensibility: The supply-demand method incorporating the simplest distribution 
scheme produces qualitatively similar results for merges with different number of 
upstream links. When a merge has U > 2 upstream links, the method Equation 
13.161 can be easily extended as 

q = min{^[i, A,^4, 

Convergence of the merge model: The discrete LWR model is considered as a 
good approximation to the continuous LWR model since it converges as Ax 
while Ax /At is constant. Although analytical convergence analysis of the 
merge model with the simplest distribution scheme has not yet been performed, 
numerical results in section 13.21 do show that it is convergent in the Li norm. 

Consistency of the merge model with the LWR model: Here we conceptually 
consider the consistency of the merge model Equation l3.16l with the LWR model 
for a link with multiple lanes. In the LWR model for a multi-lane link, all lanes 
are assumed to be identical; i.e., given the same initial and boundary condi- 
tions for each lane, flows at the same location on each lane are identical and 
the link's flow-rate or density at a location is simply the number of lanes times 
the flow-rate or density at the location of each lane, respectively, i.e., the LWR 
model does not model lane-changing and traffic in different lanes is treated as 
the same. Therefore a multi-lane link can be considered as an artificial merge: 
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for a boundary inside the link, we separate its upstream part into two links with 
identical flow characteristics, while keep the downstream part intact. Next we 
check if traffic dynamics of this artificial merge is indeed the same as those of 
the original link. Assuming the two upstream links of the artificial merge have 
N and M lanes, respectively, traffic demand of each lane is D, and traffic sup- 
ply of the downstream link is S^- Since the lanes are identical in the upstream 
links, traffic demands for the upstream links are Di = ND, D2 = MD; from 
Equation I3.16| we have 

q = mm{{N + M)D,Sd}, 
N 

M 

Hence, as expected, flow from each upstream lane is min{D, S^/ (N+M)}, which 
is the same as the original boundary flow computed from the LWR model. 

The above analyses indicate that the merge model with the simplest distribu- 
tion scheme is well-defined and qualitatively sound, although the ultimate test of its 
validity rests on empirical validation. 

3.4 Numerical simulations 

In this section, we present our numerical studies of the discrete kinematic wave model 
of merges using the simple distribution scheme. Here we apply Godunov's method 
discussed in Subsection 13. 2. II for each link, and the supply-demand method is used to 
find flows through link boundaries and merging boundaries. In particular, the simple 
distribution scheme is used for computing fluxes through the merge. The resulted 
numerical solution method is described as follows: in each cell. Equation 13.41 is used 
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to update traffic density; we compute flows through link boundaries with Equation 
13 .91 flows through a merge are computed from Equation 13.171 

In the numerical studies, we introduce a unit time r = 5 s and a unit length / = 
0.028 km. Here we study a merge formed by a two-lane mainline freeway and a one- 
lane on-ramp. The three branches of the merge have the same length, L = 400/ = 11.2 
km. The upstream mainline freeway, the on-ramp, and the downstream mainline 
freeway are labeled as links Ui, U2, and d, respectively. The simulation starts from 
t = and ends at t = 500r = 41.7 min. In the following numerical simulations, we 
partition each link into cells and the time interval into K steps, with N/K = 1/10 
always; e.g., if A^ = 50 and K = 500, the cell length is Ax = 81 and the length of 
each time step At = It. 

For both the mainline freeway and the on-ramp, we use the triangular funda- 
mental diagram; i.e., the flow-density relationships are triangular. For the main- 
line freeway, the free flow speed is Vf^m =65 mph=5.1877 1/t; the jam density is 
Pj^m = 2pj=360 veh/km, where pj=180 veh/km/lane is the jam density of a single 
lane; and the critical density pc,m = ^■'^Pj,m = 0.4pj=72 veh/km. Therefore, the 
speed- and flow-density relationships can be written as follows: 



For the on-ramp, the free flow speed Vf^r =35 mph=2.7934 1/t; the jam density is pj; 
and the critical density p^^r = 0.2pj. Similarly, we can have the following speed- and 
flow-density relationships: 
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6 




Figure 3.5: The triangular fundamental diagrams for the mainline freeway and the 
on-ramp 



Qr{p) = pVrip) 

Pj—pc,r 

The above relationships are depicted in Figure 



Vf^rP, < p < Per; 

^/,r(Pi -P)> Pc,r<P<Pj- 



Since |A*,m| < t^/,m = 5.1877Z/r, where A^,,m(p) = Vm{p) + pVI^{p) i s the character- 



istic speed of Equation 13.31 on the mainline freeway, we find the CFL fiCourant et al. 



19281) condition number 



A.,^ -^ < 0.65 <1. 
Ax 
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Since the characteristic speed on the on-ramp is smaller than that on the mainline 
freeway, the CFL condition is also satisfied for the on-ramp. 

3.4.1 Simulation of merging traffic without control 

In this subsection, we study the following merging traffic. Initially, the mainline 
freeway carries a constant flow: traffic densities on the upstream and downstream 
freeway are the same, = = 0.36pj, which is under-critical. After time t = 0, 
a constant flow with density Pu2 = 0.175pj arrives at the on-ramp, and the on-ramp 
remains uncontrolled. In our simulation, the Riemann boundary condition is imposed 
for the upstream boundaries of link ui and U2 and the downstream boundary of link 
d; i.e. the spatial derivatives of traffic density at these boundaries are assumed to be 
zero. 

After partitioning each of the three links into = 500 cells and discretizing the 

time duration into K = 5000 steps, we obtain simulation results as shown in Figure 

13.61 Figure 6(a) illustrates the evolution of traflic density on the freeway upstream 

of the merge: at time t = Or, traffic density is uniformly at pA = 0.36pj; after the 

arrival of the on-ramp flow, freeway traffic immediately upstream of the merging 

point^ becomes congested and reaches a new state ps = 0.7394pj; then a shock wave 

forms and travels upstream in a constant speed Si ~ — 0.61/r = —7.6 mph. Figure 

6(b) shows the evolution of traffic on the freeway downstream of the merge: initially, 

traffic density is also uniformly at pA- After t = 0, traffic immediately downstream of 

the merging point reaches capacity flow at density pc = OApj, and a contact wave^ 

appears since pc and pa are both on the free-flow side of the triangular fundamental 
^Traffic density at the merging point is multiple-valued. 

■^A contact wave is a flow/density discontinuity traveling at the same speed of trafBc on both 
sides of it. 
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diagram. It travels downstream at the speed S2 = 5.21/t = Vf^m- Figure 6(c) shows 
that a backward shock wave also forms on the on-ramp: traffic upstream of the shock 
has density pn = 0.175pj and downstream of the shock density pe = 0.3697pj. This 
shock travels at S3 ~ — 0.251 /r = —3.1 mph. The shock waves and the contact wave 
observed on the three branches are shown in Figure 6(d) on the p — q plane, in which 
line AB represents the shock wave on the freeway upstream of the merge and the 
slope of AB is its speed; line CA represents the contact wave on on the freeway 
downstream of the merge and its travel speed is the slope of CA, which is also the 
free flow speed; finally line DE represents the shock wave on the on-ramp and the 
slope of DE is its travel speed. 

Comparing the congested states on the upstream freeway (point B in Figure 
I3.6r d) ) and that on the on-ramp (point E in Figure l3^ d)). we have the following 
observations: 1) the ratio of the outflows from the upstream branches (freeway: qb = 
1.6349pj7/r=5933 veh/hr and on-ramp: qe = 0.4402pj//r=1597 veh/hr) are not 
proportional to the lane ratio, owing to the different geometric and flow characteristics 
of the two upstream branches, as reflected in their respective fundamental diagrams 
and 2). as long as the freeway downstream of the merge is not congested, and the 
total demand is greater than the total supply at the merging point, traffic states 
surrounding the merge, pB, Pe, and pc, are constant (i.e., stationary) states regardless 
of the initial traffic conditions. This unique characteristics of the merge model with 
our suggested distribution scheme offers a way to validate the model and the fairness 
assumption. 



CHAPTER 3. MERGING TRAFFIC MODEL 



76 



(a) 



5C0t 



4dOT 
I 



3C0t 



2d0x 



1C0t 



Ox 



100 



200 
x/l 



300 



400 




400 



(b) 



0.41 
0.4 
^-0.39 
°^0.38 
0.37 
0.36 



0.35 



Pc 


^ 




25x 

Ox 


1 50x TS^y 

\ V Pa 





400 



500 



2.5 
2 

'1.5 
. 1 
0.5 








0.5 



600 

x/l 

(d) 



700 



800 



c 




A / 










\x Q (p) 













1 

p/p 



1.5 



Figure 3.6: Simulation of merging traffic without control 

3.4.2 Simulation of merging traffic when the on-ramp is con- 
trolled 

In this subsection, we have the same initial/boundary conditions as in the preceding 
subsection, but with the on-ramp controlled by a ramp meter. For simplicity we take 
a constant metering rate r = 0.3445pj//r=1250 veh/hr. The simulation results are 
shown in Figure IT71 From Figure ITTT a). we can see that a backward traveling shock 
wave forms on the freeway upstream of the merge, traveling at si ^ — 0.33//r = —4.1 
mph. Traffic densities besides the shock are pA = 0.36pj (upstream) and pb = 
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Figure 3.7: Simulation of merging traffic with on-ramp control 

0.6278pj (downstream), respectively. Figure ITTT b) shows the evolution of traffic on 
the freeway downstream of the merge, which is identical to the case without ramp 
control. From Figure ITTT c). we can see that a backward traveling shock wave also 
forms on the on-ramp, traveling at S3 ~ — 0.48//r = —6.0 mph. On the ramp, 
traffic densities besides the shock are po = 0.175pj (upstream) and pe = 0.577pj 
(downstream), respectively. Again, Figure ITTT d) shows the the initial and congested 
states in the p — q plane. 

We note that with ramp control, the freeway upstream of the merge becomes less 
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congested (pb is lower) while the on-ramp becomes more congested {pE is higher). 
Furthermore, the freeway queue grows slower while the ramp queue grows faster with 
ramp control than without. Finally, the ramp control also affects the stationary states 
{PBi PCi Pe) and the distribution fractions, in favor of discharging higher flow from 
the freeway. 

3.4.3 Computation of convergence rates 

In this subsection we will check the convergence of the merge model with the distri- 
bution scheme Equation 13. 131 when the on-ramp is not controlled. Here, this is done 
by computing convergence rates of traffic density over the whole network. 

First, we compare traffic density solutions at time Tq = 500r for two different 
number of cells into which the network is partitioned and obtain their difference. 
Denote solutions as ([//^)f^]^ for 2N cells and {U[^)fLi for cells respectively, and 
define a difference vector {e"^^'^)^-^ between these two solutions as 

ef-^ = liUl^-, + Ul'')-Ur,^ = l,■■■,N. (3.18) 
Then, the relative error with respect to L^-, L"^- or L°°-norm can be computed as 

^2N-N ^ ||e2^-^||. (3.19) 
Finally, a convergence rate is obtained when we compare two relative errors: 

^2N-N 

Here the vector of U contains the densities of all three links which are weighted by 
the number of lanes of each link. 

We will compute convergence rate with the following conditions. For link ui, the 
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number of lanes a{ui) = 2, we define its initial condition as 

p{x,0) = a(Mi)(0.18 + 0.05 sin f^)pj, xe[0,L], 
v{x,0) = Vm{p{x,0)), xe[0,L]. 
For link U2, the number of lanes a{u2) = 1, we define its initial condition as 

p(x,0) = a(u2)(0.175 + 0.05 sin ^)pj, xG[0,L], 

^3.22^ 

v{x,0) = Vr{p{x,0)), xe[0,L]. 
For link d, the number of lanes a{d) = 2, we define its initial condition as 

p(x,0) = a(c/)(0.18- 0.05 sin ^)pj-, xe[L,2L], 
v{x,0) = VMx,0)), xe[L,2L]. 
In addition, we impose the Riemann boundary condition for the upstream boundaries 
of link ui and U2 and the downstream boundary of link d; i.e., p has zero derivative 
at these boundaries. 

From Table im convergence rates of the merge model with the simple distribution 
scheme are of order one in norm. The convergence rate is of order one is expected 
because the Godunov method used here is a first-order method, in which traffic density 
on a link is approximated by piece-wise constant functions. Unlike the Godunov 
discretization of the LWR model, which we know is convergent to the LWR model, 
we do not yet know what differential equations the discrete model converges to. 
Nevertheless, the numerical convergence analysis gives us comfort in applying the 
discrete merge model in the sense that we know the properties of the computed 
solutions will not change over discretization scales. 



3.5 Discussions 



In this chapter, we studied the discrete kinematic wave model of merges in the 
supply-demand framework, probed the supply-demand merge models by Daganzo 
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Table 3.1: Convergence rates of the discrete merge model 

and Lebacque, and gained a better understanding of various distribution schemes 
in merge models. More importantly, wc proposed the simplest distribution scheme, 
equivalent to the so-called "fairness" condition, in which the distribution fractions 
are proportional to the upstream demands. We demonstrated with both analytical 
discussions and numerical studies that the discrete merge model with the simplest 
distribution scheme is well-defined. 

The simple distribution scheme can also be easily extended to diverges in single- 
commodity traffic flow. In reality, a diverge is not simply a reverse of a merge: in 
a merge, vehicles advance in the same direction and once they travel into the down- 
stream link, their origins no longer affect traffic dynamics; while in a diverge, traffic 
dynamics in the upstream link are affected by the combination of vehicles with differ- 
ent destinations. However, for single-commodity traffic flows, in which vehicles have 
no predefined routes and can choose any downstream links when diverging, we can 
have similar supply-demand method for computing fluxes through a diverge: assum- 
ing that a diverge has K downstream links, the upstream link and the downstream 
links are properly discretized, the traffic demand of the upstream link is D^, and traf- 
fic supply of downstream link di {i = 1, K) is Si, the supply-demand method with 
the "fairness" condition, in which the influx into a downstream link is proportional 
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to its accommodation St, can be written as: 



q = mm{D^,Y,. S^}, 
Qi = PiQ, (3-24) 

where Qi is influx into the downstream link di. 

Further, for a general junction with J upstream links and K downstream links, 
we can compute fluxes as 

K 



qt = f3kq, (3.25) 
a, = J = 1,---,J 

where Dj is the traffic demand of the upstream link Uj, Sk is the traffic supply of the 
downstream link d^, is the out-flux from Uj, and qf the influx into d^. Although 
this model may not be proper for realistic traffic flow, it could be helpful for the study 
of other unidirectional flows. 

With a complete picture of the merge model, we can now simulate how a traffic 
disturbance propagates through a merge. In particul ar, since merges are well rec- 



1999b!), this study 



ognized as locations where congestion often initiates (jDaganzo, 
will help traffic engineers understand how congestion start and propagate at a merge. 
Moreover, with the complete discrete kinematic wave merge model, we are able to 
develop and evaluat e loca l ramp metering strategies, e.g., the ALINEA algorithm 



flPapageorgiou et a). 



19971 ). In the future, we will be also interested in the limit of 



the discrete kinematic wave model with the simple distribution scheme. In particu- 
lar, we will be interested in kinematic wave solutions for a general Riemann problem 
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of this model. Finally, and more importantly, the supply-demand method and the 
"fairness" condition needs to be checked against observed merging phenomena. 



Chapter 4 

Kinematic wave traffic flow model 
of diverging traffic 



4.1 Introduction 



Diverges have b een w ell recognized as a major type of bottlenecks (jDaganzo, 



Daganzo et al. 



1999b; 



19991 1 ■ where congestion of one downstream branch can propagate 
to the upstream branch and fu rther blocks vehicles traveling to other directions. 



Observations of diverging traffic (|Munoz and Daganzo 



2002, 



20031 ) have shown that 



vehicles to different downstream branches tend to segregate themselves at a point, 
called the actual diverging point, as far as 2 km upstream of the physical diverging 
point. In addition, traffic upstream of the actual diverging point is in the so-called 
"1-pipe" regime in the sense that vehicles that appear at the same time and location 
have the same speed even they have different destinations; while between the actual 
and physical diverging points traffic is in a "2-pipe" regime, since vehicles to different 
directions are separated and travel independently. I n certain case s, however, the 



actual and physical diverging points may be the same ([Newel], 



19991). Although it is 
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an interesting problem to determine the location of an actual diverging point, here 
we only intend to study traffic dynamics at an actual diverging point or simply a 
highway diverge; i.e., we want to know, under certain traffic conditions, how much 
traffic will flow into each downstream branch. 

Due to its simplicity in repr esentation and efficiency in computation , the k ine- 



matic wave theory developed by fjLighthill and Whitham . 



1955b 



Richards 



19561 ) has 



played an important role in understanding traffic dynamics on road networks. In the 
framework of the kinematic wave theory of diverging traffic, there have been two ma- 
jor lines of research in the framework of the kinematic wave theory. In the ffist line 
of research, the effect of divergin g traffic to the mainline freewav is incorporated by 



considering an off-ramp as a sink (jKuhne and Michalopoulo; 



19921 ). In these models. 



only part of the interaction between the diverging traffic and the main highway traffic, 
i.e., the impact of diverging traffic to the mainline traffic, is taken into account. In ad- 
dition, the characteristics of sinks are not trivial to obtain. In contrast, in the second 
line of research, the complete dynamics at diverges is studied by treating each branch 
equally. When vehicles do not have predefined routes and can choose any down- 



stream branches at a diverge, 



Holden and Risebro 



branch by solving an optimization problem, and 



19951) determined the flow to each 



Jin and Zhand (|2003d ) proposed a 



distribution scheme based on the deflnitions of traffic supply a nd demand. For more 



realistic situatio ns when vehicles have predeflned route choices. 



Daganzd ()l995ar ) and 



Lebacqud (|l996r i proposed the so-called supply-demand method for computing traffic 



flows through a diverge. In this method, each branch is partitioned into a number of 
cells and a duration of time is discretized into time steps. Then a diverge connects 
an upstream cell and a number of downstream cells, and the flow to a downstream 
cell during a time period is assumed to be proportional to the fraction of vehicles in 
the upstream cell that wish to travel to this downstream cell. In addition, the flows 
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to the downstream cells are assumed to reach their maxima subject to the constraint 
of traffic demand of the upstream cell and supplies of the downstream cells. 

Although the supply-demand method is quite straightforward in computing flows 
through a diverge, it fails to provide a system-wide picture of how a small traffic 
perturbation on one branch will propagate through the diverge. That is, it gives 
no answer to what kind of kinematic waves will initiate in diverging traffic. In this 
chapter, we intend to study the kinematic waves (more specifically, the instantaneous 
kinematic waves) at a diverge. In our study, we still partition each branch into a 
number of cells and discretize a duration of time into some time steps. Then in a 
(short) time step^, we compute the in-flow to a downstream branch by holding traffic 
conditions of vehicles to the other destinations constant in the upstream cell and, after 
obtaining all flows, update traffic conditions in all the upstream and downstream cells. 
We can see that, in this theory, dynamics of traffic traveling to different destinations 
are ffist decoupled during a short time period and then combined together in an 
alternating manner. With this decoupling, corresponding to each downstream branch, 
we can obtain a hyperbolic system of conservation laws, which is independent of 
the systems corresponding to the other downstream branches. To the best of our 
knowledge, this decoupling approach in solving complicated dynamics of diverging 
traffic is the ffist of its kind in traffic flow. 

Moreover, we w ill show that the decoupled systems are nonlinear resonant sys- 



tems in the sense of 



Isaacson and Templd f)l992f l and can be solved by combinations 



of shock, rarefaction, and standing waves. Here, these waves are considered instan- 
taneous since they only exist in a short time interval, and waves corresponding to 
different downstream branches interact with each other at the end of each time inter- 



^The length of this "short" time step At, theoreticaUy, can be infinitely small and, practically, 
can be obtained through a calibration exercise. 
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val. Thus, after a long time, the asymptotic kinematic waves arising from diverging 
traffic can be totally different from the instantaneous waves. In this chapter, we will 
focus on the instantaneous waves arising from each branch, and the interaction be- 
tween these waves and the asjTiiptotic kinematic waves will be further studied in a 
following chapter. 

The rest of this chapter is organized as follows. In Section 14.21 we first formulate 
the kinematic wave theory of diverges in the framework of the LWR model, and 
introduce the instantaneous kinematic wave theory for dynamics of diverging traffic. 
In Section r4.3[ we solve in detail for the instantaneous kinematic waves. In Section 
14. 4| we present a new definition of traffic demand for diverging vehicles and the 
supply-demand method based on this new definition. In Section 14.51 we carry out 
simulations of the new model for two cases. In Section 14. 6[ we make our conclusions 
and propose some follow-up research. 

4.2 A kinematic wave theory for diverges 

As we know, in a traffic stream, traffic dynamics can be affected by characteristics such 
as vehicle types and destinations, and drivers' behavior. Here vehicles of the same 
characteristics are considered to belong to a commodity. In this study we are only 
interested in the impact of vehicles' destinations on traffic dynamics at a highway 
diverge. Therefore, for a diverge with D downstream branches and one upstream 
branch, we can differentiate vehicles into D commodities. 

Here we consider an ideal diverge: on the upstream branch, traffic is considered 
in 1-pipe regime; once through the diverge, vehicles completely segregate. Therefore, 
we have a D-commodity flow on the upstream link and a single- commodity flow on 
each downstream link, and traffic dynamics on a road network with a diverge consists 
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of three parts: dynamics of the single-commodity flow on each downstream branch, 
dynamics of the D-commodity flow on the upstream branch, and dynamics at the 
diverge. In the following part, after reviewing the kinematic wave models of single- 
commodity and D-commodity traffic, we will focus on the study of traffic dynamics 
at the diverge. 



4.2.1 The kinematic wave theory of single-commodity traffic 
flow 



For dynamics of single - commodity traffic, the seminal LWR (jLighthill and Whitham 



1955b; 



Richards 



19561) kinematic wave theory can be applied. In the LWR model for 
traffic flow on a long crowded road, two equations are used to describe traffic dynamics 
of aggregate quantities such as traffic density p, flow-rate q and travel speed v: 

Traffic conservation: Pt + Qx = 0, (4.1) 
Fundamental diagram: q = Q{a,p). (4.2) 

In Equation 14.21 a{x) is an inhomogeneity factor depending on road characteristics 
such as the number of lanes at x. Since q = pv, we also have a speed-density relation: 
V = V{a,p) = Q{a,p)/p; for vehicular traffic, generally, v is nonincreasing and the 
equilibrium flow-density relation is assumed to be concave. From Equation 14.11 and 
Equation 14. 2t the LWR model can be written as 

Pt + Q{a,p)x = 0. (4.3) 

When the road is homogeneous; i.e., a{x) is constant, the LWR model is called ho- 
mogeneous; otherwise, when a{x) is dependent on location. Equation 14.31 is called an 
inhomogeneous LWR model. 
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The homogeneous LWR model is nothing but a scalar conservation law and has 
been well studied both analytically and numerically. It is solved by combinations 
of shock and expansion (rarefaction) waves. In contrast, the inhomogeneous LWR 
model Equation 14.31 can be considered as a nonlinear resonant system and has an 
additional type of kinematic waves - s tanding waves, which initiate and stay at the 



inhomogeneity (jjin and Zhang . 



2003br ). Therefore, for a single-commodity flow on a 
road branch, traffic dynamics can be considered as the evolution of a combination of 
shock, rarefaction, and standing waves. 



4.2.2 The kinematic wave theory of multi-commodity traffic 
flow 

For traffic on the upstream link, it is assumed that vehicles with different destinations 
have no differences in dynamics; i.e., the speed-density relation is independent of 
commodities. Let pi{x,t), Vi{x,t), and qi{x,t) {i = 1, ■ ■ ■ , D) denote, at place x and 
time t, the density, speed, and flow-rate of commodity i, respectively. In contrast, 
p{x,t), v{x,t), and q{x,t) are the aggregate density, speed, and flow-rate. Then we 
have (with {x, t) suppressed hereafter) 

D 

P = Y.P^^ (4.4) 

i=l 

V = Vi = V{a,p), i = l,---,D, (4.5) 

D 

q = (4.6) 

1=1 

where q = pv, qi = piVi {i = 1, - ■ ■ ,D), V{a, p) is the aggregate speed-density relation. 
We also call this kind of mult i- commodity traffic as additive. In additive traffic, the 
flow of each commodity is called a partial flow and the flow containing all commodities 
the total flow. 
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Therefore, in additive traffic, dynamics of each partial flow can be written as 

{p^)t + {p^V{a,p)).^ = 0, i = l,---,D. (4.7) 

Note that these D equations are coupled through the aggregate speed-density relation. 
By introduc ing the fraction of vehicles of commodity i, S,i = Pi/p, which satisfies 



be written as 



Lebacqud (|l996[ ) showed that traffic dynamics of additive traffic can also 



pt + ipV{a,p))^ 



0, 
0, 



i = 1,---,D-1. 



(4.8) 
(4.9) 



We can see that the fractions travel in the same speed as vehicles. This shows that the 
First-I n-First-Out princ iple is guaranteed in the kinematic wave theory of additive 



traffic (jLebacque . 



19961 ). In addition, we can see that, for multi-commodity traffic, 
the aggregate traffic dynamics can be described by the (inhomogeneous) LWR model 
and is still a combination of shock waves, rarefaction waves, and/or standing waves, 
and traffic dynamics of a commodity is of the similar pattern. 



4.2.3 A kinematic wave theory of diverging traffic 

To model traffic dynamics on a network with a diverge, the harder problem is to 
find the number of vehicles traveling from the upstream link to each downstream link 
during a time interval. This is an important problem to answer since flows from the 
upstream link to the downstream links are the necessary boundary conditions for the 
connected links. Here we propose a new approach to finding the approximate flows 



Daeanzo , 


1995a: 


Lebacque . 


1996) 



19961 ) in that this theory can be analytically solved by a 



combination of kinematic waves. 
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Like the supply-demand method of (jDaganzo, 



1995a : 



Lebacque . 



19961 ). the new 



kinematic wave theory is also in a discrete form: we partition each branch into a 
number of cells of length Ax, discretize the time duration into time steps of length 
At, and assume that the aggregate and partial traffic densities are constant in each 
cell at a time step. The discr etization of space-tim e plane in this model is required 
to satisfy the CFL condition ffCourant et al.L 119281) . which guarantees that the flows 
through the diverge are only dependent on traffic densities in the adjacent cells and 
are independent of traffic conditions farther upstream or downstream. 

In the discrete model, connected to the diverge with D downstream branches are 
an upstream cell and D downstream cells, denoted as downstream cell i {i = 1, ■ ■ ■ , D). 
Assume that, at t = 0, partial densities in the upstream cell are pf and the aggregate 
density p'" = J^iLi pf^ density in downstream cell i is pf . With these initial traffic 
conditions, we will show how traffic dynamics will evolve from t = till t = At. In 
addition, we will be able to compute the in-fiow to downstream cell i {i = 1, ■ ■ ■ , D), 
Qi. Then from traffic conservation, the out-fiow from the upstream cell is YliLili- 
With these flows, traffic conditions at the next time step t = At can be obtained. 

In this new model, to study the traffic dynamics of commodity i through the 
diverge from t = to t = At, we hold partial densities of vehicles of the other 
commodities in the upstream cell constant during this time interval. That is, pj", 
where j = 1, ■ ■ ■ , D and j ^ i, are all constant. Here we assume that traffic densities of 
commodities other than i are constant during the short time interval. This assumption 
is equivalent to saying that the change of traffic conditions of commodity i is, in a 
short time interval, independent of traffic dynamics in the other directions. This 
assumption appears to be a natural and reasonable way to decouple the interactions 
between traffic in different directions. This approach can be considered as an attempt 
to answer an important question: how do vehicles of different destinations segregate 
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themselves at a diverge? 

We set the diverging point at x = and vehicles travel in the positive direction 
of X-axis. Thus, x < for the upstream link and x > for the downstream link. 
Also, to simplify the notations, we rewrite pi as p and the total traffic densities of 
other commodities as k. Therefore, traffic dynamics of vehicles of commodity i in 
the upstream cell can be written as pt + {pV{a{x)^ p + k))^ = for x < 0, and 
pt + (p ^ {.o.i.^) 1 P)) X = for X > 0. Since k is held constant during the time interval 
[0, At], and /c = in downstream cell we have the following equation and conditions 
for traffic dynamics of commodity i: 

p, + {pV{a{x),p + k))^ = 0, (4.10) 

with the following initial traffic conditions 

{PL X <0, 
t = 0, 
Pr X > 0; 
(4.11) 

ki X < 0, 

k = I te [0, At]. 

[ X > 0; 

Here we extend our x-axis to both positive and negative infinity since, with the CFL 
condition, the traffic dynamics only depend on the conditions above. 

Thus, Equation 14.101 with Equation 14.111 is a description of traffic dynamics of 
commodity i during a short time interval. This new model will be shown to be solved 
by instantaneous kinematic waves in the following section. 



4.3 The instantaneous kinematic waves 



In this section, we discuss the instantaneous kinematic waves of Equation 14.101 with 
initial conditions Equation 14.111 For the purpose of exposition, we here assume that 
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both the upstream and downstream branches share the same speed-density relation; 
i.e., we exclude the road characteristics a{x). We also assume that q is concave and 
V non-increasing in p G [0,Pj], where pj is the jam density. With modifications, the 
results here are applicable to the case when the upstream and downstream branches 
have different characteristics. 

Since we assume that k is time independent, we have another conservation equa- 
tion: kt = 0. Therefore Equation 14. lUI with Equation 14. Ill can be rewritten cis ci 2 X 2 
system: 

p, + {pV{p + k)), = 0, ^^^^^ 
kt = 0, 

together with initial conditions 

p{x,t = 0) = 

k{x,t = 0) = 



PL 


X 


<0, 


PR 


X 


>0, 


kL 


X 


<0, 





X 


> 0. 



(4.13) 



Next, we will show that Equation 14.121 is indeed a nonlinear resonant system 



studied bv llsaacson and Templd (|l992l ). Consequently, we can use their technique to 
solve this system, which describes traffic dynamics of one commodity during a short 
time period. Note that here we do not restrict our time domain to be [0,ti]. This 
is because Equation 14.121 with Equation 14.131 has self-similar solutions, which means 
that the ffow f{x = 0, t) is constant for any time t > 0. 
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4.3.1 The properties of Equation 14. 121 as a nonlinear resonant 
system 

Denote the state U = {k,p), the flux vector F{U) = {0, pV{p + k)), Equation 14.121 
can be written as a hyperbolic system of conservation laws: 



Ut + FiU), = 0, 



(4.14) 



where x G (— oo, oo), t > 0. This system can be written in the following form of a 
quasi-linear system 



Ut + dF{U)U, = 0, 
where the Jacobian matrix dF{U) is 



(4.15) 



dF 





pV'ip + k) V{p + k)+ pV'ip + k) 

Thus the two eigenvalues of dF{U) are 



(4.16) 



Ao = 0, \, = V{p + k)+pV\p + k), 

and their corresponding right eigenvectors are 

V{p + k)+ pV'{p + k) 
-pV'ip + k) 

System Equation 14.141 is a non-strictly hyperbolic system, since it may happen 



— * 











1 



that 



Ai(f/,) = Ao(f/*) = 0. 



(4.17) 



When Equation 14.171 is satisfied, we say that traffic state U^, is partially critical. If 
denoting the partial fundamental diagram as Q{p;k) = pV{p + k), we can define 
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other partial quantities of the commodity that we are interested in, partial capacity 
Qmaxj-^^ and partial critical density 7(A;), as follows: 

Q{p;k) < g^"^^(A;), \fpe[0,pj-k], (4.18) 
Q--(A;) = Qi^iky,k). (4.19) 

Since l^' < and (f{pV{p))/dp'^ < 0, we have 

d^Q{p- k) 



dp^ 
dQjp; k) 
dk 



< 0, 

< 0; 



therefore, Q"^^^{k) and 7(/c) are unique for given k G [0, pj]. Thus the partially critical 
state f/* = {k^,p^) = {k^,'y{k^)). We call the collection of the partially critical states, 
r = {U^\k G [0,rj]}, the transition curve. 

Moreover, at a partially critical state t/*. Equation 14.121 is genuinely nonlinear; 



^Ai(t/.) = ^{pV{p + k))<0, (4.20) 

since q = pV{p + k) is concave in p, and dF{U) is nondegenerate; i.e., 
d 

— {p,V{p, + k,)) = p,V'{p, + k,) = -V{p, + k,) <0. (4.21) 



With conditions Equation 14. 171 - 



ear resonant system in the sense of (jlsaacson and Temple , 



auation [4.211 satisfied, Equ ation I^TT^ is a nonlin 



19921 ) , and it is guaranteed 



that, in a neighborhood of the state U^,, the Riemann problem for Equation 14. 121 with 
Eq nation 14 . 1 31 has a unique solution with a canonical structure. 
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4.3.2 The instantaneous kinematic waves of Equation 14.121 
with Equation 14.131 

The Riemann problem for Equation 14.121 with initial conditions Equation 14.131 is 
solved by two families of basic waves, associated with the two eigenvalues. The 0- 
waves with wave speed Aq = 0, also called standing waves, are the integral curves of 
Rq in ?7-space, and therefore are given by Q{p] fc)=const. Similarly, the 1-waves with 
wave speed Xi{U) are integral curves of -Ri in [/-space, and are given by A;=const. 
As shown in Figure 14.11 the 0-curve is concave and tangent to the 1-curve at the 
critical state f/*; states left to the transition wave r(A;) are undercritical (UC), and 
ove rcritical (OC) right to the tr ansition wave. 



Isaacson and Templd ()l992r ) showed that wave solutions to the Riemann problem 
are formed by no more than three basic waves, including standing waves, shock waves, 
and rarefaction waves, and these waves satisfy two entropy conditions: Lax's entropy 
condition (1972), i.e., that these waves increase their wave speeds from left (upstream) 
to right (downstream), and an additional condition that an UC state and an OC state 
can not be connected by a standing wave. With these conditions, wave solutions exist 
and are unique: when Ul is UC, wave solutions in ?7-space are shown in Figure 14. 2[ 
and in Figure when Ul is OC. Note that the solutions can also be categorized 
according to whether pji is OC or UC. Also recall that Uji always lies on A; = 0. 

In the remaining part of this subsection, we will discuss the Riemann solutions 
of Equation 14. 121 with Equation I4.13| and compute the boundary flux q{x = 0,t). 

1. When Ul is UC; i.e., pl < l{kL), we denote the intersection between Q{Ul) 
and = by Pl and p2, where pi < 7(0) < p2- Hence, we obtain three types of 
solutions when pR G [0,pi], (pi,P2), or [p2,Pj] as shown in Figure Hill 

Type 1 When p/j G [0,pi]; i.e., Q{pr;0) < Q{Ul) and Ur is UC, wave solutions 
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Figure 4.1: Integral curves in (p, fc)-space 

to the Riemann problem consist of two basic waves with the intermediate 
state, Ui = (pi, 0): the left wave {Ui, Ui) is a standing wave and the right 
wave (f/i, Ur) is a rarefaction wave with characteristic velocity Xi{U) > 0. 
As shown in Figure lOf the boundary flux q{x = 0,t > 0) = Q{Ul)- 

Type 2 When G (pi,p2); i.e., Q{pb,;0) > Q{Ul), wave solutions to the Rie- 
mann problem consist of two basic waves with the intermediate state, 
Ui = (pi, 0): the left wave {Ul, Ui) is a standing wave and the right wave 
{Ui, Ur) is a shock wave with positive wave speed s{Ui, Ur) = {Q{pr; 0) — 



CHAPTER 4. DIVERGING TRAFFIC MODEL 



97 




Figure 4.2: The Riemann solutions when Ul is \JC 

Q{Ul))/ {pr — Pi)- As shown in Figure the boundary flux q{x = 0,t > 
0) = QiUL). 

Type 3 When p/j G [p2,Pj]; i.e., Q{pr;0) < Q{Ul) and fZ/j is OC, wave solutions 
to the Riemann problem consist of two basic waves with the intermediate 
state, Ui = (pi. Pi) and Q{Ui) = Q{pr] 0): the left wave {Ul, Ui) is a shock 
wave with negative wave speed s(f/L, t/i) = {Q{Ul) — Q{Ui))/{pL — pi) and 
the right wave {Ui,Ur) is a standing wave. As shown in Figure the 
boundary flux q{x = 0,t > 0) = Q{Ul)- 



CHAPTER 4. DIVERGING TRAFFIC MODEL 



98 




Figure 4.3: The Riemann solutions when Ul is OC 

2. When Ul is OC; i.e., pL > 7(^l), we define [/* by f/* G T{k) and k^, = ki- 
Therefore, Q{U^) is the partial capacity Q"^^^{kL)- Denoting the intersection 
between Q{p;k) = Q{Ul) and (p > 7(0), /c = 0) by po and the intersections 
between Q{p; k) = Q{U^) and = by pi and p2, where pi < 7(0) < p2, we 
can obtain four types of solutions when p^ G [0,pi], (pi,p2), [p2,Po], or (po,Pj], 
as shown in Figure lOl 

Type 4 When p^ G [0,pi]; i.e., Q(p_r;0) < Q{U^) and Ur is UC, wave solutions 
to the Riemann problem consist of three basic waves with two intermedi- 
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Figure 4.4: An example for wave solutions of type 1 for Equation 14.121 with initial 
conditions Equation 14.131 

ate states, and Ui, where Q{Ui) = Q{U^) and ki = 0: the left wave 
{Ul, U^) is a rarefaction wave with non-positive characteristic velocity, the 
center wave t/i) is a standing wave, and the right wave {Ui^Un) is a 
rarefaction wave with positive characteristic velocity. As shown in Figure 
Wl\ the boundary flux q{x = 0,t > 0) = Q{U^) = Q'^^''{kL). 

Type 5 When G (pi,p2); i-e., Q{pr] 0) > Q{U^), wave solutions to the Riemann 
problem consist of three basic waves with two intermediate states, and 
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t 




p X 



Figure 4.5: An example for wave solutions of type 2 for Equation 14.121 with initial 
conditions Equation 14.131 

Ui, where Q{Ui) = Q{U^) and ki = 0: the left wave {Ul, U^) is a rarefaction 
wave with non-positive characteristic velocity, the center wave (f/*,t/i) is 
a standing wave, and the right wave (t/i, Un) is a shock wave with positive 
speed. As shown in Figure 14.81 the boundary flux q{x = 0, t > 0) = 

Type 6 When pp. e [p2,Po]; i-e., Q{Ul) < Q{pr;0) < Q{U^) and Ur is OC, wave 
solutions to the Riemann problem consist of two basic waves with the 
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Figure 4.6: An example for wave solutions of type 3 for Equation 14.121 with initial 
conditions Equation 14.131 

intermediate state Ui, where Q{Ui) = Q{Ur) and ki = ki. the left wave 
{Ul,Ui) is a rarefaction wave with negative characteristic velocity and 
the right wave {Ui,Ur) is a standing wave. As shown in Figure the 
boundary flux q{x = 0, t > 0) = Q{Ur). 

Type 7 When G (po,Pj]; i-e., Q{pr;0) < Q{Ul) and Ur is OC, wave solutions 
to the Riemann problem consist of two basic waves with the intermediate 
state f/i, where Q{Ui) = Q{Ur) and ki = ki- the left wave {Ul,Ui) is a 



CHAPTER 4. DIVERGING TRAFFIC MODEL 



102 




Figure 4.7: An example for wave solutions of type 4 for Equation 14.121 with initial 
conditions Equation 14.131 

shock wave with negative speed and the right wave (f/i, Ur) is a standing 
wave. As shown in Figure I^!TU1 the boundary flux q{x = 0,t > 0) = Q{Ur). 
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t=t„ 












u. 





t=t„ 



Figure 4.8: An example for wave solutions of type 5 for Equation 14.121 with initial 
conditions Equation 14.131 

4.4 The supply-demand method with a new defi- 
nition of traffic demand 



Based on the discussions in the previous section, we summarize solutions of boundary 
flux q{x = 0,t > 0) in Table lO 

Further, if we introduce a new definition of partial traffic demand for commodity 
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Figure 4.9: An example for wave solutions of type 6 for Equation 14.121 with initial 
conditions Equation 14.131 



1 in the upstream link as 

' Q{Vl) Ul is UC 
Qmaxj^^^^ otherwise 

the boundary flux q{x = 0,t > 0) can then be computed by 



(4.22) 



q{x = 0,t>0) = mm{Si,D{pL;kL)}, 
where the supply of the downstream link, Si, is the same as in (jDaganzo 



(4.23) 



1995a; 
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Figure 4.10: An example for wave solutions of type 7 for Equation 14. 121 with initial 
conditions Equation 14.131 



Lebacqud . Il996r i: i.e. 



Si 



(4.24) 



Q{Uu) Ur is OC 
gmaxj^g^ Otherwise 

The supply-demand method Equation 14.231 can also be used to calculate the 
boundary fluxes of other commodities. It yields the same solutions to the Riemann 
problem of Equation 14.121 with Equation 14.131 as the analytical solution method. 
Moreover, it is much simpler in computation and can be easily extended to the more 
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Solution type 


left state Ul 


right state pn 


q{x = 0,t > 0) 


1 


UC 


Q{pr;0)<Q{Ul),pr<i{0) 


Q{Ul) 


2 


UC 


Q{Pr;0) > Q{Ul) 


Q{Ul) 


3 


UC 


Q{Pr;0)<Q{Ul),Pb.>i{0) 


Q{Ur) 


4 


OC 


Q{pr;0)<Q{Ul),pr>i{0) 




5 


OC 


JXUr) < f^"'', an > ai, pR/aR < a 




6 


OC 


f{UR) > fr^ 


Q{Ur) 


7 


OC 


/{Ur) < fiUi), pR/aR <a, aR< ai 


Q{Ur) 



Table 4.1: Solutions of the boundary flux q{x = 0,t > 0) 

complicated cases when the upstream and downstream branches have different road 
characteristics. 

4.5 Numerical simulations 



In this section, we carry out numerical simulations of the instantaneous kinematic 
wave model presented in this chapter. We will study a small diverging network 
consisting of two downstream links and one upstream link: the lengths of three links 
are the same, L = 400/ = 11.2 km, with unit length / = 0.028 km; one downstream 
link, labeled as link di, has a{di) = 2 lanes, another downstream link, link d2, has 
a{d2) = 1 lane, and the upstream link, link u, has a{u) = 2 lanes. The simulation 
starts from t = and ends at t = 500r = 41.7 min, with unit time r = 5 s. In the 
following simulations, we partition each link into cells and the time interval into 
K steps, with N/K = 1/10 always; e.g., if A^ = 50 and K = 500, the cell length is 
Ax = 81 km and the length of each time step At = Ir s. 



We will use the exponential fundamental diagram (|Del Castillo and Benitez 



1995bh . 
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Figure 4.11: The Newell model of speed-density and flow-density relations 



and the parameters are given as follows: the free flow speed Vf = 5.01 /r = 0.028 km/s 
= 100.8 km/h; the jam density of a single lane pj = 180 veh/km/lane; the wave veloc- 
ity for jam traffic Cj = — l.OZ/r = -0.0056 km/s = -20.16 km/h; and the equilibrium 
speed-density relationship 



V{a,p) 



5 p 



1/t, 



where a{x) is the number of lanes at location x. The equilibrium functions V{a, p) and 
Q(a, p) are shown in Figure d. Ill in which the critical traffic densitv is a = 0.259apj. 



Here we apply the first-order Godunov method (IGodunov . 



19591), in which the 
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traffic conditions are updated for each cell based on traffic conservation, and boundary 
fluxes are computed with the supply-demand method described in Section 4. Since 
^* ^ Vf = 5//r, we find the CFL condition number 

< 0.625 <1. 

Therefore, the first-order Godunov method can solve the instantaneous kinematic 
wave model efficiently. 



4.5.1 Simulation I: A general case 

In this simulation, we study a general case of diverging traffic. Initially, the upstream 
branch carries a constant fiow with traffic density pu = l.llllpj = 200 veh/km and 
the proportion of vehicles traveling to downstream branch di is 80%; downstream link 
di is empty; traffic density on downstream link d2 is = 0.5556pj = 100 veh/km. 
In addition, we impose the Neumann boundary conditions on the boundary of this 
diverging network; i.e., the spatial derivatives of traffic densities at the boundary are 
set to be zero. 

With = 500 and K = 5000, we obtain simulation results as shown in Figure 
14.121 Fig. I4.12f a) illustrates the evolution of traffic on the upstream link: at time 
t = Or, traffic density is uniformly at pa = l-Hpj; after the beginning of diverging 
process, traffic immediately upstream of the diverging point reaches a new state pb = 
0.69pj, which keeps propagating on the upstream link; as a result, an expansion wave 
forms and travels upstream. Fig ]4.12r b) shows the evolution of traffic on downstream 
link di- initially, this link is empty; after t = 0, traffic immediately downstream of the 
diverging point reaches state pc = 0.22pj; along with the propagation of pc, on this 
link, another expansion wave forms and travels downstream; after around 200r, traffic 
density on this link is uniformly pc- Fig. I4.12f c) presents the evolution of traffic on 
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Figure 4.12: Simulation I: A general case 

downstream link initially, traffic density is pd = 0.56pj; after vehicles start to 
diverge at t = 0, traffic immediately downstream of the diverging point reaches pe = 
O.OApj] then a shock forms and travels downstream in a constant speed s = 0.38//r. 
The expansion waves and the shock wave observed on the three branches can be 
shown on the p — q plane as in Fig. I4.12f d). in which dashed line AB represents the 
expansion wave on the upstream link u, dashed line OC represents the expansion wave 
on downstream link di, and solid line DE represents the shock wave on downstream 
link d2, whose slope gives the shock speed. 
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As shown in Figure at the diverging point, there are three traffic states, rep- 
resented by B, C, and E in Fig. I4.12f d). Flow rates at these states are qb = 1.09pj//r, 
qc = 0.87 pjl/r, and qe = 0.22pj//r respectively. We can see that qb = Qc + Qe- Note 
that qb is the outflow of link u, qc inflow to link di, and qE inflow to link d2. Thus 
traffic is conserved at the diverge. Further, qc/qs = 80%, which is the proportion 
of vehicles on the upstream link traveling to di. Therefore, this is consistent with a 
general observation that diverging fl ows are proportionall y determined by the compo- 



sition of traffic on the upstream link (jPapageorgiou , 



3). 



19901 ) . This property of diverging 



flows guarantees that the composition of vehicles on the upstream link never changes 
as observed in our simulation. 



4.5.2 Simulation II: An extreme case 

In this subsection, an extreme case is studied. Initially, like in the previous simulation, 
the upstream branch carries a constant flow with traffic density pu = l.llllpj = 200 
veh/km, the proportion of vehicles on the upstream link traveling to downstream 
branch di is 80%, and downstream link di is empty. However, traffic on downstream 
link d2 is jammed; i.e., = Pj = 180 veh/km. Still, we impose the Neumann 
boundary conditions on the boundary of the diverging network and have the same 
discretization to the three links and the time duration: = 500 and K = 5000. 

Simulation results are shown in Figure I4.1HI In this simulation, traffic density 
on link d2 is uniformly pj as expected and not shown. Fig. I4.13r a) demonstrates 
traffic evolution on the upstream link u: initially traffic density is pa = l-Hpj! 
after the beginning of the diverging process, traffic density immediately upstream 
of the diverging point reaches the jammed density ps = 2pj; then, jammed traffic 
propagates upstream as a back-traveling shock at a speed s = — 0.921 /r. Fig. I4.13f b) 
shows traffic dynamics on link di. In this figure, we can observe small spikes of 
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(a) (b) 




p/Pj 



Figure 4.13: Simulation II: An extreme case 

density, which travel downstream and shrink along time. After around t = 200r, link 
di is almost empty as initially. Fig. 14.13r d) presents traffic states on the diverging 
network on the p — q plane. Here line AB represents the shock wave on the upstream 
link, and its slope gives the shock wave speed. 

From these figures, we can see that, if a downstream link is blocked, the upstream 
link and the other downstream link(s) will be blocked, as expected. However, the 
small spikes in Fig. I4.13f b) suggests there is still a small number of vehicles get 
out of the upstream link at the very beginning. Although this phenomenon seems 
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interesting, the appearing of small spikes in our model is due to numerical error of 
finite difference, caused by the finite number of cells, and are expected to disappear 
if we partition the links into fine enough cells. 

4.6 Discussions 

In this chapter, we have introduced a new approach to modeling diverging traffic 
dynamics, which is represented by instantaneous kinematic waves. We presented 
analytical solutions to these waves and found these solutions are consistent with the 
supply-demand method if the traffic demand is modified as in Section |01 With 
numerical simulations, we assert that this model satisfies two important properties: 
(1) diverging flow to a downstream link is proportional to the fraction of vehicles 
traveling to this link, and (2) the upstream link is blocked after a downstream link is 
blocked. 

The instantaneous kinematic wave theory of diverging traffic is different from 
existing models since it provides more details on diverging traffic dynamics and sheds 
more light at the aggregate level how vehicles segregate themselves at a diverge. 
In the follow-up studies, we would like to discuss the convergent kinematic wave 
solutions when all branches are partitioned into infinite number of cells. We will also 
be interested in finding the location of actual diverging point and how it is related to 
traveler's behavior, the geometry of a diverge, and traffic conditions. Finally we will 
test this model with field data and discuss possible applications. 



Chapter 5 

Kinematic wave traffic flow model 
for mixed traffic 



5.1 Background 

Vehicular traffic on fiigfiways often comprises different types of vefiicles witfi vary- 
ing driving performances. This heterogeneity affects traffic flow characteristics in 
significant ways, a fact that has long been recognized by the transportation engi- 
neering profession. For example, in the computation of flow capacity on a highway 
or at a signalized intersection, the Highway Capacity Manual recommends a series 
of adjustments to take account of the capacity reduction caused by heavy vehicles 
(i.e., trucks/buses/recreational vehicles). If one is interested in the effects of heavy 
vehicles on traffic flow over space and time, however, the Highway Capacity Manual 
procedures are not adequate. For this one needs a dynamic model for mixed traffic. 

Mixed traffic can be modeled at three different levels — microscopic, mesoscopic 
and macroscopic. It is perhaps most straightforward to model mixed traffic on a 
microscopic level — one simply endow, at one extreme, each individual vehicle with 
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different performance and behavior characteristics. Many commercially available sim- 
ulation packages, such as CORSIM, PARAMICS, and VISSIM allow the specification 
of multiple vehicle classes. Major challenges arise when one models mixed traffic on 
a mesoscopic level, mainly due to the correlation between various probability distri- 
butions of vehicular speeds. Nevertheless, a number of me s oscopic models of mixed 



traffic have been developed in recent years (jHelbine 



1997 



Hoogendoorn and Bow . 



20001 ). Aggregation of mesoscopic models of mixed traffic through expectation op- 
erations lead to multi-class traffic flow models in the macroscopic level. There is, 
however, another approach to develop macroscopic mixed traffic flow models. This 
is the approach of continuum modeling. It is this approach that we shall adopt in 
developing a traffic flow model for mixed traffic. 

In the continuum description of traffic flow, vehicular traffic is described as a 
special kind of fluid that are characterized by its concentration (density, p), mean 
velocity (f), and vehicle flux (flow rate q), all are functions of space (x) and time 
(t). The starting point of any continuum model of traffic flow is the conservation of 
vehicles 

d 

— J p{x,t)dx = q{xi,t) - q{x2,t), (5.1) 

and the relation between flow, density, and mean velocity q = kv. 

Equation 15.11 is an integral form of traffic conservation. When the road segment 
[xi, 0:2] shrinks to a point in space, one obtains the familiar differential form of traffic 
conservation: 

Pt + qx = 0, or pt + {pv)^ = 0. (5.2) 

If one introduces a relation between vehicle concentration and traffic speed v = 
V{p), one obtains the classic kinematic wave model developed by Lighthill, Whitham 
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( Lighthill and Whithaml . Il955br i and Richards (jRichardsl . 119561 1 



Pt + {pV{p)), = Q, pV[p)^Q{p). 



(5.3) 



The classic kinematic wave model of Lighthill and Whitham was formulated for 
homogeneous flows on a long crowded road. It does not consider the effects of perfor- 
mance differences among different types of vehicles. Recently, Daganzo extended the 
theory to treat a freeway system with two types of lanes, special la nes and regular 



lanes, and two types of vehicles 



Daganzo et al, 



Daganzo 



prior ity vehicles and regular vehicles (jDaganzo 



1997 



20021 ) . Priority vehicles are allowed to travel on either 



regular lanes or special lanes, whereas regular vehicles can only travel on the regular 
lane. The two types of vehicles in Daganzo's special lane model have different vehicle 
performances in free-flow traffic, but are indistinguishable in heavy traffic, where both 
types of vehicles travel at the same speed. 

In this chapter, we extend the kinematic wave model to vehicular traffic with 
a mixture of vehicle types. In the mixed flow each vehicle type is conserved and 
travels at the group velocity, but the differences among vehicle types are accounted 
for in determining the states of the collective flow. This model can be used to study 
traffic evolution on long crowded highways where low performance vehicles entrap 
high performance ones. It can also give a more accurate description of the I-pipe 
state in Daganzo's special lane model. 

The remaining parts of the chapter are organized as follows. In Section 15.21 we 
give the extended KW model and its basic properties. In Section 15.31 we analyze the 
Riemann problem for this model. In Section l^^ we propose a fundamental diagram of 
mixed traffic and discuss its properties. In Section l531 we provide numerical examples 
and in Section we conclude the chapter. 
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5.2 The extended KW model for mixed traffic 

Let us assume that there are i = l,---,n types of vehicles in the traffic stream 
(n > 2), each type has concentration pi{x,t) and velocity Vi{x,t). By conservation of 
each vehicle type, we have 

{pi)t + (piVi)^ = 0, i = 1,- ■ ■ ,n; (5.4) 

As in the development of the classic KW model, we postulate that equilibrium 
relations exist between vehicular speeds and traffic densities: 

Vi = Vi{pi,---,pn), (5.5) 
with Vi{0) = Vfi, the free-flow speed of each vehicle type and dp.Vi < 0, i = 

1,- j = 1,- 

Equations 15.41 and 15.51 are the general governing equations of mixed traffic flow 
without special lanes. Note that if one adopts Vi = V{Y2^=i Pi), one recovers the I-pipe 
state in Daganzo's special lane model. In this chapter we study a special case of the 
general equations for mixed traffic in which we consider two types of vehicles — one 
represents passenger cars (pi) and the other represents heavy vehicles such as trucks 
(P2), and two traffic flow regimes — free-flow and congested traffic. 

When traffic is light and there are adequate opportunities for passing, different 
classes of vehicles would travel at their own free- flow speeds Vfi. The traffic flow in 
this case can be described by 

{Pi)t + Vfi{pi)x = 0, i = l,2; (5.6) 

By deflning p = Yl^=i Pii ^'^d f / = we can use 

Pt + (pvf). = 0, (5.7) 
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to approximately model the average behavior of light traffic. 

When traffic concentration reaches a critical value pc, passing opportunities di- 
minish and vehicles of lower performance (e.g. trucks) start to entrap vehicles of 
higher performance (e.g., passenger cars). Under such conditions it is assumed that 
the various classes of vehicles are completely mixed and move at the group velocity 
V. That is, mixed traffic flow in this regime is described by 

iPi)t + (PiV), = 0, t = 1,2. (5.8) 

Through the definition of a proper average free-flow speed and the selection of a 
suitable critical density, we combine Equation I5.fil (for free- flow traffic) and Equation 
21 (for congested traffic) into one modeling equation: 



PlVipl,p2) 
p2V{pl,p2) 



0, (5.9) 



where 

y{pi,p2) 



7lPl+72P2<l 



K(pi,P2), 7iPi +72P2 > 1 
Here 71 and 72 are parameters that determine the critical density in [pi, P2) coor- 
dinates. K is a two-dimensional speed-density relation for congested traffic. It is 
understood that dp-V^ < 0, i = 1,2. 

Equation 15. 91 is a system of conservation laws with characteristic velocities: 

\i = V + piVi + P2V2, A2 = V{pi, P2). 

Here we used a special notation for partial derivatives of V with respect to pi and p2 
: dp^V = Vi and dp^V = Vg- Owing to the nature of V{pi, P2), we have Ai < A2 = V, 
that is, both characteristics travel no faster than average traffic. In fact, the second 
characteristic travels at precisely the speed of traffic. When the free-flow speeds of 
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both types of vehicles are identical, the extended KW model preserves the anisotropic 
property of the KW model. Otherwise, the extended model is not anisotrop ic in lig ht 
traffic (the nature of this violation of anisotropy is explained in detail in (|Zhanei )). 
Moreover, it can be shown that 



P2 



+ V 



PlJt VPl/x 



P2 



0, 



from which one obtains ^ = 0, that is, the level curves of in the t — x 
plane coincide with vehicle trajectories. The separation of level curves therefore 
implies first-in-first-out traffic fiow behavior between vehicle classes. 

Furthermore, the corresponding eigenvectors of the fiow Jacobian matrix are 











ri = 1 




1 — 1 






\ 1 ) 







and the Riemann invariants (w, 2), defined as \jw •ri =0, \jz • r2 = 0, are 



w 



Pi 



They are used here to obtain the expansion wave solutions of a Riemann problem 



(see next section. 



refer to ( Whitham 



or m ore details on Riemann problems and Riemann invariants. 



\mM)- 

It can be shown that the first characteristic field is nonlinear and the second 
characteristic field is linearly degenerate. We therefore have both shock and smooth 
expansion waves in the first field and contact waves (or slips) in the second field. 
We shall derive the expressions for these waves related to Riemann data in the next 
section. 
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5.3 The Riemann problem and basic wave solu- 
tions 

In this section we discuss the solutions of the extended KW model, Equation 15.91 
with the following so-called Riemann data: 



p(x,0) 



a; < 



P 




(5.10) 



To solve the above Riemann problem, we first study the right (downstream) states 
that can be connected to the left (upstream) states by an elementary wave, i.e., a 
smooth expansi on (r arefaction) wave , a con t act, or a shock (readers are referred to 



Zhand (|200n 



2002 



2001br i for a more 



(jLeVeauel . 120021 ) and lAw and Rasclel ()2000D : 

detailed discussion of Riemann problems related to systems of conservation laws in 
general and traffic flow in particular). Throughout the remaining sections, we assume 
that Vfi = Vf2 = Vf. This assumption ensures that our proposed V{pi,p2) function is 
continuous over the entire feasible (pi, P2) region. The Riemann problem of Equation 
15.91 with discontinuous V{pi,p2) is more involved and will be discussed elsewhere. 
Nevertheless, the analysis of the model. Equation 15. 9[ with Vfi = Vf2 = Vf still 
reveals many key features of mixed traffic flow. 

The 1-expansion waves: An upstream state p' can be connected to a downstream 
state p*" by a 1-expansion wave if and only if the downstream state satisfles 

w(p') = w{p''), p' > p^ 



I.e., 



P2 P2 



P\ 



P\ 



(5.11) 



This means that in the p— plane the two states are on a ray from the origin. Clearly 
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across an expansion wave traffic composition does not change, that is, vehicles observe 
the first-in-first-out rule. 

The contact waves: 

A contact wave is a slip that separates two traffic regions of different traffic 
densities and vehicle compositions but the same travel speed. That is, 

V{p') = V{p^). (5.12) 

In the p— plane, all states on a level curve of V{p) are connected by a contact wave. 

The shock waves: 

The shock waves in the extended KW model are given by the jump condition: 

s{p{-p\) = p[V{p[,p',)-p\V{p\,pl) (5.13) 
S{p'2-Pl) = P2V{p{.p'2)-plV{p\,pl) (5.14) 

After elimination of s from the equations and some algebraic manipulations one ob- 
tains 

{P\P2 - P\P2){V{p[,p',) - V{p\,pl)) = 0. 

Two possibilities exist: V{p\,p2) — V{p\,p2) — which gives the contact waves that 
we have already discussed, or 

PiP2-PiP2 = 0, (5.15) 

this gives the downstream states p^ that can be connected to the upstream state 
by a shock. Note that all these states also fall on a ray originating from the origin 
of the p— plane. This implies that across a shock vehicle composition also does not 
change, that is, vehicles observe first-in-first-out rule. Moreover, we have the following 
entropy conditions 

p' < 
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to ensure the stability of the shock. 

Now we can state the procedure to solve a Riemann problem for the extended KW 
model. Note that for any state p\ the two curves/lines given by Equations 15. Ill and 
15. 121 divide the feasible p— plane into four regions (FigureE^- If the downstream state 
falls on any of these two curves/lines, it can be connected to the upstream state by 
an elementary wave. If it falls on any of the four regions, however, an intermediate 
state is generated on the line given by Equation 15. IH which is connected with 
the upstream state by a 1-wave (i.e., an expansion or shock wave) and with the 
downstream state by a contact (Figure 15.11 Figure 15.21 shows a few examples of 
Riemann solutions. 



5.4 Fundamental diagrams for mixed traffic 



We propose the following p — V rela tion, which can 



3e der ived from Following 



model under steady-state conditions (|Zhang and Kim 
traffic flow model. 



20001 1 , to be used in the mixed 



V 



"'ZZT^ ili + r,Vf{)p^ + {h + r,Vf,)p,<l 

P1T1+P2T2 



{h + riVfi)pi + {k + T2Vf2)p2 > 1 

This fundamental diagram is shown in Figure EISl 

We call the above relation the extended speed-density relation for the triangular 
fundamental diagram. The parameters are: free flow speeds for both types of vehicles 
Vfi and Vf2, effective vehicle lengths for type-1 and type-2 vehicles h and I2, and 
response times of type-1 and type-2 vehicles ri and T2. The last two parameters 
capture, to a certain degree, the acceleration/deceleration differences between the 
two classes of vehicles. 

The capacity of mixed flow depends on vehicle composition. For example, in the 
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2jam^ 



1-shock curve/line , p2/p^=constant 




contact curve/line, V=constant 



1 -expansion curve/line, 
p2/p2=constant 



Figure 5.1: Phase diagram for determining elementary and simple waves 



case of Vfi = Vf2 = Vf, let ^ = p < oo, then the critical densities of the proposed 



speed-density relation are 

Pic = 



1 



P2c = PPlc 



ih +Pk) +P'T2) 
Note that when p = 0, i.e., there are no type-2 vehicles in the traffic stream, we 
recover the critical density for type-1 vehicles 

1 



Pic 



h + VfTi 

and when p = oo, i.e., no type-1 vehicles present in the traffic stream, we can switch 
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the positions of p in relation to Vs and r's in the above formulas and replace it with 
0. Again we recover the critical density for type-2 vehicles 



p 



P2c 



I2 + VfT2 



For any vehicle mixture (i.e. < p < 00), the capacity flow is ^ 



l2+VfT2 ' h+VfTl J ' 



5.5 Numerical solution method and simulations 



5.5.1 The Godunov method 



We a pproximate the mixed traffic KW model with the Godunov method (|Godunov 



Pl,i ^ Pl,i Pl,i+l/2^ yPl,i+l/2^ P2,i+l/2) ~ Pl,i-l/2^ yPl,i-l/2^ P2,i-l/2) _ 

At ^ a;^^ ~ ' 

P2,i ~ P2.i P2,i+l/2^ \Pl,i+l/2^ P2,i+l/2) ~ P2,i-l/2^ \Pl,i-l/2^ P2,i-l/2l _ „ 
At ^ A^ ~ ' 

in which p-^ ^ is the average of pk over cell i at time tf i.e., p;^ ^ = J^^^^J'/^^^^ pk{x, tj) dx/ Ax, 
and pI\_i/2 is the average over time interval [tj,tj+i] at the boundary Xi-1/2 between 
cells i and i — i.e., p*ki_i/2 = /t=t^ Pkixi-i/2,t) dt/At. Given {pi, P2) at tj, we hence 
can compute traffic states at the following time tj+i. 

In the above equations, the boundary average P*ki-\i2 '^^'^ t)e found by solving 
the Riemann problem for the extended KW model. Equation 15.91 with the following 
initial condition (p' = (pi,,_i, p^2,*-i) and p" = {p^i, pi^,)) 

p\ if X - Xj_i/2 < 0, 
p\ if X - Xj_i/2 > 0. 

As shown in section 2, the Riemann solutions consist of a shock or rarefaction wave 
with an intermediate state p™ and a contact wave. Since all the waves are self-similar. 
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Shock 


s riven in Equation 15. 181 


Pl,i-l/2 




s > 


p[ 


P'. 


s < 


pT 


P2 



Table 5.1: Shock wave solutions in mixed traffic 

we have p*-'_-^^^=p{xi^i/2,t)= const for all t > 0, which is determined by the shock 
or rarefaction wave connecting and p™ since the contact wave has non-negative 
velocity and is not involved. 

From Equations 15.111 and I5.15| we have 

P2 P2 



pT' 



(5.16) 



and from Equation 15.121 



(5.17) 



Combining Equations 15.161 and I5.17( we can find the intermediate state p"*, from 
which we can compute P*''_ii2 described in the following cases. 

Case 1 When < p^, they are connected by a shock, and the shock speed is 

^ ^ p\v{p{,p^,)-p'^v{pr.p^) 

p{-pT 

In this case, solutions of P*''_ii2 summarized in the Table IHTTl 



(5.18) 



Case 2 When p' > p"^, they are connected by a rarefaction wave, in which the charac- 
teristic velocity is Ai(pi,p2), and Ai(p') < Ai(p™'). If Ai(p') > 0, P*''_ii2 is the 
same as the left state p'; if Ai(p'") < 0, it is the same as the intermediate state 



p™. Otherwise, Pj_]^/2 satisfies 



Al(Pi^j_l/2,p£_i/2) 
^2,1-1/2/ Pl,i-l/2 



p\Ip\ 



(5.19) 
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Rarefaction 


Ai 


Pl,i-l/2 


^2,1-1/2 


Ai(p') >0 




pi 


Ai(p™) < 


pT 


P2 


o.w. 


riven in Equation 15.191 



Table 5.2: Rarefaction wave solutions in mixed traffic 

which implies that p*''_i/2 maximizes the total flow (pi + p2)V{pi, P2) along the 
line P2/P1 = p\Ip\. 

In this case, therefore, solutions of p*^_i/2 summarized in Table 
5.5.2 Numerical simulations 

In our simulations, we will use the extended triangular fundamental diagram (Figure 
15. 3|) . in which the parameter values are: free flow speed for both types of vehicles 
Vfi = Vf2 = Vf = 65 mph = 95.3333 ft/sec, effective vehicle lengths for type-1 and 
type-2 vehicles: h = 20 ft, /2 = 40 ft, and response times of type-1 and type-2 vehicles: 
Ti = 1.5 s, r2 = 3 s. Therefore, we have pijam = l/^i = 0-05 veh/ft and p2,jam = 
I/I2 = 0.025 veh/ft. Moreover, since h/h = ^1/^2, we have that in the extended 
triangular fundamental diagram V{pi, P2) is a function of pi/pijam + P2/ P2,jam- Thus, 
as we will see later, the evolution pattern of travel speed is the same as that of 

Pll Pl,jam + P2/ P2 Jam,' 

We will conduct numerical simulations for a ring road, whose length L = 2000/i=40,000 
ft, during a time interval from t = to T = 100ri=150 s. In order to apply 
the Godunov method, we partition the ring road into = 1000 cells with length 
Ax = L/N=AO ft, and discretize the time interval to M = A^/2 steps with the 
duration of each time step At = T/M=0.3 s. Since the CFL condition number 
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(jCourant et all Il928h 



max{|Ai|, |A2|}|^ < = 0.715 < 1, 

the Godunov method yields convergent solutions. 

For the numerical simulations, we use the following global perturbation as initial 
traffic conditions, 

pi(x,0) = (0.2 + 0.16 sin(27rx/L))pi,,-,^, 

P2(a;, 0) = (0.15 + 0.1 sin(27rx/L))p2,jam, 
in which the density of the 1-type vehicles, i.e., has higher average but smaller oscil- 
lation. 

With these conditions, the solutions of the mixed traffic flow model are depicted 
as contour plots and shown in Figure The horizontal axis in each of the sub flgures 
represents space and the vertical axis time. These flgures depict the traffic conditions 
(speed and density). As shown in the contour plots of v and pi/ pijam + P2I P2,jami 
there are roughly two traffic regions along the ring road initially. In one region 
(/i + TiVf)pi + {I2 + T2f/)p2 > 1, waves initiated from this region travel backward in 
the same speed, which is Ai(pi,p2) = -{pih + P2^2)/(pit-i + P2T"2) = -Ii/ti = -40/3 
ft/sec. In another region, (Zi + TiVf)pi + {I2 + T2f/)p2 < 1, waves initiated from this 
region travel forward at free-flow speed. Two waves separate the two regions: an 
expansion wave on the left and a shock wave on the right. The shock wave travels 
forward initially but eventually travel at — /i/ri as traffic density increases to critical 
density in the second (free-flow) region. The patterns of the solutions of pi and 
P2, however, are not the same as that of v because the change of the overall traffic 
conditions affect each vehicle class differently. 

The contour plot of P2/P1 is given in Figure from which we can see that the 
level curves do not intersect. Remembering that the level curves of P2/ pi coincide with 
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vehicle trajectories, the solutions shown here indicate that under the given conditions 
first-in-first-out rule is respected by the mixed flow model. From this figure we can 
also see the expansion and shock waves as they move through traffic, which are marked 
by changes in the slopes of the level curves. 



5.6 Concluding remarks 

In this chapter we extend the Lighthill-Whitham- Richards kinematic wave traffic flow 
model to describe traffic with different types of vehicles, where all types of vehicles 
are completely mixed and travel at the same group velocity. A study of such a model 
with two vehicle classes (e.g., passenger cars and trucks) shows that,when both classes 
of traffic have identical free-flow speeds, the model 1) satisfies first-in-first-out rule, 
2) is anisotropic, and 3) has the usual shock and expansion waves, and a family of 
contact waves. Different compositions of vehicle classes in this model propagate along 
contact waves. Such models can be used to study traffic evolution on long crowded 
highways where low performance vehicles entrap high performance ones. 
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Figure 5.2: Wave solutions to the Riemann problem: Shock + Contact wave (left) 
and Expansion wave + Contact wave (right) (In the bottom figures, thick (dashed) 
lines are characteristics, and lines with arrows are vehicles' trajectories.) 
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Figure 5.3: The extended triangular fundamental diagram 



CHAPTER 6. 



MULTI-COMMODITY NETWORK TRAFFIC MODEL 130 




Figure 5.4: Contour plots of solutions on the x — t space with the extended triangular 
fundamental diagram 
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Figure 5.5: Contour plot of P2/P1 on the x — t space with the extended triangular 
fundamental diagram 



Chapter 6 

Kinematic wave simulation model 
for multi-commodity network 
traffic flow 

6.1 Introduction 

Recurrent or non-recurrent traffic congestion in many major metropolitan areas liave 
seriously deteriorated the mobility, convenience, and productivity of highway trans- 
portation. To tackle the congestion problem, traffic engineers and scientists are facing 
a number of challenges, including: 1) the evaluation of the performance of a road 
network (e.g. total travel time or the level of service (LOS)), 2) the prediction of 
occurrence of incidents, 3) the development of traffic control schemes (e.g. ramp me- 
tering) and management strategies (e.g. traffic guidance) in Advanced Transportation 
Management and Information Systems (ATMIS) or other Intelligent Transportation 
Systems (ITS) , and 4) the estimation of traffic demand associated with a given ori- 
gin/destination (i.e. 0/D information). As we all know, a fundamental and essential 
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issue in addressing all these challenges is the understanding of traffic dynamics on a 
road network, i.e., the evolution of traffic on a road network under initial and bound- 
ary conditions, for which traffic flow models of road networks play an important role. 

Among many traffic flow models, the kinematic wave models are advantageous 
in studying traffic dynamics in a large-scale road network, due to their inherit com- 
pliance with ITS applications, theoretical rigor, and computational efficiency. Practi- 
cally, people are more interested in aggregate-level traffic conditions such as average 
travel speeds, densities, flow-rates, and travel times, which are directly concerned or 
can be easily derived in the kinematic wave models. Theoretically, the evolution of 
traffic conditions on a link can be studied as hyperbolic conservation laws, given the 
fundamental diagram, which deflnes a functional relationship between density and 
flow-rate or travel speed, and traffic conservation, which means that the change of 
the number of vehicles in a section in a time interval equals to the number of vehi- 
cles entering the section minus the number of vehicles leaving it. Computationally, 
the kinematic wave models can be solved with the Godunov method, in which in- 



flux and out-flux of a eel 



supply-demand method (|Daganzo 



can be computed through 



dnematic wave solutions or the 



1995a; 



Lebacqud . 



19961) . 



In literature, the following approaches can be considered as typical when modeling 



Vai]gha.n et al 



netwo rk traffic flow in the framework of kinematic wave theory. First 
(Il984[ ) studied network traffic flow with two continuous equations: a "local equation" , 
which ensures traffic conservation and is consistent with the traditional LWR model on 
a link at the aggregate level, and a "history equation", which computes the trajectory 
of each vehicle at the disaggregate level. Since the trajectories of all vehicles are not 
always required for many applicati ons, this model tend s to waste signiflcant amount of 



computational resources. Second, 



JavakrishnanI ()l991[ ) introduced a discrete network 



flow model, in which each link is partitioned into a number of cells, vehicles adjacent 
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to each other and with the same 0/D and path are called a "macroparticle" , and 
the position of a macroparticle at each time step is determined by its travel speed 
and the cell's length. However, with the given mechanism, this model may not be 
consistent with the LWR model since traffic conservation can be violated. Third, 



Daeanzo 


( 


1995a) 


1 

Model (IDaeanzo 



1994^ . which is a numerical approximation of the LWR model with 



a special funda mental diagram — a triangle. In this discrete model, macroparticles 



in the sense of 



JavakrishnanI (jl99l|) in a cell are ordered according to their waiting 



times and moved to the downstream cell when their waiting times are greater than a 
threshold minimum waiting time, which is computed from traffic flow at the aggregate 
level. However, the determination of the threshold mi nimum waiting tim e is quite 



tedious. Fourth, in the KWave98 simulation pack age (ILenonard II . 



based on the simplifled kinematic wave theory by iNewel] 



19981 ). which is 



19931 1. vehicles on a link 



are considered to be ordered as a queue. Then, with the in-queue and out-queue of 
the link determined, one can easily update the link queue. However, in-queues and 
out-queues at typical highway junctions such as merges and diverges create difficulties 
for this model. 

In many applications, such as dynamic traffic assignment (DTA), First-In-First- 
Out (FIFO) principle is a key concern. Consequently, all the aforementioned models 
order veh icles but in diffe r ent fa shions. Vehicl es are ordered a ccord ing to their trajec- 



tories in 



flDaganzn 



Vaughan et al. 



198^, locations in (iJavakrishnanl. 



1991 



1995a[ ). and positions in a link queue in (jLenonard II , 



), waiting times in 



19981 ). That is, these 



simulation models track either vehicles' trajectories, or positions, or waiting times, 
or que ue orders. ^ In t hese models, therefore, traffic conditions are also considered 



Javakrishnan . 



19911) , although the introduction of macroparticles save a certain amount of 



memory in computation, this saving is diminished when vehicles of different classes are evenly mixed. 
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at the vehicle level. As a result, the computational efficiency of the kinematic wave 
models is not fully utilized. 

In order to fully explore the computational efficiency inherited in the kinematic 
wave theory, we propose a new network traffic flow simulation model, in which traffic 
dynamics are studied at the aggregate level with commodity specified densities and 
flow-rates. Traffic of a commodity has the same characteristics, which can be vehicle 
type, destination, path, or any classification criterion. We hereafter refer to this model 
as multi-commodity kinematic wave (MCKW) simulation model. In this dissertation 
research, we start with a simple traffic system, where vehicles are categorized into 
multiple commodities according to their paths and no differentiations are made in 
vehicle types, driver classes, or lane types, such as high-occupancy-vehicles (HOV's) 
or HOV lanes. In the MCKW simulation platform, besides link characteristics such as 
free flow speed, capacity, and number of lanes, traffic dynamics are highly related to 
geometrical characteristics of a road network including link inhomogeneities, merges, 
diverges, and other junctions. Although it does not give a complete, detailed picture 
of traffic dynamics, this model is still of importance for many applications, in which 
the difference in vehicles, drivers, or lanes are negligible not interested or can be 
averaged out without major loss of accuracy. 

In the MCKW simulation platform, the Godunov-type approximation is applied, 
and fluxes through boundaries inside a link and junctions such as merges and diverges 
are computed systematically in the framework of the supply-demand method. In addi- 
tion, since vehicles are categorized by vehicles' origins, destinations, and paths, traffic 
dynamics are also differentiated for different commodities in the MCKW simulation 
platform. We will show that FIFO is observed in the commodity specific kinematic 
wave theory. 

In the rest of this chapter, we will discuss theories underlying the MCKW model 
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in Section and the data, network, and program structures in the MCKW model 
in Section 16 .31 In Section 16 .41 we will discuss the process of output from the MCKW 
model to obtain interested data of a road network, such as travel times. In Section 
16.51 we carry out some numerical simulations of a simple road network. In Section IHSl 
we will draw some conclusions and provide further discussions about the calibration 
of the MCKW model. 



6.2 Underlying theories of the MCKW simulation 
model 

In the multi-commodity kinematic wave (MCKW) simulation of network traffic flow, 
dynamics of total traffic, i.e., the evolution of traffic conditions of all commodities, are 
studied at the aggregate level and governed by the kinematic wave theories. These 
theories have been studied in the previous chapters for fundamental network com- 
ponents such as inhomogeneous links, merges, and diverges. They provide the un- 
derlying algorithms and building blocks for the MCKW simulation platform. At the 
disaggregate level, traffic of each commodity in the form of proportions is studied 
with its proportion, and First-In-First-Out property on a link will be discussed. 



6.2.1 Kinematic wave theories at the aggregate level 



In the MCKW model, we use the discrete form of the kinematic wave t 



can be obtained through the first-order^ Godunov method (|Godunov 



leories, which 



19591 ) of the 



continuous versions. In the discrete form, each link is partitioned into cells, of 
equal length or not, and the time inter val is discretize d into K time steps. Then, we 



second-order method was discussed in (jDaganzo . 



1999a 1 
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obtain the Godunov-type finite difference equation for total flow in cell i from time 
step j to time step j + 1 as follows: 



At 



Pi ^ h-i/2 



f3* 

Ji+1/2 



Ax 



0, 



(6.1) 



rom time step ?" to t ime step j + 1 



1928^ condition. In 



where Ax is the length of cell i, At is the time 
and the choice of is governed by the CFL ()Courant et al. 
Equation 16. H is the average of traffic density p in cell i at time step j, similarly 
p]'*'^ is the average of p at time step j + 1; //^i/2 ^^-^ through the upstream 

boundary of cell i from time step j to time step j + 1, and similarly 
downstream boundary flux. Given traffic conditions at time step j, we can calculate 
the traffic density in cell i at time step j ' + 1 as 



P. 



Ax 



ri ^ Wi-l/2 



1/2)- 



(6.2) 



Defining = pi Ax as the number of vehicles in cell i at time step j, N^^^ = 
pi~^^Ax as the number of vehicles at time step j + 1, ~ fiPi-1/2)^^ 

number of vehicles flowing into cell i from time step j to j + 1, and F.^^^^ as the 
number of vehicles flowing out of cell i, Equation 16.21 can be written as: 



1/2' 



(6.3) 



which is in the form of traffic conservation. 

Given the init i al and boundary conditions, we will use the supply-demand method 



( Daganzo 



Lebacquf . 



19961 ) for computing fluxes through cell boundaries: 



^i-1/2 fi-i/2- ■'■^ ^ general road network, there are the following types of boundaries: 



1/: 



boundaries inside a link, merges, diverges, and more complicated intersections. 

1. When the boundary at Xj_i/2 is a boundary inside a link, whose upstream cell 
is denoted as u and downstream cell d, we follow the supply-demand method 
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(6.4) 



discussed in (|Daganzd . Il995al : iLebacqud . Il996l : Ijin and Zhand. l2003br i. I.e., if 

we define the upstream demand as 

Q{Uu), when Uu is under-critical 
Qmaa-^ when Uu is over-critical 

and define the downstream supply as 

Q'J^"'^ , when Ud is under-critical 
Q{Ud), when Ud is over-critical 
then the boundary flux can be simply computed as 



Sd 



(6.5) 



fj* 

Ji-l 



/2 



minlZ^^, Sd}, 



(6.6) 



where Ud and U^ are traffic conditions including density p and road inhomogene- 
ity a at jth t ime step, of the downs tream and upstream cells, respectively. As 



discussed in (j.Tin and Zhang 



2nn3b( ). this method is consistent with analytical 



solutions of the Riemann problem for inhomogeneous roadway. 

2. When Xi-1/2 is a merging junction with P upstream merging cells, which are 
denoted a.s Up {p = 1, ■ ■ ■ , P), and a downstream cell d. The demand of upstream 
cell Up, Dp, is defined in Equation 16.41 and the supply of the downstream cell. 



Sd, is defined in 



fjjin and Zhang . 



]quat 



ion 16.51 Then, we apply the simplest distribution scheme 



2003d ) and compute the boundary fluxes as 



Ji-l/2,( 

fj* 
i- 



'^^HY.p=iDp,Sd}, 



p 



(6.7) 



' j— 1/2, Up 

where f^'*, , is the in-fiow of the downstream cell, and f^'*, ,„ is the out-fiow 

•'2— i/z,a ' •'2— 1/2, tip 

of upstream cell Up. 

In addition, if an upstream cell, e.g. Up, is signalized and denote r as the 
proportion of green light in a cycle (i.e. green ratio), then we can apply the 
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control 



ed traffic demand of Up, minlrO"^"^, -P p^' the supply-demand method 



above (jDaganzg 



1995a; 



Jin and Zhang 



2003d l. Note that r can be a continuous 



function when considering the average effect or a piece-wise constant function 
when the simulation interval At is smaller than a signal cycle. 

3. When Xi-1/2 is a diverging junction with P downstream cells, which are de- 
noted as dp iv = 1, • • ■ , -P), and an upstream cell u. In the model proposed in 



fjjin and Zhane . 



20021 ) ■ we introduced a new definition of partial traffic demand 



of vehicles traveling to Dp in cell u as follows, 



(6.8) 



Q{pp] pp) pp is UC 
Q""^^{pp) otherwise 

where pp is equal to the density of vehicles not traveling to dp, and at critical 
density Q{p] Pp) reaches its maximum. The traffic supply for dp, Sp, is defined 
by Equation 16.51 Then, the boundary flux to dp, f-li/2dp^ '^^^ computed by 



fj* 

Ji-l 



/2,dp 



Tsim.{Sp,Dp}, 



(6.9) 



and the out-flow of u, f-_i/2 u-> ^^'^ these flows, 

p 

J i-l/2,u J i-l/2,dp ■ 

p=l 



(6.10) 



Another model of traffic diverging to D downstream links we will implemen t in 



the MCKW simulation was proposed in (jDaganzo 



1995al : 



Lebacquf 



im): 



fi-i/2,u = ^i^d=i{Du,Sd/^d}, 

d=l,---,D, 



fj* 

J i 



i/2,u ^, ^, 

— c .fj* ^ _ 1 . . . n 

'j-l/2,(i 

where is the proportion of commodity d in total traffic, and here is the 
demand of the upstream cell as defined in Equation 16.41 
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When vehicles have no predefined route choice an d can choose everv down stream 



link at a diverge, we use the model proposed in (jjin and Zhang . 



2003d): 



Ji-l/2,d c;Ji_i/2,u^ 



d 



D. 



(6.12) 



4. For intersections with two or more upstream and downstream links, we can 
combine the merge and diverge models together. Note that only the computa- 
tion of demands and supplies may change, and the supply-demand method is 
still the same. 

For example, when we combine the supply-demand methods in Equation 16 . 71 and 
Equation l6.11l for an intersection with U upstream branches and D downstream 
branches, we can compute fiuxes by 

mmt,{E:^=i Sal i ^h'^fA }, 



fj* 

Ji-l/2 



Ji-l/24 ~ D •'j-1/2' a— 



1/2,0 
1/2,1 



(6.13) 



fj* 

Ji-l/2,u 



Du 



J 'I 



J* 



U 



l,---,^7, 



where ^u,d is the proportion of traffic heading downstream link d in upstream link 
u, f-*i/2 total fiux through the boundary, f-*^/2 d heading downstream 

link d, and fi*i/2 u from upstream link u. In this model, the intersection 
is considered as a combination of a merge with U u pstream branche s and a 



diverge with D downstream links in the fashion of (jDaganzo 



1995ah . Note 



that the merge model. Equation 16. 7[ and the diverge model. Equation 16. Ill are 
specific cases of Equation 16.131 



6.2.2 Commodity-based kinematic wave theories 

In the MCKW simulation platform, commodities are differentiated by origin/destination 
or path. We assume that a road network has P' origin/destination (OD) pairs and P 
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paths (P > P'). When vehicles have predefined paths, we then have a P-commodity 
traffic fiow on the road network and label vehicles taking pth path as p-commodity. 
When vehicles of an 0/D have no predefined paths, we have P'-commodity traffic 
ffow. 

In the kinematic wave theories of multi-commodity traffic, we denote total traffic 
density, travel speed, and fiow-rate respectively by p, f , and g, which are all functions 
of location x and time t. In contrast, these quantities for p-commodity vehicles 
are Pp, fp, and qp respectively. The fundamental diagram of total traffic defines a 
functional relationship between density and travel speed or fiow-rate: q = Q{a,p) 
and V = V{a,p) = Q{a,p)/p, where a{x) stands for road inhomogeneities at location 
X such as changes in the number of lanes, curvature, and fre e fiow speed. Further, we 



assume traffic on all links is additive in the following sense (j.Tin and Zhang . 

p 

P = J2pp^ 



V 



p=i 

Vp = 

p 

p=i 



V{a,p), 



p=l, 



(6.14) 
(6.15) 
(6.16) 



The kinematic wave theory of additive multi-commo dity traffic on a link can be 



described by the following theory (|,Tin and Zhangl . 1200811 ) , 

pt + Qia,p)^ = 0, 

ipp)t+ ippV{a,p))^ = 0, p=l,-- 
If denoting the local proportion of p-cor nmodity (p 



P. 



(6.17) 



have the following advection equations (jLebacque 



P) by ^p = Ppl p, we then 



{Qt + V{a,p){Q, 



0, 



P 



P. 



(6.18) 



From Equation 16.181 we can see that proportions of all commodities travel forward 
in a link along with vehicles in traffic fiow, as the change of C,p in material space. 
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{^p)t + V{a, p){^p)x, equals to zero. This is also true for all kinds of junctions, in par- 
ticular diverges, in their supply- demand models in the preceding subsection ^ There- 
fore, Equation 16. 181 also means that the profile of proportions coincides with vehicles' 
trajectories on a link. That is, if two or more commodities initially completely are 
divided by an interface, this interface will move forward along with vehicles on both 
sides of, and these commodities will never mix. Since each single vehicle can consid- 
ered as a commodity, all vehicles' trajectories keep disjoint in the commodity-based 
kinematic wave models. Therefore, FIFO principle is respected in this continuous 
model. 

In the previous subsection, we studied the discrete kinematic wave theory for total 
traffic. Here, we will present the discrete kinematic wave theory for each commodity. 
Given traffic conditions of p-commodity at time step j, i.e., j in all cells, we can 
calculate the traffic density of p-commodity in cell i at time step j ' + 1 as 

At 



i+i 

Pp,i 



(6.19) 



where /p*_i/2 in-flux of p-commodity through the upstream boundary of cell i 

during time steps j and j + 1, and /p*_,_i/2 out-flux. Furthermore, since the profile of 
the proportion of a commodity always travels forward at traffic speed, the proportion 
of a commodity in out-flux of cell i compar ed to all commo dities is equal to the 



proportion of the commodity in the cell. I.e. (|Lebacque 



19961) 



■fp,i+l/2 ' Pi fi+1/2 



p. 



(6.20) 



This is true for cells right upstream of rnerging iunctions fiJin and Zhane , 



diverging junctions (|Papageorgiou . 



2001a|). 



199G; 



Daganzo 



1995a; 



1 

flJin and Zhane 


2003c 


) and 


LebacquJ. 


1996; 


Jin and Zhane 



^That traffic is anisotropic is believed to regulate this property. 
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During time steps j and j + 1 at a boundary Xi+1/2, which has U upstream cells 
and D downstream cells, if we know the out-flux from upstream cell u {u = 1, ■ ■ ■ ,U), 
fp*u i+1/2 (P ~ ■ ■ ■ 5 -P)) "we can obtain the in-fiux of downstream cell d{d = 1, ■ ■ ■ , Z)), 

fp*di+i/2 (P = ■ ■ ■ 5 -O)' from traffic conservation in p-commodity: 

u D 

fp,u,i+l/2 = X/ fp,d,i+l/2- (6.21) 
u=l d=l 

However, when p-commodity vehicles can take more than one downstream cells, we 
have 

p 

X^/p,d,j+i/2 ^ fd,i+i/2- (6.22) 
p=i 

Note that, in Equation l6.171 the kinematic wave solutions are determined by those 
of total traffic, which are obtained by the first-order convergent Godunov method. 
Also from Equation I6.18[ we can see that Equation 16.201 yields an up- wind method 
for in Equation 16.191 Therefore, the discrete model for the commodity-based 
kinematic wave model. Equation 16. 171 converges in first order to continuous version, 
whose solutions observe FIFO principle. That is, in numerical solutions, error in 
travel time of any vehicle is in the order of At. That is, in the MCKW simulation, 
FIFO is accurate to the order of At and Ax. Therefore, when we decrease At, this 
approximation becomes more accurate. 

6.3 Network structure, data structure, and pro- 
gram flow-charts in the MCKW simulation plat- 
form 

In this section, we will discuss the programming details of the MCKW simulation of 
an illustrative road network. 
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Figure 6.1: A demonstration road network 

6.3.1 Network structure 

For the purposes of exposition, a simple traffic network, shown in Figure Ifj.lf is 
considered, and these discussions can be extended to more general road networks. 

The road network in Figure 16.11 where the arrows show traffic direction, consists 
of one origin/destination (0/D) pair and four links. In this network, there are two 
paths. Furthermore, we assume vehicles have predefined paths. ^ Thus, traffic flow 
on this road network consists of two commodities. 

In the MCKW simulation, origins and destinations have the same data structure 
as regular links and are treated as links. For a road network with nuin_origin origins, 
num_od origins and destinations, and num_link links, all links are numbered: origins 
from ^ to nuni_origin-l, the number of origins; destinations from num_origin to 
num_od-l; and regular links from nuni_od to num_links-l. In the sample network, links 
are numbered as shown in Figure IU?T1 

With these numbers, the paths or commodities are denoted as follows: commodity 

^When vehicles have no predefined paths, these discussions are also applicable. 
^The numbering starts from rather than 1 according to C/C++ conventions. 
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takes links 0, 2, 3, 5, and 1, and commodity 1 takes links 0, 2, 4, 5, and 1. This 
is equivalent to saying that links 0, 2, 5, and 1 carry 0-commodity and 1-commodity 
flows, link 3 carries only 0-commodity flow, and link 4 carries only 1-commodity 
flow. Further, the network structure and traflic flow direction is represented by the 
upstream and downstream links of each link. For example, the upstream links of link 
5 and the downstream links of link 2 are links 3 and 4, respectively. 

In the MCKW simulation, each link is partitioned into a number of cells. ^ Since 
fluxes through cell boundaries are important in computation in the kinematic wave 
theories, cell boundaries are also included in the structure of a link. ^ In the MCKW 
simulation platform, cells and boundaries are ordered according to traffic direction: 
adjacent cells and boundaries are either upstream or downstream to a cell. 

In the MCKW simulation, therefore, network can be constructed if we know the 
structures of all links, commodities on a link, and the upstream and downstream 
links of all links. That is, junctions are not used to store network structure although 
they are also numbered in Figure IHTTl The representation of network structure in the 
MCKW simulation largely simplifies the data structure, in which only links are used. 

6.3.2 Data structure 

In the MCKW simulation platform, the structure and characteristics of a road network 
and traffic conditions are all represented by links as well as cells and boundaries inside 
a link. Therefore, the major data structure is linkType, through which network and 
traffic conditions are dealt with, together with sub-structures for cells, cellType, 
and boundaries, boundaryType. The data structures are shown in Figure 16.21 and 

explained in detail as below. 

^Origins and destinations have only one cell. 

''Note that there are one more cell boundaries than cells. 
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arrivalCurves [ ' 
departureCurves [ ] 



numUpLinks 
upLinks [ ] 
numDownLinks 



downlinks 



typeUpJunction 
typeDownJunction 
downBoundary ■ 
downBoundary ■ 



boundaryType 



I 



oundaryTyp 



T 



cellType 



It 



flux 
proportion [ ] 



1 



H [ 



boundaryType 



Figure 6.2: Data structure in the MCKW mckw platform 
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In Figure ins the data structure for a link, linkType, is shown in the left box. 
Its fields are explained as follows: 

downCell is an address pointing to the furthest downstream cell of a link. 
upCell is an address pointing to the furthest upstream cell of a link. 
numCell is the number of cells in a link. 

numCommodity is the number of commodities traveling on the link. 

commodity [ ] is an array of commodities, whose length is numCommodity. Com- 
modities are ordered increasingly for all links. 

arrivalCurves [ ] is an array of accumulative flow entering the link. Its length is 
numCommodity+1. The first numCommodity entries store the cumulative flows of 
corresponding commodities, and the last entry stores total cumulative flow. 

departureCurves [ ] is the same as arrivalCurves [ ] except that we consider 
exiting flows. 

numUpLinks is the number of upstream links adjacent to a link. 

upLinks [ ] stores all the adjacent upstream links. Its length is numUpLinks. 

numDownLinks is the number of downstream links adjacent to a link. 

downLinks [ ] stores all the adjacent downstream links. Its length is numDownLinks. 

typeUpJunction denotes the type of the upstream junction incident to a link. Here, 
type stands for a linear junction connecting one upstream link and one down- 
stream link, type 1 for a merging junction, type 2 for a diverging junction 
as described by Equation 16.111 type 3 for a diverging junction by Equation 
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16.81 -Equation I6.im and type 4 for a diverging junction by Equation 16.121 For 
different types of junctions, traffic flow models are different, as shown in Section 

typeDown Junction denotes the type of the downstream junction incident to a link. 
The definition of junction types are the same as in typeUpJunction. 

downBoundary is an address pointing to the furthest downstream boundary of a link. 

upBoundary is an address pointing to the furthest upstream boundary of a link. 

In these cells, four fields in dark green are pointers with no physical meaning, two 
fields in cyan are time-dependent quantities, and the rest in blue represent quantities 
that determine network structure and are time invariant. 

As exhibited by the top right box of Figure 16.21 cellType, has the following 
fields, which characterize a cell. 



length is the length of a cell. 

numlane is the number of lanes of a cell. 

typeFD denotes the type of fundamen tal diag r ams i n a cell. Type stands for the 



triangular fundament al diagram (INewelll. 



1993 ) . For other types of fundamen- 
I 1 1 '-^ — — 1 



tal di agrams, refer to ffDel Castillo and Benite^ . I1995bt IKerner and Konhausei , 



1994^ . With the number of lanes considered, we can have the fundamental 



diagram for the cell, 
ff speed is the cell free flow speed, 
jamdensity is the cell jam density of each lane. 
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cridensity is the cell critical density of each lane, at which traffic flow reaches 
capacity. 

capacity is the cell maximum flow-rate of a lane, 
density is the cell total density. 

proportion [ ] is an array of proportions of commodities in a cell. Its length is 
numCommodity of the corresponding link. 

demand is the cell traffic demand during a time interval, as defined in Equation l(j.41 

supply is the cell traffic supply, as defined in Equation 16.51 

next is the address pointing to the upstream cell. 

As in linkType, the dark green field is a pointer, seven blue fields are for time- 
independent quantities, which denote major characteristics of a cell, and four cyan 
fields are for time- dependent quantities. Note that the number of lanes, free flow 
speed, and, therefore, the critical density and capacity may change when accidents 
occur. Besides, if t here are s i gnals on the boundarie s of a ce ll, the demand and supply 



may be restricted (jDaganzo 



1995a; 



Jin and Zhang 



2003d) 



The data structure for cell boundaries, boundaryType, is illustrated by the bot- 
tom right box in Figure lfi.2l and has the following fields: 

flux is the flux through the boundary during a time interval. 

proportion [ ] stores the proportions of all commodities in the total flux. Its length 
is numCommodity of the corresponding link. 

next is the address pointing to the upstream boundary. 
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As we can see, both flux and proportion [ ] are time dependent. 

In the MCKW simulation, cells and boundaries are ordered in the direction op- 
posite to traffic's. However, it is also straightforward to order them in the traffic 
direction. Besides, ordering will not affect computation efficiency significantly. 

6.3.3 Program flow-chart 

The program flow chart in the MCKW simulation platform is shown in Figure 16.31 
The modules in the program are explained below in the same order as they appear 
in the chart. 

1. Create network. Network structure and characteristics are created. That is, 
we provide values for the blue fields of each link and cell in Figure EISl We also 
assign all locations of the pointers. 

2. Initialize traffic. Traffic density and proportions of all commodities are ini- 
tialized for each cell. One typical initialization is to set network empty; i.e., 
traffic density of each cell is zero. 

3. Compute supply /demand. Given traffic density of a cell, we are able to 
compute traffic demand and supply according to Equation 16.41 and Equation 

4. Resolve boundary conditions. Several types of boundary conditions can be 
used. The first and most important type of boundary conditions is conditions 
for origins and destinations. In the MCKW simulation, we use traffic demand, 
specified for all commodities, for origins and traffic supply for destinations. This 
is different from the boundary conditions used in the previous chapters, where 
the Dirichlet, Neumann, and periodic boundary conditions are generally used. 
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Figure 6.3: The program flow chart in the MCKW simulation 
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However, from those boundary conditions, we can easily compute the supply 
and demand according to Equation 16.41 and Equation 16.51 Moreover, these 
boundary conditions are different from the 0/D flow matrix in dynamic loading 
studies in the sense that we cannot determine in-flow, which is also affected 
by current traffic conditions on a link incident to an origin. Therefore, with 
the same pattern in origin demand, we may have totally different 0/D flows. 
This observation suggests a criteria for evaluating the level of service of a road 
network: the amount of time for loading a number of vehicles. In addition, the 
effect of assignment algorithms can be studied through the proportions of origin 
demands. Second, signals and accidents are considered as boundary conditions 
in the MCKW simulation: signals acting at cell boundaries put a con straint on 



suppli e s of the downstream ce lls and demands of the upstream cells (IDaganzo 



■Tin and Zhan 



2nn8d ). and accidents will change the number of lanes 
and free flow speed in a cell. From the discussions above, we can see that the 
influence of incidents and accidents can be studied through imposing different 
boundary conditions. 

5. Compute link flows. From Equation 16. 6t we are able to compute fluxes 
through boundaries inside a link. From the FIFO principle of traffic flow in the 
kinematic wave theories, the proportion of a commodity in fluxes is equal to 
that in density in the upstream cells. 

6. Compute junction flows. From Equation 16.71 Equation 16. 81 - Equation lU.lOl 
Equation 16.111 and Equation 16.121 we can compute fluxes through different 
types of junctions. The proportions of different commodities can be obtained 
from the FIFO property and traffic conservation of each commodity. Equation 
16.211 Since links share junctions, we only need compute junction fluxes for a 
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set of links. For example, for the road network in Figure we only compute 
junction fluxes at the upstream and downstream junctions of links 2 and 5. 

7. Update traffic conditions. Traffic densities of all commodities in a cell can 
be updated by Equation 16.191 However, one has to be careful when computing 
the proportions since total density, as a divider, may be very small. Although 
the proportions may not be accurate for very small densities, it rarely matters. 

6.4 Cumulative flow, travel time, and other prop- 
erties of a road network 

In the MCKW simulation, we keep track of the change of traffic densities of all 
cells and fluxes through all boundaries. Besides, these quantities are specified for 
commodities. In this section, we will discuss how to obtain other traffic information 
from these quantities. 

6.4.1 Cumulative flow and vehicle identity 

Cumulative flow at a boundary Xi_i/2 from time to to t, N{xi^i/2', [to? A)^ the total 
number of vehicles passing the spot during the time interval. If the flux is f*{xi-i/2, s) 
at time s, then we have 



Correspondingly, the discrete cumulative flow, N{xi^i/2] [Jo, J]), which is from time 
steps Jq to J, is deflned as 




(6.23) 



j-i 




(6.24) 



j=Jo 
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where //^^/a during time steps j and j + 1. 

A cu rve of cumulativ e flow versus time is a lso known as a Newell-curve or simply 



N-curve (jDaganzo 



1994^ , since iNewelll (|1993^ developed a simplified version of the 
LWR kinematic wave theory based on this concept. 

From the definition of cumulative flow, we can see that an N-curve is non- 
decreasing in time. Further, it is increasing when passing flow is not zero. 

Although densities and fluxes are quantities at the aggregate level, the MCKW 
model is capable of tracking traffic information at the disaggregate level. This can 
be done also with cumulative flows: a vehicle passing a cell boundary at a time step 
can be labeled by the corresponding cumulative flow. If all cumulative flows are 
synchronized; for example, when the initial traffic in a road network is empty, then 
the same cumulative flow of a commodity refer to the same vehicle. This fact is due 
to the FIFO property in all commodities. ^ 

Therefore, in the MCKW simulation, with curves of cumulative flows as a bridge 
between the aggregate and disaggregate quantities, we are able to keep track of vehicle 
trajectories, accurate to the order of Ax and At, from cumulative flows at all cell 
boundaries. Further, with finer partition of each link, we can obtain more detailed 



information at the disaggregate level. 



6.4.2 Travel time 

For a vehicle, which can be identified by its commodity cumulative flow number under 
FIFO, its travel time across a link or from the origin to the destination can be inferred 
from N-curves. For example, when we know its arrival and departure times to a link 
from the corresponding N-curves, we can easily compute its travel time across the 
link. 

®When type 4 diverge appears, this has to be checked. 




Figure 6.4: Cumulative flows and travel time 

This can be demonstrated in Figure lfj.4[ In this figure, the left curve is the N- 
curve at location xi, and the right curve at X2. These two curves are synchronized in 
the sense that the vehicles between xi and X2 at t = are not counted in N{x2', [0, t]). 
Therefore, from FIFO, we can see that the Nq vehicle on the left N-curve is the same 
as the A^o vehicle on the right N-curve. Then, from the curve, we know that the times 
of the A^o passing xi and X2 are ti and t2 respectively. Thus, its travel time from xi 
to X2 is t2 — ti. 

In Figure 16.41 the left N-curve reaches a maximum at some time and stop in- 
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creasing after that. This means that no flow passes Xi after that time. The right 
N-curve has the same pattern. In such cases, one has to be cautious when computing 
travel time for the last vehicle, identified by the maximum cumulative flow, which 
corresponds to multiple values in time. Rigorously, therefore, the time for a vehicle 
Nq passing a location x, where the N-curve is N{x; [to,t]), can be defined by 

T{No; x) = min{s|when N{x] [to, s]) = Nq}. (6.25) 

s 

Further, the travel time for the A^o vehicle from xi to X2 is 

T{No;[xuX2]) = T{No;x2)-T{No;xi). (6.26) 

With the definition of passing time in Equation 16. 251 at x, the vehicle identity A'o 
has a one-to-one relationship with its passing time T{No]x). Therefore, the passing 
time can be considered as another identity of a vehicle. For a vehicle Nq, if we know 
its passing time at any location in a road network, we then obtain its trajectory. 

From the travel times of individual vehicles, we are able to compute the total 
travel time between two locations, in particular between an 0/D pair, as follows: 

N2 

Ti[N^,N2];[xi,X2]) = TiM;[xuX2]), (6.27) 

M=Ni 

where A^i is the first vehicle and N2 the last. We can see that, in Figure IHiH the total 
travel time is equal to the area between the two N-curves. Then the average travel 
time for each vehicle will be 

I\2 — -/Vi l\2 — jVi 

Moreover, for a road network, we can integrate travel times for all 0/D pairs 
and, therefore, obtain the total travel time and the average travel time for the whole 
road network. These quantities are important indicators of the performance of a road 



CHAPTER 6. 



MULTI-COMMODITY NETWORK TRAFFIC MODEL 157 



network. Besides, we consider the loading time for an amount of flow to be released 
from an origin as another performance indicator. 

Hence, the MCKW simulation platform can be applied to evaluate traffic manage- 
ment and control strategies, such as route assignment and ramp metering algorithms. 
These applications will be discussed in the next chapter. 



6.5 Numerical simulations 

In this section, we investigate the properties of the MCKW simulation model through 
numerical simulations. We will show the evolution of traffic and examine the conver- 
gence of solutions. In this section, the diverge connecting links 2, 3, and 4 is modeled 
by Equation 16. Ill 



6.5.1 Simulation set-up 

For these simulations, the network has the structure as shown in Figure lUTTl In this 
network, links 2, 3, and 5 have the same length, 20 miles, and the length of link 4 is 



40 miles (not drawn to proportion); link 2 has three lanes, and th e other 



two lanes; all links have the same triangular fundamental diagram (jNewell . 



inks has 



19931) : 



Q{a,p) = { -A^- ^g_29) 

^^^Vf{apj - p), ape < p < apj] 

where p is the total density of all lanes, a the number of lanes, the jam density pj=180 

vpmpl, the critical density pc=36 vpmpl, the free flow speed f/=65 mph, the capacity 

of each lane Qc = Pcf/=2340 vphpl, and the corresponding shock wave speed of jam 

traffic is Cj = —pc/{pj — Pc)vf ~ —17 mph. 

Initially, the road network is empty. Boundary conditions are defined as follows. 

Traffic supply at the destination is always 2gc- At origin, traffic demand at the origin 
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is 3qc during [0, 6.0] and zero after that, and the proportion of commodity 0, which 
takes link 3 instead of 4, is always ^ = 70%. 

Links 2, 3, and 5 are partitioned into N cells each, and link 4 into 2N cells, with 
each cell of the same length. Ax = 20/ N miles. The total simulation time of 8.4 
hours is divided into K time steps, with the lengt h of a time step At = 8.4/i^'. In our 



simulations, we set N/K = 1/30. Thus the CFL fjCourant et al.l .l 



larger than VfAt/dx = 0.91, which is valid for Godunov method (jGodunovl . Il959[ ). 



9281) nurnber is no 



6.5.2 Traffic patterns on the road network 

We let N = 400 and K = 12000. Hence Ax = 0.05 miles=80 meters, and At = 0.0007 
hours=2.52 seconds. Here the sizes of road cells and time steps are relatively small, 
in order for us to obtain results closer to those of the kinematic wave theories with 
the Godunov method. 

The contour plots of the solutions are shown in Figure 16. 51 From these figures, 
we can divide the evolution of traffic dynamics on the road network into three stages. 

In the first stage starting from 0, vehicles embark link 2 with the free fiow speed, 
prevail the link in its critical density 3pc, and arrive junction 1 at ti = 20/65 hr. At 
the diverge, junction 1, fiuxes are computed from Equation 16.111 out-fiux of link 2 is 

. 2q, 2q 20 

f2,out = mm{3qc,—,—} = Y(lc, 

which is slightly smaller than the in-fiux of link 2, in-fiux of link 3 is /s = O.7/2 = 2qc, 
which is its capacity, and in-fiux of link 4 is /4 = O.3/2 = |gc, which is less than 
half of its capacity. After ti, two streams of free fiow form on links 3 and 4, and a 
backward travelling shock wave forms on Link 2, and the shock wave speed is 

1 
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Figure 6.5: Contour plots of network traffic flow 

At t2 = ti + 20/vf, the first vehicle on link 3 reaches junction 2, which is a merge. At 
t2, the first vehicle on link 4 is half way back since the length of link 4 is double of 
link 3's. From the merge traffic flow model, Equation 16. 7| we have the in-fiux of link 

5 as 

/5,i„ = min{2gc, 2gc} = 2gc, 

which is also the out-fiux of link 3. After t2, the proportion of commodity on link 
5 is 1, as we can see on the bottom right figure. 

The second stage starts at ts = ti + 40/f/ = 60/65 hr, when the first vehicle 
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on link 4 reaches junction 2. After that, the in-flux of link 5 is still 2gc, but the 
proportion of commodity reduces to 0.5714 since commodity 1 also contributes; on 
link 3, a new state forms at p = 195.4290 vpm, which is over-critical, and a shock 
wave travels upstream at the speed of \cj\ ~ 17 mph; on link 4, p = 30.8571 vpm, 
which is under-critical. At ^4 = ta + 20/|cj| = 140/65 hr, the back-traveling shock on 
link 3 hits junction 1, and the traffic supply on link 3 is reduced. Therefore, the out- 
flux of link 2 is further reduced, and link 2 becomes more congested, as shown in the 
top left picture. This also reduces traffic flow on link 4, and the reduced flow reaches 
junction 2 at ts = t4 + 40/f/ = 180/65 hr. After ^4, link 3 becomes less congested, 
and a rarefaction wave travels backward on it at \cj\ ~ 17 mph; the proportion of 
commodity on link 5 gets higher. From the bottom middle flgure, we can see that at 
= t4 + 20/|fj|, traffic density on link 4 swings back a little due to the back traveling 
rarefaction on link 3. This shift is transported to junction 2 at = ^5 + 40/f/ and 
oscillates traffic density on link 3 and the proportion of commodity on link 5. 

The third stage starts at ty = 6, when traffic demand from origin subsides to 
zero. After that, a shock forms on link 2 and travels forward, and propagates to link 
4 and link 3. On link 4 the shock travels at fj, and on link 3 it travels slower. This 
is why the proportion of commodity on link 5 becomes 1 before it is emptied. 

This simulation indicates that oscillation of traffic conditions can be caused by 
network merges and diverges even the initial and boundary conditions are very nice. 
The traffic flow pattern on this road network suggests that, if we keep the same 
demand from the origin, an equilibrium state will be reached after some time. This 
equilibrium state will be further investigated in the following chapter. 

The traffic patterns in Figure 1^31 can be partially observed from the top left flgure 
in Figure inini where the thicker four curves give cumulative flows for commodity 0, 
the thinner for commodity 1, and solid, dashed, dotted, and dash-dot curves are for 
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Total Number of Vehicles 


TTT 


ATT 


Commodity 


23,8585 


4.7291x10^ 


1.9822 


Commodity 1 


10,2251 


1.7372x10^ 


1.6989 



Table 6.1: Total travel time (TTT) and average travel time (ATT) for two commodi- 
ties 

cumulative flows at junction 0, 1, 2, and 3, respectively. In the bottom left figure, the 
solid lines are cumulative flows of commodity at origin/destination, i.e., junction 
and 3, the dashed curve is the number of commodity vehicles in the network 
at a time, the dashed line shows the average number, the dotted curve is the travel 
time of a commodity vehicle identified by its cumulative flow, and the dotted line 
is the average travel time. The bottom right figure has the same curves and lines 
as the bottom left figure, except that they are for commodity 1. Here travel time of 
individual vehicle is computed by Equation 16.261 total travel time by Equation I6.27[ 
and average travel time by Equation 16.281 Although the formulas were developed 
for link travel times, it is valid for a network as long as vehicles of each commodity 
observe FIFO principle. Total (TTT) and average (ATT) travel times are listed in 
Table 16.11 (unit=hours) . 

Thus, vehicles that take link 4 on average use shorter time, which is still longer 
than the free flow travel time, 80/65 hr. Obviously, if all vehicles at origin decide to 
take link 3, the travel time, 60/65 will be the shortest possible travel time between the 
origin and destination. Therefore, this assignment fraction, 70%, is not an optimum 
one. More detailed analysis of the influence of the assignment fraction on travel times 
will be engaged in the following chapter. 
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Commodities and 1 at All Junctions 




Commodity at Origin/Destination 



Commodity 1 at Origin/Destination 




2 4 6 

time (hours) 




2 4 

time (hours) 



Figure 6.6: N-curves and travel times of each commodity in the road network 



6.5.3 Convergence of the MCKW simulation model 

In this subsection, we study the convergence of the MCKW simulation platform with 
increasing number of cells. Here we use the same road network, initial and boundary 
conditions as in the preceding subsection. 

As we show in the previous chapters, the discrete forms of kinematic wave theories 
converge to the i r cont inuous counterparts as we partition each link into finer cells 



( Jin and Zhang 



2nn8d ). In those studies, we generally check convergence in terms of 



traffic densities. Here we intend to show convergence in average travel times of both 
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Commodity 


AT c\r\r\ 

N=zOO 


AT A C\C\ 

N=400 


AT Cir\r\ 

N=800 


AT 1 r\r\ 

N = 1600 


AT oor\r\ 

N=3200 


Ai i 


1.98189893 


1.98215215 


1.98227240 


1.98234941 


1.98239377 


Error [10 J 




0.2532 


0.1202 


0.07700 


r\ r\ A A A 

0.0444 


Rate 






1.074 


0.d4o0 


0.795o 


Commodity 1 


N=200 


N=400 


N=800 


N=1600 


3200 


ATT 


1.69922958 


1.69892887 


1.69877593 


1.69871236 


1.69868722 


Error [10~^] 




-0.3007 


-0.1529 


-0.0636 


-0.0251 


Rate 






0.9755 


1.2664 


1.3384 



Table 6.2: Convergence rates for the MCKW simulation platform 
commodities. 

Denoting the average travel time of a commodity, T, as a function of the number 
of cells, A^; i.e., T = T{N), we can define the relative error, from A^ to 2A^, by 

^2N-N ^ \T{2N)-T{N)\. (6.30) 
Then a convergence rate is computed by 

= log2(^i]v=27v)- (6-31) 

The convergence rates of the average travel times are given in Table 16.21 

From the table, we can see that average travel times are also convergent in first 
order. Note that this convergence is different from the aforementioned traffic condi- 
tions converging to certain equilibrium state. Moreover, we can see that the results 
with A^ = 200 is already accurate enough in this case. Since the computation time of 
the MCKW simulation platform is quadrupled when A^ is doubled, in later simula- 
tions, we use Ax = 0.1 mile=160 meters and At = 0.0014 hours= 5.04 seconds with 
the same simulation period. 



CHAPTER 6. 



MULTI-COMMODITY NETWORK TRAFFIC MODEL 164 



6.6 Discussions 

In this chapter, we proposed the Mult i- Commodity Kinematic Wave (MCKW) sim- 
ulation model. In this simulation model, we integrated the kinematic wave theories 
studied in the previous three chapters to form the foundation of the algorithms, care- 
fully discussed commodity-based kinematic wave theories, and presented the data 
structure and program structure for implementation. We further demonstrated how 
to obtain cumulative flows and travel times from outputs of the MCKW simulation 
model. Simulations show that numerical results converge to FIFO solutions although 
the FIFO condition is not strictly enforced in the discrete form of commodity-based 
kinematic wave theories. 

Different from many existing simulation packages, where traffic is tracked down 
to vehicle level, the MCKW simulation concerns traffic conditions down to commodity 
level. The simulation model is designed for handling very large road networks and can 
be applied in studies of Intelligent Transportation Systems, such as dynamic traffic 
assignment, dynamic 0/D estimation, and so on. However, as pointed out earlier, the 
effects of "departure from FIFO" should be carefully considered in these applications. 

In the future, the MCKW simulation model can be enhanced in three aspects . 



Theoretically, vehicle types and special l anes can be inc o rporat ed (jDaganzo 



and nonequilibrium continuum models (jjin and Zhang 



3 



mm, 



2003dl ) may also be inte- 



grated. Numerically, parallel algorithms can be applied to improve computational 
speed since traffic conditions on different links can be updated simultaneously, and 
consumption of computer memory will be checked. Finally, for different applications, 
we also plan to design different input/out interfaces. For example, the network struc- 
ture can be imported from CIS (Geographic Information System) data, and boundary 
conditions and output can be manipulated for different applications. 



Chapter 7 

Studies of network vehicular traffic 
with kinematic wave simulations 



7.1 Introduction 

As discussed in the preceding chapter, the Mult i- Commodity Kinematic Wave (MCKW) 
simulation model is based on the kinematic wave theories of traffic dynamics at various 
road network components, including link bottlenecks, merges, and diverges. After ve- 
hicles are differentiated into commodities by their paths or origin/destination (0/D) 
pairs, we can then keep track of the evolution of traffic densities of all commodi- 
ties on a road network over time. Further, information of individual vehicles can be 
obtained through cumulative flows. Therefore, co mpared to some o ther simulation 



model s based on the kinematic wave theory (e.g.. 



Daganzd (|l995a[ l: 



Vaughan et al. 



( 1984( )). the MCKW simulation is a pure macroscopic simulation program. As such, 
the MCKW simulation model is expected to be less costly in computation but ade- 
quate for many applications in Advanced Transportation Information Systems (ATIS) 
and Advance Transportation Management Systems (ATMS). 
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In this chapter, we will explore traffic dynamics in a road network with the 
MCKW simulation model and its implications in applications. First in Section \7.2\ 
we are interested in equilibrium states of a road network under a given 0/D demand 
pattern and the influence of different assignment strategies. Then in Section 17.31 
we study the formation mechanism of periodic oscillations in a network and their 
properties. 

These simulations are studied on a small, demonstration network, where traffic 
dynamics at diverges obey Equation 16.111 The major purpose of these studies is 
to show the capabilities of the MCKW simulation model. Moreover, these studies 
can be considered as initial steps to the understanding of more complicated traffic 
phenomena in a large-scale road network and practical applications. 

7.2 Equilibrium states of a road network and pre- 
liminary examination of traffic assignment 

In Subsection 16.5.21 we show that traffic flow on a road network can reach a certain 
pattern after sufficient amount of time. We call these time-invariant traffic patterns 
as equilibrium states. 

For the equilibrium states, we are interested in the following questions: 1) how 
are the equilibrium states related to the boundary conditions? 2) what is the effect of 
changing proportions in origin demand? 3) do these equilibrium states have anything 
to do with the so-called User Equilibrium States? 
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Figure 7.1: Network for studying equilibrium state and traffic assignment 



7.2.1 The simulated network 

In this section, we study a demonstration network shown in Figure 17.11 In this 
network, links 2, 3, and 5 have the same length, 20 miles, and the length of link 4 is 
40 miles; link 2 has three lanes, and the other links has two lanes; all links have the 
same fundamental diagram: 



Q{a,p) 



VfP, < p < apc, 

^Vf{apj - p), ape < p < apj- 

where p is the total density of all lanes, a the number of lanes, the jam density 



Pj=180 vpmpl, the critical density Pc=36 vpmpl, the free flow speed f ^=65 mph, and 
the capacity of each lane qc = pcf/=2340 vphpl. 

Initially, the road network is empty. Traffic supply at the destination is always 
the capacity of two lanes. For studying equilibrium states, we assume traffic demand 
at the origin is always the capacity of link 2, 3qc- 

Links 2, 3, and 5 are partitioned into 200 cells each, and link 4 into 400 cells. 
Therefore, all cells have the same length, 0.1 miles. The total simulation time of 8.4 
hours is divided into 6000 time step s, with the length of a time step At = 0.0014 



hours=5.04 seconds. Thus the CFL ( Courant et al 



VfAt/dx = 0.91, which is valid for Godunov method (|Godunov 



1928ri number is n o bigger than 



19591 ). As we show 



in subsection I6.5.3| the simulation results are numerically convergent. 



CHAPTER 7. 



STUDIES OF NETWORK VEHICULAR TRAFFIC 



168 



In the following simulations, we will study the equilibrium states and the perfor- 
mance of the road network for a proportion of commodity 0, ^. 

7.2.2 Equilibrium states 

When ^ = 0.5, and commodity 1 vehicles have the same proportion, 1 — ^ = 0.5, 
solutions of traffic dynamics on the road network are shown in Figure 17.21 The traffic 
pattern evolves as follows. From to = to ti = 20/65 hr ^ , the first vehicle traverses 
link 2 and reaches junction 1, where fluxes can be computed from Equation 16.111 as 

f2,out = min{3gc,2gc/0.5,2gc/0.5} = 3gc, 

/3,m = I-SQ'c, 
/4,m — 1-5Q'c; 

where subscript in denotes in-flux, and out out-flux. After ti, links 2, 3, and 4 all 
carry free flow. At ^2 = ^1 + 20/65, traffic on link 3 arrives at the merge, junction 2, 
where from Equation 16.71 fluxes are the following: 





= min 




= 2gc, 


f4,out 


= 0. 



At ^3 = ^2 + 20/65, vehicles on link 4 also reaches junction 2, and fluxes become 

/5,i„ = min{2gc + 2gc, 2gc} = 2gc, 

fA,out ~ Qc- 



^Hereafter, the unit of time is always hour except if otherwise mentioned. 
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Thus, after ts, shock waves form on both links 3 and 4 and travel upstream in the 
speed of —Vf/A. At ti = + 80/65, the shock wave on link 3 reaches junction 1, and 
fluxes become 

f2,out = min{3gc,gc/0.5,2gc/0.5} = 2gc, 

fi,in ~ Qc- 

After this, link 3 equilibrates at = qc and ps = 1.2pj, a shock travels backward on 
link 2 at f//4, and a shock travels forward at the free flow speed. At = t4 + 20/65, 
the two shock waves on link 4 meet at the middle of the link and form a stable zero- 
speed shock wave, where the upstream half is at 0.2pj and the downstream half at 
1.2pj. After ^6=^4 + 80/65, traffic on link 2 is uniformly at q2 = 2qc and p2 = lApj. 

Therefore, the whole road network reaches an equilibrium state after tg- In the 
equilibrium state, links 2 and 5 both carry flow-rate 2qc, which is the capacity of 
the road network, links 3 and 4 carry flow rate q^. corresponding to the proportion 
of each commodity. Besides, links 2 and 3 are congested, and link 4 has a zero- 
speed shock, which, however, is unstable in the sense that a small oscillation in the 
upstream or downstream flow will make the zero-shock disappear. The instability of 
this equilibrium state on the road network can also been seen later when we consider 
^ away from 0.5. 

When ^ = 0.6, the contour plots of traffic densities on four links are given in 
Figure 17.31 The traffic evolution pattern is similar to that in subsection 16.5.21 Here 
we directly go into the discussion of the equilibrium state, in which, from Figure [721 
we can see that all links carry uniform flows. The equilibrium flow-rate is determined 
by the network bottleneck, link 5, and is 2qc. From observations for ^ = 0.5, the equi- 
librium density on link 5 is always OApj, and traffic density on link 2 is lApj, which 
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T otal De nsi ty oj^i'i^gp: :!l^|^j;^'^T CommodityOonLinkS 




Figure 7.2: Solutions when ^ = 0.5 

is brought by a back-travehng shock wave. Then, determined by traffic proportions, 
the flow-rate on links 3 and 4 are 1.2gc and O.Sqc, respectively. The remaining task is 
to determine the traffic densities on links 3 and 4. 

We denote traffic demand and supply on these links as D^, S3, D4, and 5*4. Then 
S3) — (1.2gc, 25c) when link 3 is under-critical, and (D3, S3) — (2gc, l-2gc) other- 
wise. Similarly, (£)4, S4) — {O.Sqc, when link 4 is under-critical, and (D4, S^) — 
{2qc, O.Sqc) otherwise. Since at junction 1 

f2,out = min{3gc, 5*3/0.6, 5*4/0.4} = 2gc, 
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Figure 7.3: Solutions when ^ = 0.6 

links 3 and 4 cannot be under-critical at the same time. Moreover, since at junction 
2 

f3,out = -^^^2qc = 1.2gc, 
•^4,o«i = D^+D4, = 0.8gc, 

links 4 cannot be over-critical at the same time. Therefore, the only possible case is 
that link 3 is over-critical and link 4 under-critical. Thus, the equilibrium densities 
on links 3 and 4 are ps = l.OApj and = 0.32pj. The solutions in Figure 17751 support 
these solutions. 
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< ^ < 0.5 


i = 0.5 


0.5 < ^ < 1 


Link 2 


(1.4,2) 


(1.4,2) 


(1.4,2) 


Link 3 


(0.4^,20 


(1.2,1) 


(2-1.6^,20 


Link 4 


(0.4 + 1.6^,2(1-0) 


(0.2,1) and (1.2,1) 


(0.4(1-0,2(1-0) 


Link 5 


(0.4,2) 


(0.4,2) 


(0.4,2) 



Table 7.1: Equilibrium density and flow-rate v.s. ^ 



However, D3 = 2qc and D4 = O.Sqc do not satisfy Equation 17.1 



there is an intermediate state at the downstream boundary of link 4 (jjin et al 



exactly. In fact . 



20021 ). 



At the intermediate state, also under-critical, the demand is D4 = |gc, which satisfies 
Equation 17.11 and p4 = ^Pj- This intermediate state, theoretically, only exists at a 
point, but is stable since its in- flux and out-flux is equal. T his kind of intermediate 



state unlikely exists on a single link (jJin and Zhane 



2003bl ). 



When ^ = 0.4, the solutions are shown in Figure 17.41 Although the evolution 
process is different from that for C, = 0.6, the equilibrium state is exactly the same as 
before if links 3 and 4 are switched. 

For a general proportion G [0, 1], the equilibrium states are listed in Table 1711 
In this table, we omit the tenuous (unstable) intermediate states, on link 3 when 
< ^ < 0.5 and on link 4 when 0.5 < ^ < 1. 



7.2.3 Travel times at equilibrium states 

For equilibrium states of different proportion C,, Table 17.21 shows travel speed and 
travel time on each link as well as average travel times of commodity (ATTq), 
commodity 1 (ATTi), and all commodities {ATT = ^ATTq + (1 — O^^^i)- In the 
table, travel times during the establishment of equilibrium states are not considered. 
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Figure 7.4: Solutions when ^ = 0.4 



The relationship between these travel times and proportion ^ is illustrated in Figure 
1731 

From the table and figure, we can see the following properties of travel times. 
First, the average travel time of all commodities, which can be considered the perfor- 
mance function of the whole road network, consists of two pieces of lines and jumps 
when ^ = 0.5. This discontinuity asserts the fact that the equilibrium state when 
^ = 0.5 is not stable. Third, the minimum average travel time ATT attains its min- 
imum when ^ = 1; i.e., when all vehicles take the shorter path. Third, the travel 
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[v/Vf, I /hours) 


< ^ < 0.5 


^ = 0.5 


0.5 < ^ < 1 


Link 2 


(|, 1.0769) 


(|, 1.0769) 


(|, 1.0769) 


Link 6 


(1, 0.3077) 


(i, 1.8462) 


/ f 20-16f\ 
'^5-45' 135 


Link 4 


/ 1-5 8(l+45)\ 
vi+45' 13(1-0 


(1 or g, 2.1538) 


(1, 0.6154) 


Link 5 


(1,0.3077) 


(1,0.3077) 


(1,0.3077) 


ATTo 


1.6923 


3.2308 


2^^*^^ + 1.3846 


ATTi 


|l±Sl+i,3846 


3.5384 


2 


ATT 


2 + 2.7692^ 


3.3846 


3.5385- 1.8462^ 



Table 7.2: Equilibrium speed and travel times v.s. ^ 

times of commodities and 1 are equal when C, = 5/6. However, ATT is not at its 
minimum when ^ = 5/6. 



7.2.4 Discussions 

In the previous subsections, we examine equilibrium states for the road network in 
Figure 17.11 and its performance at these states against the proportion ^. The de- 
termination of ^ is really part of a traffic assignment process. The refore, thi s study 




provides another angle for checking the user-equilibrium property (jWardrop . 
For this network, it is shown to have two user equilibrium states, one at ^ = 5/6 
and the other at ^ = 1, with the latter has smaller path travel times. This is a clear 
evidence that the inclusion of physical queues destroys the uniqueness property of 
traffic ass ignment with ot her commonly used models, in which queues do not take 



up space (jWardrop 



19521 ). Note that, in our discussions, the performance function 



is the average path travel time, which is different from the link performance function 



used in static traffic assignment studies (|Sheffi . 



19841) . 
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Figure 7.5: Travel times at equilibrium states v.s. proportion ^ 



Although the connection between equilibrium states from a dynamic point of 
view and the user-equilibrium is still premature, these studies shed light on solving 
dynamic traffic assignment problems, in whic h the choice of a proper link performance 



function is still in debate (jDaganzg 



1995b| ). From the MCKW simulation, we can 



obtain average origin/destination (0/D) travel times under an assignment strategy 
and use these travel times to evaluate the strategy. To make the performance function 
of an 0/D flow more accurate, one can include the loading time at the origin, which 
is defined as the average waiting time before a vehicle enters the road network. If 
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departure (or loading) flow-rate at the origin is /(t) (t G [0,T]), then the average 
loading time is 



LT 



(7.2) 



In the examples above, since the loading flow-rate is the same for all C,, there is no 
difference between the loading times. 



7.3 The formation and structure of periodic oscil- 
lations in the kinematic wave model of road 



networks 



In fjjin and Zhane . 



2003a ) it is shown that there are three types of basic kinematic 



waves on a road link: shock (decelerating) waves, rarefaction (accelerating) waves. 



and transition (standing) waves. Further, for a single merge or diverge, we sti l 



2003c; 



Jin et al. 



have 



2002; 



these three types of ^ waves o n each branch ()Jin and Zhane . 
Jin and Zhana . i2002t . i2001a) . In this section, we show with MCKW simulation an 
interesting type of solutions, periodic oscillations, which can be observed in real traffic. 
We show that these solutions can exist in a small road network with a diverge and a 
merge and check their formation and structure against network characteristics. 



7.3.1 Network for studying periodic solutions 

In this section, we study periodic solutions on a road network shown in Figure 17.61 
In this network, lengths of links 2, 3, 4, and 5 are L2 = 10 miles, L3 = 1 mile, L4 = 2 
miles, and L5 = 1 mile, respectively; the number of lanes of these links are 3, 1, 2, 
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Figure 7.6: Network for studying periodic oscillations 



and 2, correspondingly; all links have the same fundamental diagram for each lane: 



Q(a,p) 



Pj-Pc 



< p < apc] 
Vf{apj - p), apc< p < apj] 



where p is the total density of all lanes, a the number of lanes, the jam density 
Pj = 180 vpmpl, the critical density Pc=36 vpmpl, the free flow speed Vf=65 mph, and 
the capacity of each lane qc = pcf/=2340 vphpl. 

Initially, the road network is empty. Traffic supply at the destination is always 
capacity of two lanes. Here, traffic demand at the origin is always the capacity of link 
2, 3g,. 

Links 2, 3, 4, and 5 are partitioned into 800, 80, 160, and 80 cells respectively, 
with the length of all cells as 0.0125 miles. The total simulation time of 1.4 hours 
is divided into 8000 time steps, wit h the length o f a tim e step At = 1.75 x 10""^ 



hours=0.63 seconds. Thus the CFL (^Courant et al 



1928ri nurnber is no bigger than 



Godunov 



19591 ) . As we show in 



VfAt/dx = 0.91, which is valid for Godunov method 
subsection 16.5.31 the simulation is numerically convergent. Therefore, with this very 
fine partition of links, we are able to obtain results closer to theoretical solutions. 
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Link 3 


under-critical 


over-critical 


Link 4 


under-critical 


over-critical 


^3 














< Qc 






< 2g, 


2gc 



Table 7.3: Constraints on equilibrium states 



7.3.2 Periodic oscillations 

As in the preceding section, we also consider the equilibrium states for a proportion 
of commodity 0, ^. In equilibrium states, at junction 1, we have 

f2,out = min{3gc,5'3/^, 5*4/(1 -0}, 

f3,in = if2,outi i^-^) 
fiM = ~ C) f2,outi 

and at junction 2 

/5,i„ = m.m{Dz + L'4, 2gc}, 

fifiut = D^^D^ f5,ini 

where Di and Si are traffic demand and supply of link i {i = 2,3,4,5) respectively, 
-D2 = Sqc, and 5*5 = 2gc- Moreover, we have 

/3,m fsfiutj 
fi,in fi,out) 

and constraints defined in Table [721 where ^3 and are fiow-rates in the upstream 
cells of links 3 and 4. 

Since links 3 and 4 are initially empty, = Qc and 6*4 = 2gc- Therefore, when 
^ = or 0.5 < ^ < 1, f2,out < 2gc according to Equation 17.31 Note that the capacity 
of link 5 is 2gc- Under this circumstance, solutions are trivial. When < ^ < 0.5, 
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from Equation 17.31 we have f2,out > 2gc- This means that at least one of links 3 and 
4 has to be over-critical. When both links are over-critical, we have from Equation 
QEquation O and Table O that ^ = |. When < ^ < i, link 4 is congested in 
equilibrium state, and the solution pattern is similar to Figure 17.31 and Figure 17.41 
Refer to subsection 17.2.21 for an analysis of the equilibrium states. 

For the network in Figure 171)1 interesting solutions occur for | < < |) when 
link 3 is over-critical and link 4 under-critical. When ^ = 0.45, contour plots of traffic 
densities on the four links are shown in Figure 17.71 This figure shows the formation of 
periodic solutions as follows: (i) During to = and ti = L^/vf, capacity fiow travels 
forward on link 2. When the first vehicle reaches junction 1 at ti, we obtain from 
Equation O that /2,o«t = gc/0.45, fs^in = qc, and /4,j„ = ^g^- (h) Thus, after ti, a 
back-traveling shock wave forms on link 2, capacity fiow travels on link 3, and free fiow 
travels on link 4. At t2 = ti+L^/vf, the first vehicle on link 3 reaches junction 2. Since 
the first commodity-1 vehicle is still half way on link 4, from Equation 17.41 we obtain 
= Qc, and fi^out = 0. Therefore, there is no change in traffic patterns 
on links 2 or 3 at t2- (hi) AT ^3=^1 + L^/vf, vehicles on link 4 reach junction 2. 
From Equation 01 we have f^^in = 2qc, fz^out = 0.9gc, and f4^^out = l-l^c- (iv) After 
ta, a back-traveling shock forms on link 3 and the over-critical traffic state, whose 
fiow-rate is 0.9gc, propagates upstream. On link 4, there exists a fiimsy intermediate 
area, which is still under-critical. Thus there is no backward wave on link 4. The 
shock wave on link 3 travels at the speed jVf and hits junction 1 at ^4 = ^3 + 4^3/1'^. 
(v) At t4, we have 5*3 = 0.9gc, which yields from Equation 17.31 that f2,out = 2gc, 
f3,in = 0.9qc, and f^^in = l.lgc- Therefore, a forward shock wave forms on link 4, whose 
upstream fiow-rate is l.lgc and downstream fiow-rate ^qc- (vi) At ^5=^4 + L/^/vf, 
the new traffic state reaches junction 2. Consequently, f^^out = 2/2. Igc > 0.9gc, and 
fi,out = 2.2/2. Igc < l-l^'c- Thus the downstream of link 3 becomes less congested, and 



CHAPTER 7. 



STUDIES OF NETWORK VEHICULAR TRAFFIC 



180 



TotRl Densily on Link 2 Density of Commodity D on Linit 3 




Figure 7.7: Contour plots of periodic oscillations 

a rarefaction wave travels upstream. Note that there is no back-traveling wave on link 
4. (vii) At tQ — t^+ALs/vf, the rarefaction wave reaches junction 1 and increases out- 
flux of link 2, f2,out = 2.1164gc, and f^^in = 1.164gc > 1-lg'c- (viii) When the new flow 
on link 4 reaches junction 2 ai — Iq + L^/vf, link 3 becomes more congested again. 
After another time period, {4:L3+L4)/vf, link 4 will go back to lower congestion. This 
process repeats in a periodic manner with the period T — 2(4L3 -|- L4)/v/=0.1846 
hours. After a number of periods, traffic conditions on links 2, 3, and 4 reach a stable 
pattern, a period of which is explained as follows. 
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Figure 7.8: Density, flow-rate, and speed at the end of link 2, start of link 3, and 
start of link 4 

The formation of periodic solutions can also be observed at certain locations, 
such as the end of link 2, start of link 3, and start of link 4, as shown in Equation 
17.81 From Figure FTTl we can see that, most of the time, wave speeds on links 2, 3, 
and 4 are —\vf, —\vf, and Vf, respectively. Thus, from curves in Figure FTHl we can 
derive traffic pattern at other locations. 

In a period, traffic conditions evolve in the following five stages, (i) Assuming the 
period starts at to = when the downstream of link 3 becomes less congested at ps. 
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which is always over-critical, at to, link 4 has a smaller density /94. (ii) At ti = 4L^/vf, 
traffic density on link 4 starts to increase to p^. (iii) After t2 = ti + L^/vf, traffic 
density on link 3 increased to p^. This finishes the half period, (iv) At = t2+4L^/vf, 
traffic density on link 3 decreases to ps. (v) At t4 = T, traffic density on link 3 gets 
back to ps. 

7.3.3 The structure of periodic solutions 

In half a period, solutions at the end of link 2, start of link 3, and start of link 4 in 
the (p, g)-plane are shown in Figure ITIH in which lines are for fundamental diagrams 
and dots for solution data. From the figure, we can see that solutions in each period 
is highly symmetric, with transition layer evenly distributed in the (p, g)-plane. In 
half a period, p decreases and approximately satisfies the following equation: 



where p and p are the maximum and minimum of p respectively, and a is a constant 
to be determined, and the other half satisfies 



dp 



(7.6) 



d t 



dp 



a(p-p)(p-p), 



(7.7) 



d t 



which is increasing. 



In a period, the solutions can be approximated by 




when t e [0, |), 
when t G [f ,T). 



(7.8) 



The corresponding v and q are solved by 




when t G [0, |^)^ 
when t G [f-, T) 



(7.9) 
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Figure 7.9: Periodic solutions at the end of link 2, start of link 3, and start of link 4 
in the (p, g)-plane 



q + {q-q)/ {^ + exp [-a(p -p){t- T/4)]} when t e [0, f ), 
q + {q-q)/ {l + exp[a(p-p)(t-3T/4)]} whent e [f,T). 



(7.10) 



As shown in Figure 17.101 density, flow-rate, and travel speed at the end of link 2 can 
be well approximated with a = 500, and the average of q2 in a period is equal to 
OApjVf ~ 2gc- Therefore, this pattern is stable. 

Solutions at other locations on link 2, or locations on links 3 and 4 share simi- 
lar structure as in Equation 17.81 -Equation I7.10t but the phase or amplitude may be 
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Figure 7.10: Structure of periodic solutions on link 2 

different. 

7.3.4 Discussions 

From Figure 17. 7[ we can see that we can observe several periods at the same time 
only on link 2, although traffic on links 3 and 4 are also periodic. Therefore, these 
periodic waves can be considered to exist on the link upstream to a diverge and a 
merge when the proportion satisfies a certain condition. The special road structure, 
i.e., a diverge and a merge, can be considered as the generator of periodic solutions. 
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We can also see that the period is determined by the length of links between the 
diverge and merge. In addition, possible approximation of the periodic solutions are 
provided. 

This study shows the effect of network structure on traffic dynamics. We expect 
more un-revealed solutions patterns exist in more complicated road networks. In the 
future, it will also be interesting to study the propagation of periodic solutions in road 
networks. It is not hard to check that when traffic demand decreases, the periodic 
solutions will disappear. 

Since periodic solutions have important negative impact on driver's behavior and 
vehicles' emission, this study will help to resolve such kind of solutions. 

7.4 Conclusions 

In this chapter, with the MCKW simulation model, we studied equilibrium states in 
a demonstration network and oscillations caused by a special network structure. 

In Section \7.2\ we showed that traffic dyanmics on a network approach equi- 
librium states with constant boundary conditions. Against different combination of 
traffic commodities, solutions of equilibrium states and the performance of the road 
network were carefully examined. We further discussed possible relationship between 
these equilibrium states and the user-equilibrium. Traffic equilibrium states in a road 
network can be used pursuing important concepts in many network traffic studies, 
such as traffic assignment, origin demand estimation, etc. 

In Section 17.31 we revealed the existence of periodic solutions in a road network. 
We showed the structure of periodic solutions, whose period is solely dependent on 
the network structure. This study shows that, in a road network, special network 
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structure may cause stop and go oscillatory waves. ^ The periodic oscillations are 
in fact a new type of kinematic waves initiated from simple, jump initial conditions. 
From this study, we expect more complicated kinematic waves when more complicated 
networks are considered. 

Studies in this chapter show that, with the MCKW simulation, we are able to 
study more complicated kinematic waves in a road network. This is important for 
understanding traffic dynamics as demonstrated that complicated traffic phenomena 
can be caused by network topology and 0/D demand pattern. However, other factors 
may also contribute to stop and go waves in a realistic network. 



^Another type of veh i cle clusters are s hown to exist in higher-order models 
i Kerner and Konhausen. ll994Hjin and Zhang. l2003af) . 



Chapter 8 



Conclusions 



8.1 Summary 



The ori ginal kinematic wave theory, als o known as the Lighthill-Whitham- Richards 



theory (jLighthill and Whitham . 



1955b; 



Richards 



19561 ). explains traffic dynamics 



of single origin/destination traffic on a linear road with kinematic waves, including 
decelerating (shock) waves and accelerating (expansion or rarefaction) waves. In this 
dissertation, we study traffic phenomena in a road network within the framework of 
kinematic waves. We have theoretically and numerically investigate traffic dynamics 
for link inhomogeneities, junctions, and mixed-type vehicles. Furthermore, we have 
developed a simulation platform of multi-commodity network traffic and studied some 
network traffic phenomena. 

The major results of this dissertation are as follows. 

1. In Chapter 2, we reformulate the Lighthill-Whitham- Richards model into a non- 
linear resonant system for link bottlenecks such as lane-drops. We then show 
that there is an additional type of kinematic waves, namely, standing (tran- 
sition) waves. A standing wave always stays at the inhomogeneity spot, and 
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traffic on its both sides has the same flow-rate but different velocities. We 
further solve the Riemann problem by ten combinations of shock, rarefaction, 
and transition waves. The wave solutions are summarized and shown to be 
consistent with solutions by the supply-demand method. Using numerical sim- 
ulations, we have also shown how traffic queues up in the region upstream to a 
link bottleneck. 

2. In Chapter 3, we take a closer look at existing kinematic wave models for merg- 
ing traffic within the supply-demand framework. In particular, we study the 
distribution scheme, which helps to uniquely determine traffic flow from each 
upstream branch to the downstream link. Further, we propose a "fairness" con- 
dition, under which the out-fiow from each upstream branch is proportional to 
its demand. This condition yields the simplest distribution scheme, which is the 
only valid distribution independent of downstream conditions. We also show 
that it can capture the key characteristics of a merge and leads to a merge model 
that is computationally efficient and easy to calibrate. Furthermore, we show 
that the new distribution scheme is well-defined and qualitatively sound. With 
numerical simulations, we demonstrate that the distribution scheme produce 
convergent solutions. 

3. In Chapter 4, wc propose a new kinematic wave model for highway diverges. 
Based on the assumption that, in a short time interval, diverging traffic flows 
are independent of each other, dynamics of traffic to a downstream link can be 
described by a system of non-strict conservation laws, or a nonlinear resonant 
system, whose Riemann problem can be solved by seven types of wave combi- 
nations. These waves are called as instantaneous kinematic waves. Further, we 
show that this model is equivalent to a supply-demand method with modified 
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definitions of traffic demands of diverging traffic. With numerical simulations, 
we show that this model is consistent with some existing rules on diverging 
traffic. 

4. In Chapter 5, we study traffic dynamics with mixed-type vehicles, for which we 
find an additional family of kinematic waves, contact waves. We carefully study 
the wave solutions of the Riemann problem and develop a Godunov method for 
solving the model. Using simulations, we demonstrate how mixed traffic evolves 
on a road link and show that the First-In-First-Out (FIFO) principle is always 
observed. 

5. In Chapter 6, we propose a multi-commodity kinematic wave (MCKW) model of 
network traffic ffow. In this model, we apply the models studied in the previous 
chapters for different network components and classify vehicles into a number 
of commodities according to their paths. The proportions of commodities in a 
road cell are updated based on the fact that traffic is anisotropic. We then show 
that this model yields solutions, whose departure from the FIFO principle is in 
the order of a time interval. This model is macroscopic, but is able to provide 
individual vehicles' trajectories and travel times by using cumulative curve. We 
also propose an implementation of the MCKW simulation and carefully design 
the data structure for network topology, traffic characteristics, and simulation 
algorithms. Our numerical simulations demonstrate that the numerical results 
converge to FIFO solutions. 

6. In Chapter 7, we first study equilibrium states in a road network with a single 
origin-destination (0/D) pair and two routes. We demonstrate the formation 
of an equilibrium state. After examining the relationship between equilibrium 
states and the distribution of vehicles to different routes, we find that multiple 
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equilibrium status may exist for the same 0/D flow but different route distribu- 
tions. We then show the formation of periodic osciUations for certain network 
structure and route distributions. The periodic oscillations are in fact a new 
type of kinematic wave in a road network. We then discuss their structure and 
properties. 

The kinematic wave models of network vehicular traffic studied in this disserta- 
tion are theoretically rigor, numerically reliable, and computational efficient. These 
studies help better understand traffic dynamics in a road network, in particular the 
formation and propagation of traffic congestion, and establish a solid foundation for 
applications in traffic control and management. 

8.2 Future research directions 

There are three directions in which we could extend this dissertation research: further 
investigation of the kinematic wave theories, enhancement of the MCKW simulation 
model, and applications of the MCKW simulation model to solve real-world problems. 

8.2.1 Further investigations of the kinematic wave theories 

In studies of inhomogeneous links (Chapter 2) and diverges (Chapter 4), the kinematic 
wave theories are shown to yield consistent results with the supply-demand method 
given proper deflnitions of traffic supply and demand. That is, the supply-demand 
method can be considered as the discrete form of the corresponding continuous kine- 
matic wave theory. It is also interesting to directly prove that the discrete model 
in the supply-demand method converge to the corresponding kinematic wave theory. 
For the inhomogeneous LWR model under jump initial conditions, for example, we 
can find solutions by the supply-demand method and obtain a time series for each 
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cell. Then will these series get closer to the wave solutions predicted in Chapter 2 as 
the length of a time interval keeps decreasing? The answer to this problem will help 
to understand the continuous kinematic wave model of merging traffic, whose discrete 
model is simple and clear (Chapter 3). In addition, the diverging model (Chapter 4) 
is an instantaneous approximation, also in discrete form, whose continuous counter- 
part is another interesting topic. One p robable alternat ive approach is to introduce 



appropriate fundamental diagram as in (jPaganzo 



19971) . 



The kinematic wave models studied in this dissertation are deterministic. How- 
ever, these models can be extended to include randomness in a road network. For 
example, non-recurrent traffic congestion is generally caused by accidents or work 
zones, which are reflected by a change in the number of lanes. Since the number of 
lanes has been included in the models we studied, these stochastic effects could be 
considered in the kinematic waves. As another example, weaving effects on traffic 
dynamics could also be studied in a similar fashion. 

Finally, it is of interest to study kinematic waves that arise from the traffic 
dynamics of other traffic systems, such as those with special lanes, special vehicles, 
and so on. 



8.2.2 Calibration, validation, and enhancement of the MCKW 
simulation model 

The enhancement of MCKW simulation model will be carried out in three aspects. 
First, we can incorporate new kinematic wave theories so that it is capable of simu- 
lating traffic systems with more components, such as different types of vehicles and 
special lanes. To achieve this, traffic has to be categorized into more commodities, 
and more complicated data and program structures are expected. Second, we can 
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improve the computational efficiency of the MCKW simulation for large-scale road 
networks. In this dissertation, the solution methods, either the supply-demand meth- 
ods or the Godunov methods, are of ffist order. One improvem ent can be made 



bv applving higher-order methods in the MCKW model (e.g. 



Colella and Puckett 



Daganzo 



1999a : 



20001 ) ■ which are more efficient. Another approach is to incorpo- 



rating parallel algorithms and more efficient memory management methods. Finally, 
we can design the programming of the simulation model in order to satisfy different 
requirements of different applications. For example, we need very detailed informa- 
tion when studying traffic dynamics and phenomena, but need less information in 
traffic assignment. Thus it will be better to use two different programming structure 
for these two applications. 

The MCKW simulation model, as other simulation models, is also subject to 
calibration and validation with observed traffic dynamics. A conceptual approach is 
to calibrate and validate the simulation model ffist for inhomogeneous links, merges, 
and diverges and then for a whole road network. We can check existing databases 
or collect data by ourselves in the calibration and validation process. In many data 
sets collected by loop detectors, however, volumes (the number of vehicles passing a 
detector) are not conserved over a link. This inconsistency will be fatal for calibrating 
the kinematic wave models, since the fundamental assumption in these models is 
traffic conservation. Thus, it may be necessary to seek advanced measures, such as 
vedio-taping or by satellite sensing techniques, to provide the ideal accuracy in data 
in these studies. 



8.2.3 Applications of the MCKW simulation model 

First, using the MCKW simulation model, we can better understand the formation 
and characteristics of traffic congestion. This will form a foundation to alleviate it 
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by either expansion of infrastructure or traffic control and management. 

Second, the MCKW simulation model can be applied in evaluating strategies in 
infrastructure expansion and traffic operations and management. For example, we 
could develop and evaluate on-ramp metering and arterial signal control algorithms by 
using the MCKW model, in which signals have been incorporated in the computation 
of traffic demand (see Chapter 3). In addition, in simulation-based dynamic traffic 
assignment, the MCKW model can be used as a loading model. 

Finally, the MCKW simulation model can be used to estimate the travel demands 
in a road network with observed traffic conditions at a number of locations and can 
be further applied in regional and transportation plan. 
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